DOI: https://doi.org/10.1515/cmam-2025-0143
تاريخ النشر: 2026-02-20
المؤلف: M. M. Chernyshov وآخرون
الموضوع الرئيسي: هندسة الحاسوب وتوليد الشبكات
نظرة عامة
تتناول الورقة البحثية الطرق الحسابية لحل مشاكل القيمة الحدية متعددة الأبعاد باستخدام الشبكات غير المنتظمة وتقريبات الحجم المحدود. يركز المؤلفون على دمج مثلثات ديلوني ومخططات فورونوي لإنشاء شبكات مدمجة من فورونوي-ديلوني (MVD)، حيث تتكون العقد من كل من رؤوس ديلوني ورؤوس فورونوي، مما يؤدي إلى خلايا رباعية قائمة الزاوية. يسمح هذا الإطار بالتقريب الفعال للدوال العددية والمتجهة، فضلاً عن مشغلات التدرج والتباعد الثابتة الضرورية في حساب التفاضل والتكامل المتجه.
تشمل النتائج الرئيسية بناء مخططات فورونوي من مثلثات ديلوني حادة الزاوية، وتنفيذ استراتيجيتين لتقسيم المجال: تقريبات مركزية للعقد باستخدام مثلثات ديلوني وتقريبات مركزية للخلايا باستخدام مضلعات فورونوي. يتم استخدام طريقة الحجم المحدود لتقريب كل من الحل والتدفقات، مع تعريف أحجام التحكم بواسطة خلايا فورونوي لعقد ديلوني والعكس بالعكس. كما تتناول الورقة تطوير نظائر منفصلة لمشغلات التدرج والتباعد على شبكات MVD، مع الالتزام بمنهجية التمييز المماثل. يتم التحقق من التكنولوجيا الحسابية المقترحة من خلال تجارب عددية على مشكلة انتشار-تفاعل ثنائية الأبعاد في وسائط غير متجانسة، مما يظهر دقتها وقوتها وكفاءتها.
مقدمة
تناقش مقدمة الورقة تطبيق النماذج الرياضية في تحليل العمليات الفيزيائية والكيميائية والبيولوجية من خلال أنظمة من المعادلات التفاضلية الجزئية المرتبطة (PDEs). يمكن أن تكون هذه المعادلات التفاضلية الجزئية إما بيضاوية أو بارابولية أو فرط بيضاوية، وعادة ما يتم حلها باستخدام طرق عددية في مجالات متعددة الأبعاد تتميز بأشكال هندسية معقدة. تسلط الورقة الضوء على الاستخدام الشائع للشبكات غير المنتظمة، وخاصة الشبكات مثلثية الشكل وشبكات رباعية الأبعاد، في طرق العناصر المحدودة (FEM) وطرق الحجم المحدود (FVM) لمشاكل القيمة الحدية. كما تؤكد على أهمية بناء أحجام التحكم وتأسيس قوانين الحفظ على المستوى المنفصل، خاصة من خلال استخدام الشبكات غير المنتظمة وطريقة التوازن.
يقدم المؤلفون مفهوم تكنولوجيا الشبكة المزدوجة، باستخدام كل من شبكات ديلوني وفورونوي لتعزيز قابلية حل المشاكل، خاصة في سياق تقريبات الحجم المحدود لمشكلة انتشار-تفاعل ثنائية الأبعاد في وسط غير متجانس. توضح المقدمة هيكل الورقة، التي تشمل أقسامًا حول صياغة مشكلة القيمة الحدية، ووصف أنواع الشبكات، وتطوير تقريبات المشغلين، وعرض النتائج العددية، مما يظهر في النهاية كفاءة نهج FVM المقترح على الشبكات المدمجة من فورونوي-ديلوني.
نقاش
في هذا القسم، يتناول المؤلفون مشكلة الانتشار-التفاعل الثابتة الخطية في مجال محدود \(\Omega \subset \mathbb{R}^2\)، المميزة بالمعادلة
\[
-\text{div}(K(x) \nabla u(x)) + r(x) u(x) = f(x), \quad x \in \Omega,
\]
مع شروط حدودية متجانسة من نوع ديريشلي \(u(x) = 0\) على \(\partial \Omega\). معامل الانتشار \(K(x)\) هو موتر متماثل وموجب محدد، بينما معامل التفاعل \(r(x)\) غير سالب. يستنتج المؤلفون تقديرات مسبقة، لا سيما
\[
\|\nabla u\| \leq c \kappa \|f\|,
\]
والتي ستكون حاسمة في تحليل التقريبات المنفصلة. يتم استخدام طريقة الحجم المحدود لبناء تقريبات للمشكلة، مع تقديم متغيرات متجهة مساعدة وتعريف متجه التدفق كـ \(q(x) = -K(x) \nabla u(x)\). تتضمن الطريقة إنشاء شبكة أساسية للمتغير \(u(x)\) وشبكة منفصلة للتدفق \(q(x)\)، مما يؤدي إلى نهج منهجي لتفكيك المعادلات الحاكمة.
يتناول القسم أيضًا بناء شبكات ديلوني وفورونوي، مما يؤدي إلى شبكة مدمجة من فورونوي-ديلوني (MVD). تسهل هذه الشبكة MVD تعريف دوال الشبكة العددية والمتجهة، حيث يقوم المؤلفون بإنشاء فضاءات المنتج الداخلي لهذه الدوال. يتم تقريب مشغلات التدرج والتباعد ضمن هذا الإطار، مما يضمن التناسق مع المبادئ الفيزيائية الأساسية. يختتم المؤلفون بعرض تجارب عددية تتحقق من دقة وقوة طريقة الحجم المحدود المقترحة، خاصة في سياق الوسائط غير المتجانسة، مما يظهر أن الخطأ في الحل التقريبي ينخفض مع تحسين الشبكة ويبرز تأثير عدم تجانس الوسط على دقة الحل.
DOI: https://doi.org/10.1515/cmam-2025-0143
Publication Date: 2026-02-20
Author(s): M. M. Chernyshov et al.
Primary Topic: Computational Geometry and Mesh Generation
Overview
The research paper discusses computational methods for solving multidimensional boundary value problems using irregular meshes and finite-volume approximations. The authors focus on the integration of Delaunay triangulations and Voronoi diagrams to create merged Voronoi-Delaunay (MVD) meshes, where the nodes consist of both Delaunay vertices and Voronoi vertices, resulting in orthodiagonal quadrilateral cells. This framework allows for effective approximation of scalar and vector functions, as well as the invariant gradient and divergence operators essential in vector calculus.
Key findings include the construction of Voronoi diagrams from acute-angled Delaunay triangulations, and the implementation of two domain partitioning strategies: node-centered approximations using Delaunay triangles and cell-centered approximations using Voronoi polygons. The finite volume method is employed to approximate both the solution and fluxes, with control volumes defined by Voronoi cells for Delaunay nodes and vice versa. The paper also details the development of discrete analogues of gradient and divergence operators on MVD meshes, adhering to the Mimetic Discretization methodology. The proposed computational technology is validated through numerical experiments on a two-dimensional diffusion-reaction problem in anisotropic media, demonstrating its accuracy, robustness, and efficiency.
Introduction
The introduction of the paper discusses the application of mathematical models in analyzing physical, chemical, and biological processes through systems of coupled partial differential equations (PDEs). These PDEs, which can be elliptic, parabolic, or hyperbolic, are typically solved using numerical methods in multidimensional domains characterized by complex geometries. The paper highlights the prevalent use of unstructured meshes, particularly triangular and tetrahedral meshes, in finite element methods (FEM) and finite volume methods (FVM) for boundary value problems. It emphasizes the significance of constructing control volumes and establishing conservation laws at the discrete level, particularly through the use of irregular meshes and the balance method.
The authors introduce the concept of dual mesh technology, utilizing both Delaunay and Voronoi meshes to enhance the solvability of problems, particularly in the context of finite-volume approximations for a two-dimensional diffusion-reaction problem in an anisotropic medium. The introduction outlines the structure of the paper, which includes sections on the formulation of the boundary value problem, descriptions of the mesh types, development of operator approximations, and presentation of numerical results, ultimately demonstrating the efficiency of the proposed FVM approach on merged Voronoi-Delaunay meshes.
Discussion
In this section, the authors address the linear steady-state diffusion-reaction problem in a bounded domain \(\Omega \subset \mathbb{R}^2\), characterized by the equation
\[
-\text{div}(K(x) \nabla u(x)) + r(x) u(x) = f(x), \quad x \in \Omega,
\]
with homogeneous Dirichlet boundary conditions \(u(x) = 0\) on \(\partial \Omega\). The diffusion coefficient \(K(x)\) is a symmetric, positive-definite tensor, while the reaction coefficient \(r(x)\) is nonnegative. The authors derive a priori estimates, notably
\[
\|\nabla u\| \leq c \kappa \|f\|,
\]
which will be pivotal in analyzing discrete approximations. The finite-volume method is employed to construct approximations of the problem, introducing auxiliary vector variables and defining the flux vector as \(q(x) = -K(x) \nabla u(x)\). The method involves creating a primary mesh for the variable \(u(x)\) and a separate mesh for the flux \(q(x)\), leading to a systematic approach for discretizing the governing equations.
The section further details the construction of Delaunay and Voronoi meshes, culminating in a merged Voronoi-Delaunay (MVD) mesh. This MVD mesh facilitates the definition of scalar and vector mesh functions, with the authors establishing inner product spaces for these functions. The gradient and divergence operators are approximated within this framework, ensuring consistency with the underlying physical principles. The authors conclude by presenting numerical experiments that validate the accuracy and robustness of the proposed finite-volume method, particularly in the context of anisotropic media, demonstrating that the error in the approximate solution decreases with mesh refinement and highlighting the impact of medium anisotropy on solution accuracy.