skip to main content

This content will become publicly available on March 10, 2023

Title: Passive concentration dynamics incorporated into the library IB2d, a two-dimensional implementation of the immersed boundary method
Abstract In this paper, we present an open-source software library that can be used to numerically simulate the advection and diffusion of a chemical concentration or heat density in a viscous fluid where a moving, elastic boundary drives the fluid and acts as a source or sink. The fully-coupled fluid-structure interaction problem of an elastic boundary in a viscous fluid is solved using Peskin’s immersed boundary method. The addition or removal of the concentration or heat density from the boundary is solved using an immersed boundary-like approach in which the concentration is spread from the immersed boundary to the fluid using a regularized delta function. The concentration or density over time is then described by the advection-diffusion equation and numerically solved. This functionality has been added to our software library, IB2d , which provides an easy-to-use immersed boundary method in two dimensions with full implementations in MATLAB and Python. We provide four examples that illustrate the usefulness of the method. A simple rubber band that resists stretching and absorbs and releases a chemical concentration is simulated as a first example. Complete convergence results are presented for this benchmark case. Three more biological examples are presented: (1) an oscillating row of more » cylinders, representative of an idealized appendage used for filter-feeding or sniffing, (2) an oscillating plate in a background flow is considered to study the case of heat dissipation in a vibrating leaf, and (3) a simplified model of a pulsing soft coral where carbon dioxide is taken up and oxygen is released as a byproduct from the moving tentacles. This method is applicable to a broad range of problems in the life sciences, including chemical sensing by antennae, heat dissipation in plants and other structures, the advection-diffusion of morphogens during development, filter-feeding by marine organisms, and the release of waste products from organisms in flows. « less
; ; ;
Award ID(s):
2111765 1853608
Publication Date:
Journal Name:
Bioinspiration & Biomimetics
Sponsoring Org:
National Science Foundation
More Like this
  1. Interactions between an evolving solid and inviscid flow can result insubstantial computational complexity, particularly in circumstances involving varied boundary conditions between the solid and fluid phases. Examples of such interactions include melting, sublimation, and deflagration, all of which exhibit bidirectional coupling, mass/heat transfer, and topological change of the solid-fluid interface. The diffuse interface method is a powerful technique that has been used to describe a wide range of solid-phase interface-driven phenomena. The implicit treatment of the interface eliminates the need for cumbersome interface tracking, and advances in adaptive mesh refinement have provided a way to sufficiently resolve diffuse interfaces without excessive computational cost. However, the general scale-invariant coupling of these techniques to flow solvers has been relatively unexplored. In this work, a robust method is presented for treating diffuse solid-fluid interfaces with arbitrary boundary conditions. Source terms defined over the diffuse region mimic boundary conditions at the solid-fluid interface, and it is demonstrated that the diffuse length scale has no adverse effects. To show the efficacy of the method, a one-dimensional implementation is introduced and tested for three types of boundaries: mass flux through the boundary, a moving boundary, and passive interaction of the boundary with an incident acoustic wave.more »These demonstrate expected behavior in all cases. Convergence analysis is also performed and compared against the sharp-interface solution, and linear convergence is observed. This method lays the groundwork for the extension to viscous flow, and the solution of problems involving time-varying mass-flux boundaries.« less
  2. The immersed boundary method is a widely used mixed Eulerian/Lagrangian framework for simulating the motion of elastic structures immersed in viscous fluids. In this work, we consider a poroelastic immersed boundary method in which a fluid permeates a porous, elastic structure of negligible volume fraction, and extend this method to include stress relaxation of the material. The porous viscoelastic method presented here is validated for a prescribed oscillatory shear and for an expansion driven by the motion at the boundary of a circular material by comparing numerical solutions to an analytical solution of the Maxwell model for viscoelasticity. Finally, an application of the modelling framework to cell biology is provided: passage of a cell through a microfluidic channel. We demonstrate that the rheology of the cell cytoplasm is important for capturing the transit time through a narrow channel in the presence of a pressure drop in the extracellular fluid.
  3. Nematocysts are stinging organelles used by members of the phylum Cnidaria (e.g., jellyfish, anemones, hydrozoans) for a variety of important functions including capturing prey and defense. Nematocysts are the fastest-known accelerating structures in the animal world. The small scale (microns) coupled with rapid acceleration (in excess of 5 million g) present significant challenges in imaging that prevent detailed descriptions of their kinematics. The immersed boundary method was used to numerically simulate the dynamics of a barb-like structure accelerating a short distance across Reynolds numbers ranging from 0.9–900 towards a passive elastic target in two dimensions. Results indicate that acceleration followed by coasting at lower Reynolds numbers is not sufficient for a nematocyst to reach its target. The nematocyst’s barb-like projectile requires high accelerations in order to transition to the inertial regime and overcome the viscous damping effects normally encountered at small cellular scales. The longer the barb is in the inertial regime, the higher the final velocity of the projectile when it touches its target. We find the size of the target prey does not dramatically affect the barb’s approach for large enough values of the Reynolds number, however longer barbs are able to accelerate a larger amount of surroundingmore »fluid, which in turn allows the barb to remain in the inertial regime for a longer period of time. Since the final velocity is proportional to the force available for piercing the membrane of the prey, high accelerations that allow the system to persist in the inertial regime have implications for the nematocyst’s ability to puncture surfaces such as cellular membranes or even crustacean cuticle.« less
  4. In the phase-field description of moving contact line problems, the two-phase system can be described by free energies, and the constitutive relations can be derived based on the assumption of energy dissipation. In this work we propose a novel boundary condition for contact angle hysteresis by exploring wall energy relaxation, which allows the system to be in non-equilibrium at the contact line. Our method captures pinning, advancing and receding automatically without the explicit knowledge of contact line velocity and contact angle. The microscopic dynamic contact angle is computed as part of the solution instead of being imposed. Furthermore, the formulation satisfies a dissipative energy law, where the dissipation terms all have their physical origin. Based on the energy law, we develop an implicit finite element method that is second order in time. The numerical scheme is proven to be unconditionally energy stable for matched density and zero contact angle hysteresis, and is numerically verified to be energy dissipative for a broader range of parameters. We benchmark our method by computing pinned drops and moving interfaces in the plane Poiseuille flow. When the contact line moves, its dynamics agrees with the Cox theory. In the test case of oscillating drops, themore »contact line transitions smoothly between pinning, advancing and receding. Our method can be directly applied to three-dimensional problems as demonstrated by the test case of sliding drops on an inclined wall.« less
  5. Proper heart morphogenesis requires a delicate balance between hemodynamic forces, myocardial activity, morphogen gradients, and epigenetic signaling, all of which are coupled with genetic regulatory networks. Recently both in vivo and in silico studies have tried to better understand hemodynamics at varying stages of veretebrate cardiogenesis. In particular, the intracardial hemodynamics during the onset of trabeculation is notably complex—the inertial and viscous fluid forces are approximately equal at this stage and small perturbations in morphology, scale, and steadiness of the flow can lead to significant changes in bulk flow structures, shear stress distributions, and chemical morphogen gradients. The immersed boundary method was used to numerically simulate fluid flow through simplified two-dimensional and stationary trabeculated ventricles of 72, 80, and 120 h post fertilization wild type zebrafish embryos and ErbB2-inhibited embryos at seven days post fertilization. A 2D idealized trabeculated ventricular model was also used to map the bifurcations in flow structure that occur as a result of the unsteadiness of flow, trabeculae height, and fluid scale ( R e ). Vortex formation occurred in intertrabecular regions for biologically relevant parameter spaces, wherein flow velocities increased. This indicates that trabecular morphology may alter intracardial flow patterns and hence ventricular shear stressesmore »and morphogen gradients. A potential implication of this work is that the onset of vortical (disturbed) flows can upregulate Notch1 expression in endothelial cells in vivo and hence impacts chamber morphogenesis, valvulogenesis, and the formation of the trabeculae themselves. Our results also highlight the sensitivity of cardiac flow patterns to changes in morphology and blood rheology, motivating efforts to obtain spatially and temporally resolved chamber geometries and kinematics as well as the careful measurement of the embryonic blood rheology. The results also suggest that there may be significant changes in shear signalling due to morphological and mechanical variation across individuals and species.« less