DOI: https://doi.org/10.1002/fld.70080
تاريخ النشر: 2026-05-15
المؤلف: Dmitry Ammosov وآخرون
الموضوع الرئيسي: النمذجة الرياضية المتقدمة في الهندسة
نظرة عامة
في هذه الورقة، يقوم المؤلفون بتطوير نماذج المرونة المسامية متعددة المستويات باستخدام طريقة تجانس متعددة المستويات لمعالجة التحديات التي تطرحها الخصائص ذات التباين العالي في الوسائط المرنة المسامية. غالبًا ما تفشل تقنيات التجانس التقليدية في تقديم حلول دقيقة بسبب غياب المعلمات الماكروسكوبية المناسبة. يتم استخدام النهج متعدد المستويات، الذي يتضمن تعريف حالات أو مستويين متوسطين متعددة، للتغلب على هذه القيود. يقوم المؤلفون بتوسيع العمل السابق من خلال اشتقاق معادلات متعددة المستويات العامة بدقة من خلال صياغة مشاكل خلايا مقيدة مرتبطة في المناطق المفرطة العينة، مما يسمح بأخذ تأثيرات التجانس المختلفة في الاعتبار.
تقدم الدراسة توسعًا شاملاً متعدد المستويات للحقول الدقيقة وتشتق نموذج مرونة مسامية متعدد المستويات عام قابل للتطبيق على عدد غير محدد من المستويات. بالإضافة إلى ذلك، يتم تقديم نموذجين مبسطين، حيث يكون النموذج المبسط الثاني حالة محددة من نموذج المرونة المسامية متعدد الشبكات القياسي. تظهر التجارب العددية التي أجريت على وسائط غير متجانسة مختلفة أن النماذج متعددة المستويات المقترحة تحقق دقة عالية في تقريب الحلول المتوسطة المرجعية. ومع ذلك، يظهر النموذج المبسط الثاني أخطاء أكبر، خاصة في التقاط خصائص الحل المرجعي في السيناريوهات غير المتجانسة المعقدة.
مقدمة
تناقش مقدمة ورقة البحث أهمية نظرية المرونة المسامية في فهم تدفقات العمليات الميكانيكية المترابطة في الوسائط المسامية، وهو أمر حاسم للتطبيقات مثل التنبؤ بانخفاض الأرض، وتحليل استقرار آبار الحفر، والزلازل الناتجة عن السوائل. اقترحها في الأصل بيوت، تتضمن النظرية نظامًا من المعادلات التفاضلية الجزئية المترابطة التي تصف الإزاحات الميكانيكية وضغط السوائل، مع تأثر الربط بقوى الجسم الميكانيكية وتأثيرات الانفعال الحجمي. ومع ذلك، تواجه المحاكاة العددية لنماذج المرونة المسامية تحديات بسبب الطبيعة متعددة المقاييس للوسائط المسامية، مما يتطلب شبكات حسابية دقيقة تزيد من التكاليف الحاسوبية.
لمعالجة هذه التحديات، تسلط الورقة الضوء على استخدام طرق التجانس لحساب الخصائص الفعالة للوسائط غير المتجانسة، مما يسمح بتطوير نماذج ماكروسكوبية. بينما تكون هذه الطرق فعالة، قد لا تحقق نتائج دقيقة في الحالات ذات الخصائص ذات التباين العالي، مما يستدعي الحاجة إلى نهج متعدد المستويات يحدد حالات متوسطة متعددة. يتم تقديم نظرية المرونة المسامية متعددة الشبكات (MPET) كامتداد لنظرية الوسائط المسامية المتعددة، وهي ذات صلة خاصة في الأنسجة البيولوجية ذات خصائص قنوات السوائل المتغيرة. تقترح الورقة طريقة تجانس متعددة المستويات تصيغ مشاكل خلايا مقيدة في المناطق المفرطة العينة لاشتقاق نماذج متعددة المستويات، والتي يتم التحقق منها من خلال التجارب العددية التي تظهر دقة عالية. يتم توضيح هيكل الورقة، مع تخصيص الأقسام اللاحقة للمقدمات، وتطبيق الطريقة، والأمثلة العددية، وملخص النتائج.
مناقشة
في هذا القسم، يناقش المؤلفون الجوانب الأساسية لنماذج المرونة المسامية، مع التركيز على كل من النموذج القياسي ونموذج المرونة المسامية متعدد الشبكات. يتميز النموذج القياسي، الذي اقترحه بيوت في الأصل، بنظام من المعادلات التفاضلية الجزئية المترابطة التي تحكم مجال متجه الإزاحة \( \mathbf{u} \) ومجال الضغط \( p \). يتم تمثيل المعادلات على النحو التالي:
\[
-\nabla \cdot \sigma(\mathbf{u}) + \alpha \nabla p = f, \quad x \in \Omega,
\]
\[
\alpha \frac{\partial \nabla \cdot \mathbf{u}}{\partial t} + S \frac{\partial p}{\partial t} – \nabla \cdot (\kappa \nabla p) = g, \quad x \in \Omega,
\]
حيث \( \sigma(\mathbf{u}) \) هو موتر الإجهاد، \( \alpha \) هو معامل بيوت، \( S \) هو معامل التخزين، و \( \kappa \) هو معامل النفاذية. يبرز المؤلفون التحديات المتمثلة في تطبيق هذا النموذج على الوسائط المسامية متعددة المقاييس، خاصة في السياقات البيولوجية، حيث تكون الخصائص ذات التباين العالي شائعة. لمعالجة هذه التحديات، يقدمون نموذج المرونة المسامية متعدد الشبكات (MPET)، الذي يتضمن عدة مقصورات سائلة ومعدلات نقل السوائل بين الشبكات، ممثلة على النحو التالي:
\[
-\nabla \cdot \sigma(\mathbf{u}) + \sum_i \alpha_i \nabla p_i = f, \quad x \in \Omega,
\]
\[
\alpha_i \frac{\partial \nabla \cdot \mathbf{u}}{\partial t} + S_i \frac{\partial p_i}{\partial t} + \sum_{j \neq i} S_{ij} \frac{\partial p_j}{\partial t} – \nabla \cdot (\kappa_i \nabla p_i) = g_i, \quad x \in \Omega.
\]
يشير المؤلفون إلى أن معاملات تخزين المسامية المتقاطعة \( S_{ij} \) غالبًا ما يتم تجاهلها بسبب نقص البيانات التجريبية. يقترحون اشتقاق نماذج المرونة المسامية متعددة المستويات باستخدام طريقة تجانس متعددة المستويات، والتي سيتم توضيحها في الأقسام اللاحقة. يهدف هذا النهج إلى تقديم تمثيل أكثر دقة للتفاعلات المعقدة في الوسائط المسامية متعددة المقاييس، خاصة في الأنسجة البيولوجية.
DOI: https://doi.org/10.1002/fld.70080
Publication Date: 2026-05-15
Author(s): Dmitry Ammosov et al.
Primary Topic: Advanced Mathematical Modeling in Engineering
Overview
In this paper, the authors develop multicontinuum poroelasticity models using a multicontinuum homogenization method to address the challenges posed by high-contrast properties in poroelastic media. Traditional homogenization techniques often fail to yield accurate solutions due to the absence of appropriate macroscopic parameters. The multicontinuum approach, which involves defining multiple average states or continua, is employed to overcome these limitations. The authors extend previous work by rigorously deriving generalized multicontinuum equations through the formulation of coupled constraint cell problems in oversampled regions, allowing for the consideration of various homogenized effects.
The study presents a comprehensive multicontinuum expansion of fine-scale fields and derives a general multicontinuum poroelastic model applicable to an arbitrary number of continua. Additionally, two simplified models are introduced, with the second simplified model being a specific case of the standard multiple-network poroelasticity model. Numerical experiments conducted on different heterogeneous media demonstrate that the proposed multicontinuum models achieve high accuracy in approximating reference average solutions. However, the second simplified model exhibits greater errors, particularly in capturing the characteristics of the reference solution in complex heterogeneous scenarios.
Introduction
The introduction of the research paper discusses the significance of poroelasticity theory in understanding coupled flow and mechanical processes in porous media, which is crucial for applications such as land subsidence prediction, borehole stability analysis, and fluid-induced seismicity. Originally proposed by Biot, the theory involves a system of coupled partial differential equations that describe mechanical displacements and fluid pressure, with the coupling influenced by mechanical body forces and volumetric strain effects. However, the numerical simulation of poroelasticity models faces challenges due to the multiscale nature of porous media, necessitating fine computational grids that increase computational costs.
To address these challenges, the paper highlights the use of homogenization methods to compute effective properties for heterogeneous media, allowing for the development of macroscopic models. While these methods are effective, they may not yield accurate results in cases with high-contrast properties, prompting the need for multicontinuum approaches that define multiple average states. The multiple-network poroelastic theory (MPET) is introduced as an extension of multiple porous media theory, particularly relevant in biological tissues with varying fluid channel properties. The paper proposes a multicontinuum homogenization method that formulates constraint cell problems in oversampled regions to derive multicontinuum models, which are validated through numerical experiments demonstrating high accuracy. The structure of the paper is outlined, with subsequent sections dedicated to preliminaries, method application, numerical examples, and a summary of findings.
Discussion
In this section, the authors discuss the foundational aspects of poroelasticity models, focusing on both the standard and multiple-network poroelasticity models. The standard model, originally proposed by Biot, is characterized by a system of coupled partial differential equations governing the displacement vector field \( \mathbf{u} \) and the pressure field \( p \). The equations are represented as:
\[
-\nabla \cdot \sigma(\mathbf{u}) + \alpha \nabla p = f, \quad x \in \Omega,
\]
\[
\alpha \frac{\partial \nabla \cdot \mathbf{u}}{\partial t} + S \frac{\partial p}{\partial t} – \nabla \cdot (\kappa \nabla p) = g, \quad x \in \Omega,
\]
where \( \sigma(\mathbf{u}) \) is the stress tensor, \( \alpha \) is the Biot coefficient, \( S \) is the storage coefficient, and \( \kappa \) is the permeability coefficient. The authors highlight the challenges of applying this model to multiscale porous media, particularly in biological contexts, where high-contrast properties are prevalent. To address these challenges, they introduce the multiple-network poroelasticity model (MPET), which incorporates multiple fluid compartments and fluid transfer rates between networks, represented as:
\[
-\nabla \cdot \sigma(\mathbf{u}) + \sum_i \alpha_i \nabla p_i = f, \quad x \in \Omega,
\]
\[
\alpha_i \frac{\partial \nabla \cdot \mathbf{u}}{\partial t} + S_i \frac{\partial p_i}{\partial t} + \sum_{j \neq i} S_{ij} \frac{\partial p_j}{\partial t} – \nabla \cdot (\kappa_i \nabla p_i) = g_i, \quad x \in \Omega.
\]
The authors note that the cross-porosity storage coefficients \( S_{ij} \) are often neglected due to insufficient experimental data. They propose to derive multicontinuum poroelasticity models using a multicontinuum homogenization method, which will be elaborated in subsequent sections. This approach aims to provide a more accurate representation of the complex interactions in multiscale porous media, particularly in biological tissues.