• Documents
  • Authors
  • Tables
  • Log in
  • Sign up
  • MetaCart
  • DMCA
  • Donate

CiteSeerX logo

Advanced Search Include Citations
Advanced Search Include Citations | Disambiguate

The immersed interface method for the Navier-Stokes equations with singular forces. (2001)

by Z Li, Lai MC
Venue:Journal of Computational Physics
Add To MetaCart

Tools

Sorted by:
Results 1 - 10 of 83
Next 10 →

An immersed interface method for incompressible navier-stokes equations

by Long Lee, Randall, J. Leveque - SIAM J. Sci. Comput , 2003
"... Abstract. The method developed in this paper is motivated by Peskin’s immersed boundary (IB) method, and allows one to model the motion of flexible membranes or other structures immersed in viscous incompressible fluid using a fluid solver on a fixed Cartesian grid. The IB method uses a set of discr ..."
Abstract - Cited by 99 (3 self) - Add to MetaCart
Abstract. The method developed in this paper is motivated by Peskin’s immersed boundary (IB) method, and allows one to model the motion of flexible membranes or other structures immersed in viscous incompressible fluid using a fluid solver on a fixed Cartesian grid. The IB method uses a set of discrete delta functions to spread the entire singular force exerted by the immersed boundary to the nearby fluid grid points. Our method instead incorporates part of this force into jump conditions for the pressure, avoiding discrete dipole terms that adversely affect the accuracy near the immersed boundary. This has been implemented for the two-dimensional incompressible Navier–Stokes equations using a high-resolution finite-volume method for the advective terms and a projection method to enforce incompressibility. In the projection step, the correct jump in pressure is imposed in the course of solving the Poisson problem. This gives sharp resolution of the pressure across the interface and also gives better volume conservation than the traditional IB method. Comparisons between this method and the IB method are presented for several test problems. Numerical studies of the convergence and order of accuracy are included.

An immersed interface method for simulating the interaction of a . . .

by Sheng Xu, Z. Jane Wang , 2006
"... ..."
Abstract - Cited by 45 (4 self) - Add to MetaCart
Abstract not found

Unconditionally stable discretizations of the immersed boundary equations

by Elijah P. Newren , Aaron L. Fogelson , Robert D. Guy , Robert M. Kirby , 2007
"... ..."
Abstract - Cited by 41 (8 self) - Add to MetaCart
Abstract not found

A fast solver for the Stokes equations with distributed forces in complex geometries

by George Biros, Lexing Ying, Denis Zorin - J. Comput. Phys
"... We present a new method for the solution of the Stokes equations. The main features of our method are: (1) it can be applied to arbitrary geometries in a black-box fashion; (2) it is second order accurate; and (3) it has optimal algorithmic complexity. Our approach, to which we refer as the Embedded ..."
Abstract - Cited by 41 (10 self) - Add to MetaCart
We present a new method for the solution of the Stokes equations. The main features of our method are: (1) it can be applied to arbitrary geometries in a black-box fashion; (2) it is second order accurate; and (3) it has optimal algorithmic complexity. Our approach, to which we refer as the Embedded Boundary Integral method, is based on Anita Mayo’s work for the Poisson’s equation: “The Fast Solution of Poisson’s and the Biharmonic Equations on Irregular Regions”, SIAM Journal on Numerical Analysis, 21 (1984), pp. 285–299. We embed the domain in a rectangular domain, for which fast solvers are available, and we impose the boundary conditions as interface (jump) conditions on the velocities and tractions. We use an indirect boundary integral formulation for the homogeneous Stokes equations to compute the jumps. The resulting equations are discretized by Nyström’s method. The rectangular domain problem is discretized by finite elements for a velocitypressure formulation with equal order interpolation bilinear elements (£¥ ¤-£¥ ¤). Stabilization is used to circumvent the ¦¨§�©������� � condition for the pressure space. For the integral equations, fast matrix vector multiplications are achieved via an ���¨���� � algorithm based on a block representation of the discrete integral operator, combined with (kernel independent) singular value decomposition to sparsify low-rank blocks. The regular grid solver is a Krylov method (Conjugate Residuals) combined with an optimal two-level Schwartz-preconditioner. For the integral equation we use GMRES. We have tested our algorithm on several numerical examples and we have observed optimal convergence rates. Key Words: Stokes equations, fast solvers, integral equations, double-layer potential, fast multipole methods, embedded domain methods, immersed interface methods, fictitious
(Show Context)

Citation Context

...nuous coefficients and singular forces. It has been successfully applied to moving boundary problems, for example for the Stokes problem with elastic interfaces [32] and for the Navier-Stokes problem =-=[35]-=-. If the singular forces are known then the jumps are known and TESCs can be computed explicitly. For discontinuous coefficients IIM modifies the stencils for points close to the boundary in order to ...

An immersed interface method for viscous incompressible flows involving rigid and flexible boundaries

