معادلة إمدن شاندراسيخار

تصغير|200بك|يسار
في الفيزياء الفلكية، تُعد **معادلة إمدن شاندراسيخار** شكلًا لا بعدي لمعادلة بواسون لتوزيع كثافة كرة غازية متساوية في درجة الحرارة ومتماثلة كرويًا تخضع لقوة الجاذبية الخاصة بها، سُميت نسبةً لروبرت إمدن وسابرامانين شاندراسيخار. قُدمت المعادلة لأول مرة من قبل روبرت إمدن في عام 1907. المعادلة هي:

\frac{1}{\xi^2} \frac{d}{d\xi}\left(\xi^2 \frac{d\psi}{d\xi}\right)= e^{-\psi}

حيث \xi هو نصف القطر اللا بعدي و\psi ترتبط بكثافة كرة الغاز التي تُعطى بالعلاقة \rho=\rho_c e^{-\psi}، حيث \rho_c هي كثافة الغاز في المركز. المعادلة ليس لها حل واضح معروف. إذا استُخدم سائل بوليتروبي بدلًا من سائل متساوي الحرارة، يمكن اشتقاق معادلة لين إمدن. عادةً ما يُستخدم افتراض تساوي الحرارة لوصف لبّ النجم. تُحل المعادلة بالشروط الأولية التالية:

\psi =0, \quad \frac{d\psi}{d\xi} =0 \quad \text{at} \quad \xi=0.

تظهر المعادلة في فروع أخرى من الفيزياء أيضًا، على سبيل المثال تظهر نفس المعادلة في نظرية انفجار فرانك كامينيتسكي للوعاء الكروي. دُرست النسخة النسبية لهذا النموذج متساوي الحرارة المتماثل كرويًا من قبل سوبرامنن شاندراسيخار عام 1972.

## الاشتقاق
بالنسبة للنجم الغازي متساوي الحرارة، ينتج الضغط p عن الضغط الحركي والضغط الإشعاعي

p = \rho\frac{k_B}{W H} T + \frac{4\sigma}{3c} T^4

حيث

* \rho هي الكثافة
* k_B هو ثابت بولتزمان
* W هو متوسط الوزن الجزيئي
* H هي كتلة البروتون
* T هي درجة حرارة النجم
* \sigma هو ثابت ستيفان بولتزمان
* c هي سرعة الضوء

تتطلب معادلة توازن النجم توازنًا بين قوة الضغط وقوة الجاذبية

\frac{1}{r^2} \frac{d}{dr} \left(\frac{r^2}{\rho}\frac{dp}{dr}\right)= - 4\pi G \rho

حيث r هو نصف القطر المُقاس من المركز، وG هو ثابت الجاذبية. يُعاد كتابة المعادلة كما يلي

\frac{k_B T}{WH}\frac{1}{r^2} \frac{d}{dr} \left(r^2\frac{d \ln \rho}{dr} \right) = - 4\pi G \rho

باستخدام التحويل التالي

\psi = \ln \frac{\rho_c}{\rho}, \quad \xi = r \left(\frac{4\pi G \rho_c W H}{k_B T}\right)^{1/2}

حيث \rho_c هي الكثافة المركزية للنجم، تُصبح المعادلة

\frac{1}{\xi^2} \frac{d}{d\xi}\left(\xi^2 \frac{d\psi}{d\xi}\right)= e^{-\psi}

الشروط الحدية هي

\psi =0, \quad \frac{d\psi}{d\xi} =0 \quad \text{at} \quad \xi=0

لأن \xi\ll 1، يصبح الحل

\psi = \frac{\xi^2}{6} - \frac{\xi^4}{120} + \frac{\xi^6}{1890} + \cdots 

## محددات النموذج
بافتراض أن الكرة متساوية الحرارة تتمتع ببعض السلبيات. على الرغم من أن الكثافة التي نحصل عليها كحل لهذه الكرة الغازية متساوية الحرارة تنخفض ابتداءً من المركز، فهي تنخفض ببطء شديد بدرجة تحول من تمتع الكرة بسطح مُحدد جيدًا وكتلة محدودة. يمكن اثبات أنه عندما \xi \gg 1،

\frac{\rho}{\rho_c}=e^{-\psi}=\frac{2}{\xi^2} \left[1+\frac{A}{\xi^{1/2}} \cos\left(\frac{\sqrt 7}{2}\ln \xi + \delta\right) + O(\xi^{-1})\right]

حيث A و \delta هي ثوابت يمكن معرفة قيمتها بالحل العددي. يؤدي سلوك الكثافة هذا إلى زيادة الكتلة مع زيادة نصف القطر. بالتالي، يكون النموذج عادةً صالحًا لوصف لبّ النجم، حيث تكون درجة الحرارة ثابتةً تقريبًا.

## الحل المفرد
باستخدام التحويل x=1/\xi، تُحول المعادلة إلى

x^4 \frac{d^2\psi}{dx^2} = e^{-\psi}

تمتلك المعادلة حلًا مفردًا يُعطى كما يلي

e^{-\psi_s} = 2x^2, \quad \text{or} \quad -\psi_s = 2 \ln x+ \ln 2

لذلك، يمكن استخدام متغير جديد -\psi = 2 \ln x + z، حيث يمكن اشتقاق معادلة z،

\frac{d^2z}{dt^2}-\frac{dz}{dt}+ e^z -2 =0, \quad \text{where} \quad t=\ln x

يمكن تخفيض هذه المعادلة إلى معادلة من الدرجة الأولى من خلال استخدام العلاقة التالية

y=\frac{dz}{dt} = \xi \frac{d\psi}{d\xi} - 2

لينتج

y\frac{dy}{dz} - y + e^z- 2 = 0

## التخفيض
هناك تخفيض آخر يُنسب إلى إدوارد آرثر ميلن. دعونا نُعرف العلاقة التالية

u = \frac{\xi e^{-\psi}}{d\psi/d\xi}, \quad v = \xi \frac{d\psi}{d\xi} 
ثم

\frac{u}{v}\frac{dv}{du} = \frac{u-1}{u+v-3}

## الخصائص

* إذا كان \psi(\xi) هو حل لمعادلة إمدن شاندراسيخار، فإن \psi(A\xi)-2\ln A هو أيضًا حل للمعادلة، حيث A هو ثابت اعتباطي.
* حلول معادلة إمدن شاندراسيخار التي تكون محدودة في نقطة الأصل هي بالضرورة d\psi/d\xi=0 عندما \xi=0.

## انظر أيضًا

* معادلة لين-إمدن

## مراجع

تصنيف:الحركيات النجمية
تصنيف:جريان الموائع
تصنيف:معادلات تفاضلية عادية
تصنيف:معادلات فيزيائية