Analytical Expression for Spherically Symmetric Photoacoustic Sources: A Unified General Solution (Theoretical Analysis and Derivation)

إليك ملخص تقني مفصل للورقة البحثية بعنوان "تعبير تحليلي للمصادر الصوتية الضوئية ذات التماثل الكروي: حل عام موحد" من قبل لي وآخرون.

1. بيان المشكلة

يعتمد التصوير الصوتي الضوئي (PAI) على حل معادلة الموجة الصوتية الضوئية لإعادة بناء الصور أو تحليل الإشارات. وبينما تُعد الطرق العددية (مثل طريقة الفروق المحدودة في المجال الزمني) شائعة، إلا أنها قد تكون مكلفة حاسوبياً لعمليات المحاكاة الأمامية، خاصة عند تحليل التماثل الكروي الذي يشيع في الأنسجة البيولوجية والهياكل الدقيقة.
تعالج هذه الورقة الحاجة إلى حل تحليلي شامل وموحد للضغط الصوتي المكاني والزماني الناتج عن مصادر صوتية ضوئية ذات توزيعات ضغط أولية ذات تماثل كروي. الهدف هو اشتقاق صيغة عامة قابلة للتطبيق على أي توزيع أولي متماثل كروياً $p_0(r)$ وتوفير حلول مغلقة محددة لأنواع التوزيع الشائعة (المنتظم، الغاوسي، الأسي، والقانون القوي) لتسهيل تصميم الأنظمة السريع وتحليل الإشارات.

2. المنهجية

يشتق المؤلفون الحل بدءاً من معادلة الموجة الصوتية الضوئية الأساسية تحت فرضية التسخين اللحظي. تسير المنهجية عبر الخطوات المنطقية التالية:

• الصياغة الأساسية: تبدأ المشكلة بمعادلة الموجة غير المتجانسة:
  $$\left(\nabla^2 - \frac{1}{v_s^2}\frac{\partial^2}{\partial t^2}\right)p(\mathbf{r}, t) = -\frac{\beta}{C_p}\frac{\partial H(\mathbf{r}, t)}{\partial t}$$
  باستخدام دالة غرين (Green's function)، يتم التعبير عن الحل كتكامل يتضمن توزيع الضغط الأولي $p_0(\mathbf{r}')$ ودالة دلتا ديراك التي تمثل تأخير انتشار الموجة.

• اختزال التماثل الكروي: بافتراض أن $p_0(\mathbf{r}') = p_0(r')$، يقوم المؤلفون بتحويل التكامل المكاني ثلاثي الأبعاد إلى تكامل شعاعي أحادي البعد.

  • استخدموا خاصية القياس لدالة دلتا: $\delta(t - |\mathbf{r}-\mathbf{r}'|/v_s) = v_s \delta(|\mathbf{r}-\mathbf{r}'| - v_s t)$.
  • من خلال توظيف الإحداثيات الكروية وقانون جيب التمام، تم إجراء التكامل الزاوي تحليلياً. أدى ذلك إلى اختزال تكامل الحجم إلى تكامل خطي على الإحداثي الشعاعي $r'$.