by D. V. Le, B. C. Khoo, J. Peraire - J. Comp. Phys , 2006
"... We present an immersed interface method for the incompressible Navier-Stokes equations capable of handling rigid immersed boundaries. The immersed boundary is represented by a set of Lagrangian control points. In order to guarantee that the no-slip condition on the boundary is satisfied, singular fo ..."
Abstract - Cited by 37 (3 self) - Add to MetaCart
We present an immersed interface method for the incompressible Navier-Stokes equations capable of handling rigid immersed boundaries. The immersed boundary is represented by a set of Lagrangian control points. In order to guarantee that the no-slip condition on the boundary is satisfied, singular forces are applied on the fluid. The forces are related to the jumps in pressure and the jumps in the derivatives of both pressure and velocity, and are interpolated using cubic splines. The strength of the singular forces is determined by solving a small system of equations iteratively at each time step. The Navier-Stokes equations are discretized on a staggered Cartesian grid by a second order accurate projection method for pressure and velocity. Keywords: Immersed interface method, Navier-Stokes equations, Cartesian grid method, finite difference, fast Poisson solvers, irregular domains.
(Show Context)

Citation Context

...i [23] for solving elliptic equations, and later extended to Stokes flow with elastic boundaries or surface tension [24]. The method was developed further for the Navier-Stokes equations in Li et al. =-=[25]-=-, Lee [22] and Le et al. [19] for problems with flexible boundaries. The method was also used by Calhoun [6] and Li et al. [26] for solving the two-dimensional streamfunction-vorticity equations on ir...

A Cartesian grid method for solving the two-dimensional streamfunction-vorticity equations in irregular regions

by Donna Calhoun - J. Comput. Phys
"... We describe a method for solving the two-dimensional Navier–Stokes equations in irregular physical domains. Our method is based on an underlying uniform Cartesian grid and second-order finite-difference/finite-volume discretizations of the streamfunction-vorticity equations. Geometry representing st ..."
Abstract - Cited by 32 (2 self) - Add to MetaCart
We describe a method for solving the two-dimensional Navier–Stokes equations in irregular physical domains. Our method is based on an underlying uniform Cartesian grid and second-order finite-difference/finite-volume discretizations of the streamfunction-vorticity equations. Geometry representing stationary solid obstacles in the flow domain is embedded in the Cartesian grid and special discretizations near the embedded boundary ensure the accuracy of the solution in the cut cells. Along the embedded boundary, we determine a distribution of vorticity sources needed to impose the no-slip flow conditions. This distribution appears as a right-hand-side term in the discretized fluid equations, and so we can use fast solvers to solve the linear systems that arise. To handle the advective terms, we use the high-resolution algorithms in CLAWPACK. We show that our Stokes solver is second-order accurate for steady state solutions and that our full Navier–Stokes solver is between first- and second-order accurate and reproduces results from well-studied benchmark problems in viscous fluid flow. Finally, we demonstrate the robustness of our code on flow in
(Show Context)

Citation Context

...apply our finite-difference scheme to any of these methods. For one example of how the immersed interface method can be applied to a primitive variable formulation of the Navier–Stokes equations, see =-=[29]-=-. APPENDIX A.1. Derivation of Coefficients for the Immersed Interface Method Here we provide formulas for the coefficients needed to construct the matrix S in (33) and for extending the solution ω to ...

Numerical Methods for Fluid-Structure Interaction -- A Review

by Gene Hou, Jin Wang, Anita Layton , 2012
"... The interactions between incompressible fluid flows and immersed struc-tures are nonlinear multi-physics phenomena that have applications to a wide range of scientific and engineering disciplines. In this article, we review representative numerical methods based on conforming and non-conforming me ..."
Abstract - Cited by 21 (0 self) - Add to MetaCart
The interactions between incompressible fluid flows and immersed struc-tures are nonlinear multi-physics phenomena that have applications to a wide range of scientific and engineering disciplines. In this article, we review representative numerical methods based on conforming and non-conforming meshes that are currently avail-able for computing fluid-structure interaction problems, with an emphasis on some of the recent developments in the field. A goal is to categorize the selected methods and assess their accuracy and efficiency. We discuss challenges faced by researchers in this field, and we emphasize the importance of interdisciplinary effort for advancing the study in fluid-structure interactions.

Numerical treatment of two-dimensional interfaces for acoustic and elastic waves

