Ivan C. Christov

Ziyu Wang, Ivan C. Christov

Soft-walled microchannels arise in many applications, ranging from organ-on-a-chip platforms to soft-robotic actuators. However, despite extensive research on their static and dynamic response, the potential failure of these devices has not been addressed. To this end, we explore fluid--structure interaction in microchannels whose compliant top wall is governed by a nonlocal mechanical theory capable of simulating both deformation and material failure. We develop a one-dimensional model by coupling viscous flow under the lubrication approximation to a state-based peridynamic formulation of an Euler--Bernoulli beam. The peridynamic formulation enables the wall to be modeled as a genuinely nonlocal beam, and the integral form of its equation of motion remains valid whether the deformation field is smooth or contains discontinuities. Through the proposed computational model, we explore the steady and time-dependent behaviors of this fluid--peridynamic structure interaction. We rationalize the wave and damping dynamics observed in the simulations through a dispersion (linearized) analysis of the coupled system, finding that, with increasing nonlocal influence, wave propagation exhibits a clear departure from classical behavior, characterized by a gradual suppression of the phase velocity. The main contribution of our study is to outline the potential failure scenarios of the microchannel's soft wall under the hydrodynamic load of the flow. Specifically, we find a dividing curve in the space spanned by the dimensionless Strouhal number (quantifying unsteady inertia of the beam) and the compliance number (quantifying the strength of the fluid--structure coupling) separating scenarios of potential failure during transient conditions from potential failure at the steady load.

Aditya A. Ghodgaonkar, Ivan C. Christov

We discuss the numerical solution of nonlinear parabolic partial differential equations, exhibiting finite speed of propagation, via a strongly implicit finite-difference scheme with formal truncation error $\mathcal{O}\left[(Δx)^2 + (Δt)^2 \right]$. Our application of interest is the spreading of viscous gravity currents in the study of which these type of differential equations arise. Viscous gravity currents are low Reynolds number (viscous forces dominate inertial forces) flow phenomena in which a dense, viscous fluid displaces a lighter (usually immiscible) fluid. The fluids may be confined by the sidewalls of a channel or propagate in an unconfined two-dimensional (or axisymmetric three-dimensional) geometry. Under the lubrication approximation, the mathematical description of the spreading of these fluids reduces to solving the so-called thin-film equation for the current's shape $h(x,t)$. To solve such nonlinear parabolic equations we propose a finite-difference scheme based on the Crank--Nicolson idea. We implement the scheme for problems involving a single spatial coordinate (i.e., two-dimensional, axisymmetric or spherically-symmetric three-dimensional currents) on an equispaced but staggered grid. We benchmark the scheme against analytical solutions and highlight its strong numerical stability by specifically considering the spreading of non-Newtonian power-law fluids in a variable-width confined channel-like geometry (a "Hele-Shaw cell") subject to a given mass conservation/balance constraint. We show that this constraint can be implemented by re-expressing it as nonlinear flux boundary conditions on the domain's endpoints. Then, we show numerically that the scheme achieves its full second-order accuracy in space and time. We also highlight through numerical simulations how the proposed scheme accurately respects the mass conservation/balance constraint.