• اشتقاق التعبير الموحد:

  • تم تحديد حدود التكامل بناءً على الشرط $|r - v_s t| \le r' \le r + v_s t$.
  • بتطبيق قاعدة ليبنيز (Leibniz's rule) لتفاضل تكامل بحدود متغيرة، اشتق المؤلفون تعبيراً تحليلياً موحداً للضغط $p(r, t)$ لأي $p_0(r)$ تعسفي متماثل كروياً:
    $$p(r, t) = \frac{1}{2r} \left[ (r + v_s t)p_0(r + v_s t) + (r - v_s t)p_0(|r - v_s t|) \right]$$

• اشتقاق حالات محددة: تم تطبيق الصيغة الموحدة على أربعة أنواع من توزيعات الضغط الأولي:

  1. الكرة المنتظمة (Uniform Sphere): توزيع دالة خطوة ضمن نصف قطر $a_0$.
  2. التوزيع الغاوسي (Gaussian Distribution): $p_0(r) \propto e^{-r^2/2\sigma^2}$.
  3. التوزيع الأسي (Exponential Distribution): $p_0(r) \propto e^{-r/a}$.
  4. توزيع القانون القوي (Power-Law Distribution): $p_0(r) \propto (r^2 + a^2)^{-\nu}$.

• تقريب المجال البعيد: تحلل الورقة سلوك هذه الحلول في نظام المجال البعيد ($r \gg$ أبعاد المصدر)، حيث يظهر أن الضغط يهيمن عليه الحد الذي يتضمن $(r - v_s t)$، والذي يمثل الموجة المنتشرة نحو الخارج.

3. المساهمات الرئيسية

• الحل العام الموحد: يوفر اشتقاق المعادلة (15) إطاراً تحليلياً واحداً يعمل مع أي ملف تعريف للضغط الأولي متماثل كروياً، مما يلغي الحاجة لإعادة اشتقاق التكاملات لكل شكل توزيع جديد.
• حلول مغلقة للنماذج الشائعة: توفر الورقة تعبيرات تحليلية صريحة وجاهزة للاستخدام لأكثر أربعة نماذج مصادر شيوعاً في التصوير الصوتي الضوئي (المنتظم، الغاوسي، الأسي، والقانون القوي)، بما في ذلك التعريفات الجزئية للكرة المنتظمة (داخل مقابل خارج نقاط المراقبة).
• تبسيط المجال البعيد: يقدم المؤلفون تقريبات مبسطة للمجال البعيد لهذه التوزيعات، وهي ضرورية لفهم خصائص الإشارة عند مسافات كبيرة (مثل سيناريوهات التصوير السريري).
• التنفيذ مفتوح المصدر: يوفر المؤلفون مستودع GitHub (SlingBAG Ultra) يحتوي على كود لمحاكاة أمامية فائقة السرعة بناءً على هذه النماذج التحليلية، مما يتيح النمذجة السريعة والتحقق من الصحة.

4. النتالئج

• الكرة المنتظمة: يكشف الحل عن سلوكيات زمنية متميزة اعتماداً على ما إذا كان المراقب داخل أو خارج الكرة. بالنسبة لمراقب خارج الكرة ($r > a_0$)، تكون الإشارة غير صفرية فقط خلال النافذة الزمنية $r-a_0 \le v_s t \le r+a_0$، مع اعتماد خطي على $(r - v_s t)$.
• المصدر الغاوسي: تحافظ موجة الضغط الناتجة على ملف تعريف غاوسي ولكنها تُعدل بواسطة التوهين الهندسي $1/r$ وحد خطي $(r-v_s t)$.
• التوزيع الأسي والقانون القوي: تظهر التعبيرات المشتقة كيف تترجم خصائص الاضمحلال المحددة للضغط الأولي (الأسي مقابل القانون القوي) مباشرة إلى الشكل الزمني للإشارة الصوتية المستلمة.
• سلوك المجال البعيد: في المجال البعيد، يصبح مكون الموجة "المتقاربة نحو الداخل" (الحد الأول في المعادلة الموحدة) ضئيلاً بسبب الانتشار الهندسي والمسافة، مما يترك مكون "الانتشار نحو الخارج" كإشارة مهيمنة.

5. الأهمية

• تصميم وتحسين الأنظمة: تسمح التعبيرات التحليلية للباحثين بمحاكاة الإشارات الصوتية الضوئية بسرعة دون العبء الحسابي للمحركات العددية. وهذا أمر حيوي لتحسين وضع الكواشف، ومتطلبات عرض النطاق الترددي، وخوارزميات إعادة البناء.
• تحليل الإشارات: من خلال توفير حلول دقيقة، تساعد الورقة في فهم الشكل الجوهري للإشارات الصوتية الضوئية الناتجة عن أنواع مختلفة من الأنسجة (التي يتم نمذجتها بتوزيعات $p_0$ مختلفة)، مما يحسن دقة التصوير الصوتي الضوئي الكمي.
• القيمة التعليمية والنظرية: يعمل الاشتقاق التفصيلي كمرجع صارم للمعالجة الرياضية لانتشار الموجات متماثلة التناظر الكروي في التصوير الصوتي الضوئي.
• قابلية التكرار: يضمن توفر كود "SlingBAG Ultra" إمكانية تطبيق هذه النتائج النظرية فوراً من قبل المجتمع العلمي لمهام النمذجة الأمامية.

باختصار، تجسر هذه الورقة الفجوة بين فيزياء الموجات الأساسية والتطبيقات الهندسية العملية في التصوير الصوتي الضوئي من خلال توفير مجموعة أدوات تحليلية قوية، موحدة، وفعالة حاسوبياً للمصادر ذات التماثل الكروي.