by Bruno Lombard - J. Comput. Phys
"... We present a numerical method to take into account 2D arbitrary-shaped inter-faces in classical finite-difference schemes, on a uniform Cartesian grid. This work extends the “Explicit Simplified Interface Method ” (ESIM), previously proposed in 1D (2001, J. Comput. Phys. 168, pp. 227-248). The physi ..."
Abstract - Cited by 19 (10 self) - Add to MetaCart
We present a numerical method to take into account 2D arbitrary-shaped inter-faces in classical finite-difference schemes, on a uniform Cartesian grid. This work extends the “Explicit Simplified Interface Method ” (ESIM), previously proposed in 1D (2001, J. Comput. Phys. 168, pp. 227-248). The physical problem under study concerns the linear hyperbolic systems of acoustics and elastodynamics, with stationary interfaces. Our method maintains, near the interfaces, properties of the schemes in homogeneous medium, such as the order of accuracy and the stability limit. Moreover, it enforces the numerical solution to satisfy the exact interface conditions. Lastly, it provides subcell geometrical features of the interface inside the meshing. The ESIM can be coupled automatically with a wide class of nu-merical schemes (Lax-Wendroff, flux-limiter schemes,...) for a negligible additional computational cost. Throughout the paper, we focus on the challenging case of an interface between a fluid and an elastic solid. In numerical experiments, we provide comparisons between numerical solutions and original analytic solutions, showing the efficiency of the method.
(Show Context)

Citation Context

... the advection or the Maxwell equations [4] with interfaces. 29 The extension of the ESIM to nonlinear problems is a challenging project. We think to CFD problems, such as the Navier-Stokes equations =-=[15]-=- or the Euler equations with moving interfaces. We have especially in mind the numerical treatment of material interfaces in multicomponent flows. The ”Ghost Fluid Method” of Fedkiw and al. [5] is eff...

A level-set method for interfacial flows with surfactant

by Jian-jun Xu , Zhilin Li , John Lowengrub , Hongkai Zhao , 2006
"... ... drop deformations and more complex drop–drop interactions compared to the analogous cases for clean drops. The effects of surfactant are found to be most significant in flows with multiple drops. To our knowledge, this is the first time that the level-set method has been used to simulate fluid i ..."
Abstract - Cited by 17 (1 self) - Add to MetaCart
... drop deformations and more complex drop–drop interactions compared to the analogous cases for clean drops. The effects of surfactant are found to be most significant in flows with multiple drops. To our knowledge, this is the first time that the level-set method has been used to simulate fluid interfaces with surfactant.
(Show Context)

Citation Context

...s described briefly in section 3.1. The extension of this approach to flows with non-zero Reynolds numbers (but matched viscosities and densities) is straightforward using the techniques described in =-=[27, 37, 30]-=- for the IIM. The extension to flows with variable viscosity has also been recently developed [36]. Finally, the nondimensional form of the surfactant equation is identical to Eq. (9) with Ds replaced...

A new high-order immersed interface method for solving elliptic equations with imbedded interface of discontinuity

by Xiaolin Zhong - J. Comput. Phys , 2007
"... Abstract This paper presents a new high-order immersed interface method for elliptic equations with imbedded interface of discontinuity. Compared with the original second-order immersed interface method of [R.J. LeVeque, Z. Li. The immersed interface method for elliptic equations with discontinuous ..."
Abstract - Cited by 13 (1 self) - Add to MetaCart
Abstract This paper presents a new high-order immersed interface method for elliptic equations with imbedded interface of discontinuity. Compared with the original second-order immersed interface method of [R.J. LeVeque, Z. Li. The immersed interface method for elliptic equations with discontinuous coefficients and singular sources. SIAM J. Numer. Anal. 31 (1994) 1001-25], the new method achieves arbitrarily high-order accuracy for derivatives at an irregular grid point by imposing only two physical jump conditions together with a wider set of grid stencils. The new interface difference formulas are expressed in a general explicit form so that they can be applied to different multi-dimensional problems without any modification. The new interface algorithms of up to O(h 4 ) accuracy have been derived and tested on several one and twodimensional elliptic equations with imbedded interface. Compared to the standard second-order immersed interface method, the test results show that the new fourth-order immersed interface method leads to a significant improvement in accuracy of the numerical solutions. The proposed method has potential advantages in the application to two-phase flow because of its high-order accuracy and simplicity in applications.
(Show Context)

Citation Context

...,1,8,2]. Among these developments, Wiegmann and Bube [55,56] developed an explicit-jump immersed interface method for the special cases where the explicit jump conditions of variables and derivatives ([u], [ux], [uxx], etc.) are known. Though in a simpler explicit form, the explicit-jump immersed interface method is not applicable to the general jump conditions given by (2). The immersed interface methods have been applied to the Stokes flow with elastic boundaries or surface tension [25], Hele–Shaw flow [17], incompressible flow based on the Navier–Stokes equations with singular source terms [29,5,23], and nonlinear problems in magneto-rheological fluids [19]. Despite these applications, the immersed interface methods are often difficult to apply to complex two or three-dimensional two-phase flow problems. In order to maintain a second-order accuracy, it is necessary to obtain jump conditions at the interface for flow variables and their first and second derivatives. For the Navier–Stokes equations with an interface of discontinuity, it is easy to derive the physical jump conditions for flow variables and their first derivatives across the interface. But it is difficult to obtain jump cond...

Powered by: Apache Solr
  • About CiteSeerX
  • Submit and Index Documents
  • Privacy Policy
  • Help
  • Data
  • Source
  • Contact Us

Developed at and hosted by The College of Information Sciences and Technology

© 2007-2019 The Pennsylvania State University