[论文解读] Stochastic Eulerian-Lagrangian Methods for Fluid-Structure Interactions with Thermal Fluctuations and Shear Boundary Conditions

Paul J. Atzberger|arXiv (Cornell University)|Oct 29, 2009

Rheology and Fluid Dynamics Studies参考文献 32被引用 5

一句话总结

本文提出了一种随机欧拉-拉格朗日方法（SELM），用于模拟软材料在热涨落和剪切流作用下的流固耦合行为。通过结合拉格朗日微结构动力学与欧拉流体动量方程及随机边界条件，该方法实现了对聚合物流体、囊泡悬浮液和凝胶中流变响应的精确、统计一致的模拟，揭示了由微结构重排驱动的剪切变稀与结构稀释行为。

ABSTRACT

A computational approach is introduced for the study of the rheological properties of complex fluids and soft materials. The approach allows for a consistent treatment of microstructure elastic mechanics, hydrodynamic coupling, thermal fluctuations, and externally driven shear flows. A mixed description in terms of Eulerian and Lagrangian reference frames is used for the physical system. Microstructure configurations are represented in a Lagrangian reference frame. Conserved quantities, such as momentum of the fluid and microstructures, are represented in an Eulerian reference frame. The mathematical formalism couples these different descriptions using general operators subject to consistency conditions. Thermal fluctuations are taken into account in the formalism by stochastic driving fields introduced in accordance with the principles of statistical mechanics. To study the rheological responses of materials subject to shear, generalized periodic boundary conditions are developed where periodic images are shifted relative to the unit cell to induce shear. Stochastic numerical methods are developed for the formalism. As a demonstration of the methods, results are presented for the shear responses of a polymeric fluid, lipid vesicle fluid, and a gel-like material.

研究动机与目标

• 开发一个一致的计算框架，用于建模具有微结构弹性的软材料，同时考虑流体动力耦合与热涨落。
• 解决在随机驱动下，欧拉描述与拉格朗日描述耦合时违反守恒定律的挑战，尤其是动量守恒。
• 将周期性边界条件推广至剪切流诱导，采用类似Lees-Edwards的位移图像方法。
• 发展保持统计力学原理的随机数值方法，适用于惯性与过阻尼两种情形。
• 实现对剪切粘度、法向应力及粘弹性模量等流变性质的定量研究，基于微尺度动力学。

提出的方法

• 采用混合欧拉-拉格朗日格式：微结构在拉格朗日坐标系中追踪，而流体动量在欧拉空间中守恒。
• 采用广义周期性边界条件，结合位移周期性图像以诱导宏观剪切流，将Lees-Edwards边界条件扩展至流体动力系统。
• 在流体方程中引入满足涨落-耗散平衡的随机驱动场，确保统计一致性。
• 应用跳跃条件与变量变换，以在数值格式中处理剪切域边界处的速度不连续性。
• 开发高效的随机求解器，用于求解在完全解析与准稳态流体区域中的涨落流场。
• 使用应力估算器（公式63）计算宏观应力张量，该方法同时考虑n体微结构相互作用与广义边界条件。

实验结果

研究问题

• RQ1 在存在热噪声的流固系统中，如何实现欧拉描述与拉格朗日描述的一致耦合？
• RQ2 如何制定广义周期性边界条件，以诱导剪切流，同时保持动量守恒与统计力学原理？
• RQ3 微结构动力学（如链拉伸、取向与键断裂）如何调控软材料的时间依赖剪切粘度？
• RQ4 热涨落如何在剪切凝胶网络中维持非零稳态剪切粘度？
• RQ5 微结构重排与键重组如何导致复杂流体中的结构稀释或剪切变稀行为？

主要发现

• 聚合物流体表现出剪切变稀行为，剪切速率升高时剪切粘度降低，这是由于聚合物链取向与张力减小所致。
• 在囊泡流体中，振荡剪切响应揭示了频率相关的粘弹性响应，其储能模量与损耗模量在宽频率范围内被量化。
• 凝胶状材料表现出三阶段流变响应：初始因链拉伸导致粘度上升，随后因链取向而下降，最终在热涨落维持下保持非零稳态粘度。
• 凝胶的时间依赖剪切粘度初始上升，源于链间键的拉伸，随后因链重取向而下降，并在热噪声引起的瞬时错位作用下稳定于有限值。
• SELM框架通过同时解析微结构力学、流体动力耦合与热涨落，成功捕捉了包括结构稀释与剪切变稀在内的复杂流变现象。
• 应力估算器（公式63）可准确计算宏观应力，涵盖n体相互作用与广义边界条件下的微尺度力。