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

CiteSeerX logo

Advanced Search Include Citations
Advanced Search Include Citations | Disambiguate

An immersed interface method for incompressible Navier–Stokes equations, (2003)

by L Lee, R J LeVeque
Venue:SIAM J. Sci. Comput.
Add To MetaCart

Tools

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

The immersed boundary method: a projection approach.

by Kunihiko Taira , Tim Colonius - J. Comput. Phys., , 2007
"... Abstract A new formulation of the immersed boundary method with a structure algebraically identical to the traditional fractional step method is presented for incompressible flow over bodies with prescribed surface motion. Like previous methods, a boundary force is applied at the immersed surface t ..."
Abstract - Cited by 59 (12 self) - Add to MetaCart
Abstract A new formulation of the immersed boundary method with a structure algebraically identical to the traditional fractional step method is presented for incompressible flow over bodies with prescribed surface motion. Like previous methods, a boundary force is applied at the immersed surface to satisfy the no-slip constraint. This extra constraint can be added to the incompressible Navier-Stokes equations by introducing regularization and interpolation operators. The current method gives prominence to the role of the boundary force acting as a Lagrange multiplier to satisfy the no-slip condition. This role is analogous to the effect of pressure on the momentum equation to satisfy the divergence-free constraint. The current immersed boundary method removes slip and non-divergence-free components of the velocity field through a projection. The boundary force is determined implicitly without any constitutive relations allowing the present formulation to use larger CFL numbers compared to some past methods. Symmetry and positive-definiteness of the system are preserved such that the conjugate gradient method can be used to solve for the flow field. Examples show that the current formulation achieves second-order temporal accuracy and better than first-order spatial accuracy in L 2 -norms for oneand two-dimensional test problems. Results from two-dimensional simulations of flows over stationary and moving cylinders are in good agreement with those from previous experimental and numerical studies.
(Show Context)

Citation Context

...onal fractional step algorithm, the entire IBM itself has not been regarded as a fractional step (projection) method, as reported here. We follow the algebraic approach of Perot [30], where the fractional step method is written as a block-LU decomposition. In the next section, we review the traditional fractional step method as it is the fundamental basis for our IBM. In Section 3, we introduce the immersed boundary projection method. This formulation is compared to previous methods in Section 4; namely the original IBM [31], the direct forcing method [28], the immersed interface method (IIM) [23], and the distributed Lagrange multiplier (DLM) method [16]. In Section 5, numerical examples are considered to assess the temporal and spatial accuracy of the current method. Flows over stationary and moving cylinders are simulated and results are compared to previous experimental and numerical studies. Section 6 summarizes the current formulation. Further details on the discretization of the immersed boundary projection method are placed in the Appendix.2. Fractional step method We consider the incompressible Navier–Stokes equationsou ot þ u ru rp þ 1 Re r2u; ð1Þ r u 0; ð2Þwhere u, ...

An adaptive, formally second order accurate version of the immersed boundary method

by Boyce E. Griffith , Richard D. Hornung , David M. McQueen , Charles S. Peskin , 2006
"... ..."
Abstract - Cited by 48 (10 self) - Add to MetaCart
Abstract not found

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

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

...f the interface (i.e. we assume that n × τ = k ≡ constant, so that n can point either towards, or outwards from, the center of curvature). The proof for expressions (8)-(11) can be found in detail in =-=[23, 24, 25, 21]-=-. We note that from expressions (8)-(11) the values of the jumps of the first and second derivatives of velocity and pressure with respect to the (x, y) coordinates are easily obtained by a simple coo...

Moving Overlapping Grids with Adaptive Mesh Refinement for High-Speed Reactive and Non-reactive Flow

