An official website of the United States government
We propose a novel solid-fluid coupling method to capture the subtle hydrophobic and hydrophilic interactions between liquid, solid, and air at their multi-phase junctions. The key component of our approach is a Lagrangian model that tackles the coupling, evolution, and equilibrium of dynamic contact lines evolving on the interface between surface-tension fluid and deformable objects. This contact-line model captures an ensemble of small-scale geometric and physical processes, including dynamic waterfront tracking, local momentum transfer and force balance, and interfacial tension calculation. On top of this contact-line model, we further developed a mesh-based level set method to evolve the three-phase T-junction on a deformable solid surface. Our dynamic contact-line model, in conjunction with its monolithic coupling system, unifies the simulation of various hydrophobic and hydrophilic solid-fluid-interaction phenomena and enables a broad range of challenging small-scale elastocapillary phenomena that were previously difficult or impractical to solve, such as the elastocapillary origami and self-assembly, dynamic contact angles of drops, capillary adhesion, as well as wetting and splashing on vibrating surfaces.
more »
« less
A Contact Proxy Splitting Method for Lagrangian Solid-Fluid Coupling
We present a robust and efficient method for simulating Lagrangian solid-fluid coupling based on a new operator splitting strategy. We use variational formulations to approximate fluid properties and solid-fluid interactions, and introduce a unified two-way coupling formulation for SPH fluids and FEM solids using interior point barrier-based frictional contact. We split the resulting optimization problem into a fluid phase and a solid-coupling phase using a novel time-splitting approach with augmentedcontact proxies, and propose efficient custom linear solvers. Our technique accounts for fluids interaction with nonlinear hyperelastic objects of different geometries and codimensions, while maintaining an algorithmically guaranteed non-penetrating criterion. Comprehensive benchmarks and experiments demonstrate the efficacy of our method.
more »
« less
- PAR ID:
- 10471686
- Publisher / Repository:
- ACM TOG
- Date Published:
- Journal Name:
- ACM Transactions on Graphics
- Volume:
- 42
- Issue:
- 4
- ISSN:
- 0730-0301
- Page Range / eLocation ID:
- 1 to 14
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We present an Eulerian vortex method based on the theory of flow maps to simulate the complex vortical motions of incompressible fluids. Central to our method is the novel incorporation of the flow-map transport equations forline elements, which, in combination with a bi-directional marching scheme for flow maps, enables the high-fidelity Eulerian advection of vorticity variables. The fundamental motivation is that, compared to impulsem, which has been recently bridged with flow maps to encouraging results, vorticityωpromises to be preferable for its numerical stability and physical interpretability. To realize the full potential of this novel formulation, we develop a new Poisson solving scheme for vorticity-to-velocity reconstruction that is both efficient and able to accurately handle the coupling near solid boundaries. We demonstrate the efficacy of our approach with a range of vortex simulation examples, including leapfrog vortices, vortex collisions, cavity flow, and the formation of complex vortical structures due to solid-fluid interactions.more » « less
-
Abstract We present a new computational fluid dynamics approach for simulating two‐phase flow in hybrid systems containing solid‐free regions and deformable porous matrices. Our approach is based on the derivation of a unique set of volume‐averaged partial differential equations that asymptotically approach the Navier‐Stokes Volume‐of‐Fluid equations in solid‐free regions and multiphase Biot Theory in porous regions. The resulting equations extend our recently developed Darcy‐Brinkman‐Biot framework to multiphase flow. Through careful consideration of interfacial dynamics (relative permeability and capillary effects) and extensive benchmarking, we show that the resulting model accurately captures the strong two‐way coupling that is often exhibited between multiple fluids and deformable porous media. Thus, it can be used to represent flow‐induced material deformation (swelling, compression) and failure (cracking, fracturing). The model's open‐source numerical implementation,hybridBiotInterFoam, effectively marks the extension of computational fluid mechanics into modeling multiscale multiphase flow in deformable porous systems. The versatility of the solver is illustrated through applications related to material failure in poroelastic coastal barriers and surface deformation due to fluid injection in poro‐visco‐plastic systems.more » « less
-
We introduce an accurate and efficient method for characterizing surface wetting and interfacial properties, such as the contact angle made by a liquid droplet on a solid surface, and the vapor–liquid surface tension of a fluid. The method makes use of molecular simulations in conjunction with the indirect umbrella sampling technique to systematically wet the surface and estimate the corresponding free energy. To illustrate the method, we study the wetting of a family of Lennard-Jones surfaces by water. For surfaces with a wide range of attractions for water, we estimate contact angles using our method, and compare them with contact angles obtained using droplet shapes. Notably, our method is able to capture the transition from partial to complete wetting as surface–water attractions are increased. Moreover, the method is straightforward to implement and is computationally efficient, providing accurate contact angle estimates in roughly 5 nanoseconds of simulation time.more » « less
-
ABSTRACT This work proposes a novel approach for coupling non‐isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid‐filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface‐capturing phase‐field fracture method and interface‐tracking thermo‐fluid‐structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase‐field approach to represent fractures and reconstruct the geometry to frame a thermo‐fluid‐structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase‐field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid‐filled fracture with the phase‐field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase‐field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1145/3592115
