An official website of the United States government
Boundary conditions between a porous solid and a fluid has been a long-standing problem in modeling porous media. For deformable poroelastic materials such as hydrogels, the question is further complicated by the elastic stress from the solid network. Recently, an interfacial permeability condition has been developed from the principle of positive energy dissipation on the hydrogel–fluid interface. Although this boundary condition has been used in flow computations and yielded reasonable predictions, it contains an interfacial permeability g as a phenomenological parameter. In this work, we use porescale models of flow into a periodic array of solid cylinders or parallel holes to determine g as a function of the pore size and porosity. This provides a means to evaluate the interfacial permeability for a wide range of poroelastic materials, including hydrogels, foams and biological tissues, to enable realistic flow simulations.
more »
« less
Evaluating impedance boundary conditions to model interfacial dynamics in acoustofluidics
We present a numerical study to investigate the efficacy of impedance boundary conditions in capturing the interfacial dynamics of a particle subjected to an acoustic field and study the concomitant time-averaged acoustic streaming and radiation force fields. While impedance boundary conditions have been utilized to represent fluid–solid interface in acoustofluidics, such models assume the solid material to be locally reactive to the acoustic waves. However, there is a limited understanding of when this assumption holds true, raising concerns about the suitability of impedance boundary conditions. Here, we systematically investigate the applicability of impedance boundary conditions by comparing the predictions of an impedance boundary approach against a fully coupled fluid–solid model. We contrast the oscillation profiles of the fluid–solid interface predicted by the two models. We consider different scatterer materials to identify the extent to which the differences in interfacial dynamics impact the time-averaged fields and highlight the divergence within the predictions of the two models. Our findings indicate that, although impedance boundary conditions can yield qualitatively similar results to the full model in certain situations, the predictions from the two models generally differ both qualitatively and quantitatively. These results underscore the importance of exercising caution when applying these boundary conditions to model general acoustofluidic systems.
more »
« less
- Award ID(s):
- 2229636
- PAR ID:
- 10570089
- Publisher / Repository:
- AIP Publishing
- Date Published:
- Journal Name:
- Physics of Fluids
- Volume:
- 36
- Issue:
- 8
- ISSN:
- 1070-6631
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
This study explores the Faraday instability as a mechanism to enhance heat transfer in two-phase systems by exciting interfacial waves through resonance. The approach is particularly applicable to reduced-gravity environments where buoyancy-driven convection is ineffective. A reduced-order model, based on a weighted residual integral boundary layer method, is used to predict interfacial dynamics and heat flux under vertical oscillations with a stabilising thermal gradient. The model employs long-wave and one-way coupling approximations to simplify the governing equations. Linear stability theory informs the oscillation parameters for subsequent nonlinear simulations, which are then qualitatively compared against experiments conducted under Earth’s gravity. Experimental results show up to a 4.5-fold enhancement in heat transfer over pure conduction. Key findings include: (i) reduced gravity lowers interfacial stability, promoting mixing and heat transfer; and (ii) oscillation-induced instability significantly improves heat transport under Earth’s gravity. Theoretical predictions qualitatively validate experimental trends in wavelength-dependent enhancement of heat transfer. Quantitative discrepancies between model and experiment are rationalised by model assumptions, such as neglecting higher-order inertial terms, idealised boundary conditions, and simplified interface dynamics. These limitations lead to underprediction of interface deflection and heat flux. Nevertheless, the study underscores the value of Faraday instability as a means to boost heat transfer in reduced gravity, with implications for thermal management in space applications.more » « less
-
We develop a model for solid objects surrounded by a fluid that accounts for the possibility of acoustic pressures giving rise to damage on the surface of the solid. The propagation of an acoustic pressure in the fluid domain is modeled by the acoustic wave equation. On the other hand, the response of the solid is described by linear elastodynamics coupled with a gradient damage model, one that is based on a cohesive-type phase-field description of fracture. The interaction between the acoustic pressure and the deformation and damage of the solid are represented by transmission conditions at the fluid-solid interface. The resulting governing equations are discretized using a finite-element/finite-difference method that pays particular attention to the spatial and temporal scales that need to be resolved. Results from model-based simulations are provided for a benchmark problem as well as for recent experiments in nanopulse lithotripsy. A parametric study is performed to illustrate how damage develops in response to the driving force (magnitude and location of the acoustic source) as a function of the fracture resistance of the solid. The results are shown to be qualitatively consistent with experimental observations for the location and size of the damage fields on the solid surface. A study of limiting cases also suggests that both the threshold for damage and the critical fracture energy are important to consider in order to capture the transition from damage initiation to complete localization. A low-cycle fatigue model is proposed that degrades the fracture resistance of the solid as a function of accumulated tensile strain energy, and it is shown to be capable of capturing damage localization in simulations of multi-pulse nano-pulse lithotripsy.more » « less
-
null (Ed.)The fluid dynamics of a bubble collapsing near an elastic or viscoelastic material is coupled with the mechanical response of the material. We apply a multiphase fluid–solid coupled computational model to simulate the collapse of an air bubble in water induced by an ultrasound shock wave, near different types of materials including metals (e.g. aluminium), polymers (e.g. polyurea), minerals (e.g. gypsum), glass and foams. We characterize the two-way fluid–material interaction by examining the fluid pressure and velocity fields, the time history of bubble shape and volume and the maximum tensile and shear stresses produced in the material. We show that the ratio of the longitudinal acoustic impedance of the material compared to that of the ambient fluid, $$Z/Z_0$$ , plays a significant role. When $$Z/Z_0<1$$ , the material reflects the compressive front of the incident shock into a tensile wave. The reflected tensile wave impinges on the bubble and decelerates its collapse. As a result, the collapse produces a liquid jet, but not necessarily a shock wave. When $$Z/Z_0>1$$ , the reflected wave is compressive and accelerates the bubble's collapse, leading to the emission of a shock wave whose amplitude increases linearly with $$\log (Z/Z_0)$$ , and can be much higher than the amplitude of the incident shock. The reflection of this emitted shock wave impinges on the bubble during its rebound. It reduces the speed of the bubble's rebound and the velocity of the liquid jet. Furthermore, we show that, for a set of materials with $$Z/Z_0\in [0.04, 10.8]$$ , the effect of acoustic impedance on the bubble's collapse time and minimum volume can be captured using phenomenological models constructed based on the solution of Rayleigh–Plesset equation.more » « less
-
Interactions between an evolving solid and inviscid flow can result in substantial 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. Two-dimensional results are presented as well 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.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0225930