by William D. Henshaw, Donald W. Schwendeman - Journal of Computational Physics , 1979
"... We consider the solution of the reactive and non-reactive Euler equations on two-dimensional domains that evolve in time. The domains are discretized using moving overlapping grids. In a typical grid construction, boundary-tted grids are used to represent moving boundaries, and these grids overlap w ..."
Abstract - Cited by 23 (6 self) - Add to MetaCart
We consider the solution of the reactive and non-reactive Euler equations on two-dimensional domains that evolve in time. The domains are discretized using moving overlapping grids. In a typical grid construction, boundary-tted grids are used to represent moving boundaries, and these grids overlap with stationary background Cartesian grids. Block-structured adaptive mesh renement (AMR) is used to re-solve ne-scale features in the ow such as shocks and detonations. Renement grids are added to base-level grids according to an estimate of the error, and these re nement grids move with their corresponding base-level grids. The numerical ap-proximation of the governing equations takes place in the parameter space of each component grid which is dened by a mapping from ( xed) parameter space to (moving) physical space. The mapped equations are solved numerically using a second-order extension of Godunov's method. The sti source term in the reac-tive case is handled using a Runge-Kutta error-control scheme. We consider cases when the boundaries move according to a prescribed function of time and when
(Show Context)

Citation Context

... class. This method applies boundary conditions by introducing a body force into the equations [30] and thus eectively smearing the interface. In the immersed interface method and related approaches =-=[31,32]-=-, the interface is kept sharp with special discretization stencils applied where the grid meets the boundary. Fictitious domain methods [33,34] impose boundary conditions through the use of constraint...

Logically Rectangular Grids and Finite Volume Methods for PDEs in Circular and Spherical Domains

by Donna A. Calhoun, Christiane Helzel, Randall, J. Leveque - In preparation; http://www.amath. washington.edu/~rjl/pubs/circles , 2005
"... Abstract. We describe a class of logically rectangular quadrilateral and hexahedral grids for solving PDEs in circular and spherical domains, including grid mappings for the circle, the surface of the sphere and the three-dimensional ball. The grids are logically rectangular and the computational do ..."
Abstract - Cited by 22 (6 self) - Add to MetaCart
Abstract. We describe a class of logically rectangular quadrilateral and hexahedral grids for solving PDEs in circular and spherical domains, including grid mappings for the circle, the surface of the sphere and the three-dimensional ball. The grids are logically rectangular and the computational domain is a single Cartesian grid. Compared to alternative approaches based on a multiblock data structure or unstructured triangulations, this approach simplifies the implementation of numerical methods and the use of adaptive refinement. A more general domain with a smooth boundary can be gridded by composing one of the mappings from this paper with another smooth mapping from the circle or sphere to the desired domain. Although these grids are highly non-orthogonal, we show that the high-resolution wave-propagation algorithm implemented in clawpack can be effectively used to approximate hyperbolic problems on these grids. Since the ratio between the largest and smallest grid is below 2 for most of our grid mappings, explicit finite volume methods such as the wave propagation algorithm do not suffer from the center or pole singularities that arise with polar or latitude-longitude grids. Numerical test calculations illustrate the potential use of these grids for a variety of applications including Euler equations, shallow water equations, and acoustics in a heterogeneous medium. Pattern formation from a reaction-diffusion equation on the sphere is also considered. All examples are implemented in the clawpack software package and full source code is available on the web, along with matlab routines for the various mappings.
(Show Context)

Citation Context

...dvanced by many other researchers as well, for hyperbolic equations (e.g., [16], [17], [19], [20]), elliptic equations (e.g., [27], [38], [39]), and incompressible Navier-Stokes equations (e.g., [2], =-=[32]-=-). In this paper, we focus our attention on simple domains in two space dimensions where the use of a mapped coordinate system is natural and often easier to work with computationally than embedded bo...

A spectral boundary integral method for flowing blood cells

