scholarly journals An Unconditionally Energy Stable Immersed Boundary Method with Application to Vesicle Dynamics

2013 ◽  
Vol 3 (3) ◽  
pp. 247-262 ◽  
Author(s):  
Wei-Fan Hu ◽  
Ming-Chih Lai
Keyword(s):  
Immersed Boundary Method ◽  
Navier Stokes ◽  
Immersed Boundary ◽  
Step Size ◽  
Time Step ◽  
Time Step Size ◽  
Boundary Method ◽  
Stretching Factor ◽  
Vesicle Dynamics ◽  
Energy Stable

AbstractWe develop an unconditionally energy stable immersed boundary method, and apply it to simulate 2D vesicle dynamics. We adopt a semi-implicit boundary forcing approach, where the stretching factor used in the forcing term can be computed from the derived evolutional equation. By using the projection method to solve the fluid equations, the pressure is decoupled and we have a symmetric positive definite system that can be solved efficiently. The method can be shown to be unconditionally stable, in the sense that the total energy is decreasing. A resulting modification benefits from this improved numerical stability, as the time step size can be significantly increased (the severe time step restriction in an explicit boundary forcing scheme is avoided). As an application, we use our scheme to simulate vesicle dynamics in Navier-Stokes flow.

scholarly journals Immersed boundary method combined with proper generalized decomposition for simulation of a flexible filament in a viscous incompressible flow

2017 ◽  
Vol 39 (2) ◽  
pp. 109-119
Author(s):  
Cuong Q. Le ◽  
H. Phan-Duc ◽  
Son H. Nguyen
Keyword(s):  
Immersed Boundary Method ◽  
Stokes Equations ◽  
Fluid Pressure ◽  
Navier Stokes ◽  
Immersed Boundary ◽  
Proper Generalized Decomposition ◽  
Time Step ◽  
Navier Stokes Equations ◽  
Boundary Method ◽  
Fluid Filament

In this paper, a combination of the Proper Generalized  Decomposition (PGD) with the Immersed Boundary method (IBM) for solving  fluid-filament interaction problem is proposed. In this combination, a  forcing term constructed by the IBM is introduced to Navier-Stokes equations  to handle the influence of the filament on the fluid flow. The PGD is  applied to solve the Poission's equation to find the fluid pressure  distribution for each time step. The numerical results are compared with  those by previous publications to illustrate the robustness and  effectiveness of the proposed method.

scholarly journals Analysis of a Decoupled Time-Stepping Scheme for Evolutionary Micropolar Fluid Flows

2016 ◽  
Vol 2016 ◽  
pp. 1-13 ◽  
Author(s):  
S. S. Ravindran
Keyword(s):  
Micropolar Fluid ◽  
Stokes Equations ◽  
Fluid Model ◽  
Navier Stokes ◽  
Optimal Order ◽  
Step Size ◽  
Time Step ◽  
Order Error ◽  
Time Step Size ◽  
Time Stepping

Micropolar fluid model consists of Navier-Stokes equations and microrotational velocity equations describing the dynamics of flows in which microstructure of fluid is important. In this paper, we propose and analyze a decoupled time-stepping algorithm for the evolutionary micropolar flow. The proposed method requires solving only one uncoupled Navier-Stokes and one microrotation subphysics problem per time step. We derive optimal order error estimates in suitable norms without assuming any stability condition or time step size restriction.

scholarly journals A Stress Mapping Immersed Boundary Method for Viscous Flows

2021 ◽  
Vol 87 (3) ◽  
pp. 1-20
Author(s):  
Karim M. Ali ◽  
Mohamed Madbouli ◽  
Hany M. Hamouda ◽  
Amr Guaily
Keyword(s):  
Immersed Boundary Method ◽  
Stokes Equations ◽  
Cross Flow ◽  
Navier Stokes ◽  
Immersed Boundary ◽  
Element Formulation ◽  
Navier Stokes Equations ◽  
Boundary Method ◽  
Background Grid ◽  
Dimensional Simulation

This work introduces an immersed boundary method for two-dimensional simulation of incompressible Navier-Stokes equations. The method uses flow field mapping on the immersed boundary and performs a contour integration to calculate immersed boundary forces. This takes into account the relative location of the immersed boundary inside the background grid elements by using inverse distance weights, and also considers the curvature of the immersed boundary edges. The governing equations of the fluid mechanics are solved using a Galerkin-Least squares finite element formulation. The model is validated against a stationary and a vertically oscillating circular cylinder in a cross flow. The results of the model show acceptable accuracy when compared to experimental and numerical results.

An Immersed Boundary Method for Compressible Flows with Complex Boundaries

2013 ◽  
Vol 477-478 ◽  
pp. 281-284
Author(s):  
Jie Yang ◽  
Song Ping Wu
Keyword(s):  
Immersed Boundary Method ◽  
Compressible Flows ◽  
Stokes Equations ◽  
Splitting Method ◽  
Linear Interpolation ◽  
High Order Accuracy ◽  
Navier Stokes ◽  
Immersed Boundary ◽  
Central Difference ◽  
Boundary Method

An immersed boundary method based on the ghost-cell approach is presented in this paper. The compressible Navier-Stokes equations are discretized using a flux-splitting method for inviscid fluxes and second-order central-difference for the viscous components. High-order accuracy is achieved by using weighted essentially non-oscillatory (WENO) and Runge-Kutta schemes. Boundary conditions are reconstructed by a serial of linear interpolation and inverse distance weighting interpolation of flow variables in fluid domain. Two classic flow problems (flow over a circular cylinder, and a NACA 0012 airfoil) are simulated using the present immersed boundary method, and the predictions show good agreement with previous computational results.

scholarly journals Stability of Finite Difference Discretizations of Multi-Physics Interface Conditions

2013 ◽  
Vol 13 (2) ◽  
pp. 386-410 ◽  
Author(s):  
Björn Sjögreen ◽  
Jeffrey W. Banks
Keyword(s):  
Stokes Equations ◽  
Normal Mode Analysis ◽  
Linear Equations ◽  
Mode Analysis ◽  
Navier Stokes ◽  
Computational Domain ◽  
Interface Conditions ◽  
Step Size ◽  
Time Step ◽  
Time Step Size

AbstractWe consider multi-physics computations where the Navier-Stokes equations of compressible fluid flow on some parts of the computational domain are coupled to the equations of elasticity on other parts of the computational domain. The different subdomains are separated by well-defined interfaces. We consider time accurate computations resolving all time scales. For such computations, explicit time stepping is very efficient. We address the issue of discrete interface conditions between the two domains of different physics that do not lead to instability, or to a significant reduction of the stable time step size. Finding such interface conditions is non-trivial.We discretize the problem with high order centered difference approximations with summation by parts boundary closure. We derive L2 stable interface conditions for the linearized one dimensional discretized problem. Furthermore, we generalize the interface conditions to the full non-linear equations and numerically demonstrate their stable and accurate performance on a simple model problem. The energy stable interface conditions derived here through symmetrization of the equations contain the interface conditions derived through normal mode analysis by Banks and Sjögreen in [8] as a special case.

Sign in / Sign up

Export Citation Format

Share Document