skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • Passive concentration dynamics incorporated into the library IB2d, a two-dimensional implementation of the immersed boundary method
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 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.  more » « less
Award ID(s):
2111765 1853608
PAR ID:
10321896
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Bioinspiration & Biomimetics
Volume:
17
Issue:
3
ISSN:
1748-3182
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; (, The Journal of the Acoustical Society of America)
    Propagation and attenuation of sound through a layered phononic crystal with viscous constituents is theoretically studied. The Navier–Stokes equation with appropriate boundary conditions is solved and the dispersion relation for sound is obtained for a periodic layered heterogeneous structure where at least one of the constituents is a viscous fluid. Simplified dispersion equations are obtained when the other component of the unit is either elastic solid, viscous fluid, or ideal fluid. The limit of low frequencies when periodic structure homogenizes and the frequencies close to the band edge when propagating Bloch wave becomes a standing wave are considered and enhanced viscous dissipation is calculated. Angular dependence of the attenuation coefficient is analyzed. It is shown that transition from dissipation in the bulk to dissipation in a narrow boundary layer occurs in the region of angles close to normal incidence. Enormously high dissipation is predicted for solid–fluid structure in the region of angles where transmission practically vanishes due to appearance of so-called “transmission zeros,” according to El Hassouani, El Boudouti, Djafari-Rouhani, and Aynaou [Phys. Rev. B 78, 174306 (2008)]. For the case when the unit cell contains a narrow layer of high viscosity fluid, the anomaly related to acoustic manifestation of Borrmann effect is explained. 
    more » « less
  2. ; (, Nonlinearity)
    Abstract We consider absorbing chemical reactions in a fluid flow modelled by the coupled advection–reaction–diffusion equations. In these systems, the interplay between chemical diffusion and fluid transportation causes the enhanced dissipation phenomenon. We show that the enhanced dissipation time scale, together with the reaction coupling strength, determines the characteristic time scale of the reaction. 
    more » « less
  3. ; (, Journal of Fluid Mechanics)
    The scattering of surface waves by structures intersecting liquid surfaces is fundamental in fluid mechanics, with prior studies exploring gravity, capillary and capillary–gravity wave interactions. This paper develops a semi-analytical framework for capillary–gravity wave scattering by a fixed, horizontally placed, semi-immersed cylindrical barrier. Assuming linearised potential flow, the problem is formulated with differential equations, conformal mapping and Fourier transforms, resulting in a compound integral equation framework solved numerically via the Nyström method. An effective-slip dynamic contact line model accounting for viscous dissipation links contact line velocity to deviations from equilibrium contact angles, with fixed and free contact lines of no dissipation as limiting cases. The framework computes transmission and reflection coefficients as functions of the Bond number, slip coefficient and barrier radius, validating energy conservation and confirming a$$90^\circ$$phase difference between transmission and reflection in specific limits. A closed-form solution for scattering by an infinitesimal barrier, derived using Fourier transforms, reveals spatial symmetry in the diffracted field, reduced transmission transitioning from gravity to capillary waves and peak contact line dissipation when the slip coefficient matches the capillary wave phase speed. This dissipation, linked to impedance matching at the contact lines, persists across a range of barrier sizes. These results advance theoretical insights into surface-tension-dominated fluid mechanics, offering a robust theoretical framework for analysing wave scattering and comparison with future experimental and numerical studies. 
    more » « less
  4. ; ; ; (, Journal of Fluid Mechanics)
    We present a systematic simulation campaign to investigate the pairwise interaction of two mobile, monodisperse particles submerged in a viscous fluid and subjected to monochromatic oscillating flows. To this end, we employ the immersed boundary method to geometrically resolve the flow around the two particles in a non-inertial reference frame. We neglect gravity to focus on fluid–particle interactions associated with particle inertia and consider particles of three different density ratios aligned along the axis of oscillation. We systematically vary the initial particle distance and the frequency based on which the particles show either attractive or repulsive behaviour by approaching or moving away from each other, respectively. This behaviour is consistently confirmed for the three density ratios investigated, although particle inertia dictates the overall magnitude of the particle dynamics. Based on this, threshold conditions for the transition from attraction to repulsion are introduced that obey the same power law for all density ratios investigated. We furthermore analyse the flow patterns by suitable averaging and decomposition of the flow fields and find competing effects of the vorticity induced by the fluid–particle interactions. Based on these flow patterns, we derive a circulation-based criterion that provides a quantitative measure to categorize the different cases. It is shown that such a criterion provides a consistent measure to distinguish the attractive and repulsive arrangements. 
    more » « less
  5. ; (, Numerical Methods for Partial Differential Equations)
    Abstract In this article we study the stability of explicit finite difference discretization of advection–diffusion equations (ADE) with arbitrary order of accuracy in the context of method of lines. The analysis first focuses on the stability of the system of ordinary differential equations that is obtained by discretizing the ADE in space and then extends to fully discretized methods in combination with explicit Runge–Kutta methods. In particular, we prove that all stable semi‐discretization of the ADE leads to a conditionally stable fully discretized method as long as the time‐integrator is at least first‐order accurate, whereas high‐order spatial discretization of the advection equation cannot yield a stable method if the temporal order is too low. In the second half of the article, the analysis and the stability results are extended to a partially dissipative wave system, which serves as a model for common practice in many fluid mechanics applications that incorporate a viscous stress in the momentum equation but no heat dissipation in the energy equation. Finally, the major theoretical predictions are verified by numerical examples. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email