scispace - formally typeset
Open AccessJournal ArticleDOI

Matched interface and boundary (MIB) method for elliptic problems with sharp-edged interfaces

Sining Yu, +2 more
- 01 Jun 2007 - 
- Vol. 224, Iss: 2, pp 729-756
Reads0
Chats0
TLDR
This work generalizes the matched interface and boundary method previously designed for solving elliptic problems with curved interfaces to the aforementioned problems, and introduces primary and secondary fictitious values to deal with sharp-edged interfaces.
About
This article is published in Journal of Computational Physics.The article was published on 2007-06-01 and is currently open access. It has received 181 citations till now. The article focuses on the topics: Singularity & Boundary (topology).

read more

Citations
More filters
Journal ArticleDOI

Biomolecular electrostatics and solvation: a computational perspective.

TL;DR: The solvation of biomolecules with a computational biophysics view toward describing the phenomenon is discussed, and the main focus lies on the computational aspect of the models.
Journal ArticleDOI

Treatment of charge singularities in implicit solvent models

TL;DR: This paper presents a novel method for solving the Poisson-Boltzmann (PB) equation based on a rigorous treatment of geometric singularities of the dielectric interface and a Green's function formulation of charge singularities, which provides reliable benchmark numerical solutions for the PB equation.
Journal ArticleDOI

MIBPB: A software package for electrostatic analysis

TL;DR: A matched interface and boundary (MIB)‐based PBE software package, the MIBPB solver, for electrostatic analysis, and further accelerates the rate of convergence of linear equation systems resulting from the M IBPB by using the Krylov subspace (KS) techniques.
Journal ArticleDOI

Poisson-Boltzmann-Nernst-Planck model.

TL;DR: There is an excellent consistency between the results predicted from the present PBNP model and those obtained from the PNP model in terms of the electrostatic potentials, ion concentration profiles, and current-voltage curves.
Journal ArticleDOI

Differential geometry based solvation model I: Eulerian formulation

TL;DR: This paper presents a differential geometry based model for the analysis and computation of the equilibrium property of solvation, and derives coupled generalized Poisson-Boltzmann equation and generalized geometric flow equation for the electrostatic potential and the construction of realistic solvent-solute boundaries.
References
More filters
Journal ArticleDOI

Immersed boundary methods

TL;DR: The term immersed boundary (IB) method is used to encompass all such methods that simulate viscous flows with immersed (or embedded) boundaries on grids that do not conform to the shape of these boundaries.
Journal ArticleDOI

Numerical analysis of blood flow in the heart

TL;DR: In this article, the authors extended previous work on the solution of the Navier-Stokes equations in the presence of moving immersed boundaries which interact with the fluid and introduced an improved numerical representation of the δ-function.
Journal ArticleDOI

A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method)

TL;DR: A new numerical method for treating interfaces in Eulerian schemes that maintains a Heaviside profile of the density with no numerical smearing along the lines of earlier work and most Lagrangian schemes is proposed.
Journal ArticleDOI

Combined Immersed-Boundary Finite-Difference Methods for Three-Dimensional Complex Flow Simulations

TL;DR: In this paper, a second-order accurate, highly efficient method is developed for simulating unsteady three-dimensional incompressible flows in complex geometries, which is achieved by using boundary body forces that allow the imposition of the boundary conditions on a given surface not coinciding with the computational grid.
Journal ArticleDOI

the immersed interface method for elliptic equations with discontinuous coefficients and singular sources

TL;DR: In this paper, the authors developed finite difference methods for elliptic equations of the form \[ abla \cdot (\beta (x)) + \kappa (x)u(x) = f(x)) in a region in one or two dimensions.
Related Papers (5)
Frequently Asked Questions (15)
Q1. What are the contributions mentioned in the paper "Matched interface and boundary (mib) method for elliptic problems with sharp-edged interfaces" ?

This work generalizes the matched interface and boundary ( MIB ) method previously designed for solving elliptic problems with curved interfaces to the aforementioned problems. The concepts of primary and secondary fictitious values are introduced to deal with sharp-edged interfaces. 

The present work provides a solution to this class of problems on the Cartesian grid by extending the marched interface and boundary ( MIB ) method [ 60,61 ] previously designed for straight or curved interfaces. Second, the standard ( higherorder ) central finite difference ( FD ) schemes are utilized so that the condition number of the discretization matrix is relatively low, and efficient linear algebraic equation solvers can be used. 

The concept of secondary fictitious values is essential for overcoming the geometric difficulties created by sharp-edged interfaces. 

One of the most challenging problems in the field is the solution of elliptic equations with sharp-edged coefficients, i.e., non-smooth interfaces. 

One feature of the MIB is that it disassociates between the discretization of the elliptic equation and the enforcement of interface jump conditions. 

For irregular interfaces, it is nature to construct a solution in the finite element method formulation [2,9,38], in particular, using the discontinuous Galerkin technique [22]. 

dramatic local mesh refinement is required in the vicinity of sharp edges [13], and leads to severe increase in computational time and memory requirement. 

The standard MIB method can always handle a C1 irregular interface even if it has large curvature, since enough auxiliary points can always be found by refining mesh. 

For the problems with smooth (say C1) interfaces, the present method can easily achieve 2nd order convergence since required auxiliary points can always be found by refining the mesh. 

To restore the order of convergence of the discretization at irregular points, the authors need to treat (bux)x and (buy)y to the designed order of convergence near the interface. 

If only one of F component, Fi (i 2 {1,2,3,4}), is the fictitious value whose expansion is to be determined and the other three are function values or have expansions in term of function values, then the jump condition of u(x,y) is enough for the determination of Fi:uþ u ¼ ½u : ð24ÞLet wþ0 ¼ ðwþ0;i 1;wþ0;i;wþ0;iþ1Þ and w 0 ¼ ðw 0;i;w 0;iþ1;w 0;iþ2Þ. 

In the present work, the authors further generalize the MIB method to allow the presence of sharp-edged interfaces, thin-layered interfaces and interfaces that intersect with the domain boundary. 

The MIB approach makes use of fictitious domains so that the standard high order central finite difference (FD) method can be applied across the interface without the loss of accuracy. 

If (three surrounding points are on the same side or have secondary fictitious value on the same side), then * call On-interface scheme 3 – else pause If ((i, j) is an irregular point and (i, j) 2 X ), use primary fictitious values and call the standard five-point FD scheme. 

This difficulty is not resolvable by the previous MIB scheme because the auxiliary points in the ‘improper’ region leads to inconsistence between the number of unknowns and the number of jump conditions.