by Hong Zhao , Amir H G Isfahani , Luke N Olson , Jonathan B Freund - J. Comp. Phys., 229:3726
"... a b s t r a c t A spectral boundary integral method for simulating large numbers of blood cells flowing in complex geometries is developed and demonstrated. The blood cells are modeled as finitedeformation elastic membranes containing a higher viscosity fluid than the surrounding plasma, but the so ..."
Abstract - Cited by 22 (1 self) - Add to MetaCart
a b s t r a c t A spectral boundary integral method for simulating large numbers of blood cells flowing in complex geometries is developed and demonstrated. The blood cells are modeled as finitedeformation elastic membranes containing a higher viscosity fluid than the surrounding plasma, but the solver itself is independent of the particular constitutive model employed for the cell membranes. The surface integrals developed for solving the viscous flow, and thereby the motion of the massless membrane, are evaluated using an OðN log NÞ particle-mesh Ewald (PME) approach. The cell shapes, which can become highly distorted under physiologic conditions, are discretized with spherical harmonics. The resolution of these global basis functions is, of course, excellent, but more importantly they facilitate an approximate de-aliasing procedure that stabilizes the simulations without adding any numerical dissipation or further restricting the permissible numerical time step. Complex geometry no-slip boundaries are included using a constraint method that is coupled into an implicit system that is solved as part of the time advancement routine. The implementation is verified against solutions for axisymmetric flows reported in the literature, and its accuracy is demonstrated by comparison against exact solutions for relaxing surface deformations. It is also used to simulate flow of blood cells at 30% volume fraction in tubes between 4.9 and 16.9 lm in diameter. For these, it is shown to reproduce the well-known non-monotonic dependence of the effective viscosity on the tube diameter.
(Show Context)

Citation Context

...heir close-range interactions. For simulation methods requiring boundary-fitted meshes, such as arbitrary Lagrangian–Eulerian methods, the mesh movement and regeneration for the complicated flow geometry is challenging and computationally expensive. Immersed boundary or immersed interface methods can avoid this by representing the changing geometry on a fixed mesh. These methods account for the cell–fluid interactions by diffusing the membrane force to the nearby fluid domain mesh points [28] or by modifying the local discrete differential operators when their stencils cross the cell membrane [21]. Both cell–cell interactions and cellular flow in small blood vessels have been simulated by the immersed boundary method in two dimensions [4,3]. The immersed finite-element method [23] and the lattice-Boltzmann method [24] similarly model the cell membranes as added body forces and have been used for studying cell aggregation and blood rheology. These fixed-mesh methods avoid re-meshing, but still require either a fine volumetric mesh or adaptive meshing to resolve close interactions between cells [39]. The mesh-based Lagrangian–Eulerian and immersed boundary or interface methods start with...

Penalty immersed boundary method for an elastic boundary with mass

by Yongsam Kim, Charles S. Peskin - Physics of Fluids
"... The immersed boundary (IB) method has been widely applied to problems involving a moving elastic boundary that is immersed in fluid and interacting with it. But most of the previous applications of the IB method have involved a massless elastic boundary and used efficient Fourier transform methods f ..."
Abstract - Cited by 22 (3 self) - Add to MetaCart
The immersed boundary (IB) method has been widely applied to problems involving a moving elastic boundary that is immersed in fluid and interacting with it. But most of the previous applications of the IB method have involved a massless elastic boundary and used efficient Fourier transform methods for the numerical solutions. Extending the method to cover the case of a massive boundary has required spreading the boundary mass out onto the fluid grid and then solving the Navier-Stokes equations with a variable mass density. The variable mass density of this previous approach makes Fourier transform methods inapplicable, and requires a multigrid solver. Here we propose a new and simple way to give mass to the elastic boundary and show that the new method can be applied to many problems for which the boundary mass is important. The new method does not spread mass to the fluid grid, retains the use of Fourier transform methodology, and is easy to implement in the context of an existing IB method code for the massless case.
(Show Context)

Citation Context

...conserva28tion” [25], which is not implemented here, since most of our applications do not involve enclosed volumes. For other approaches that also lead to better conservation of volume, see [4] and =-=[16]-=-. 5 Applications 5.1 Flapping Filament in a Flowing Soap Film The first application is a flapping flexible filament in a flowing soap film. This replicates a simulation done by Luoding Zhu [27,28] by ...

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.
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