स्मोलुचोव्स्की जमावट समीकरण

सांख्यिकीय भौतिकी में, स्मोलुचोव्स्की जमावट समीकरण एक जनसंख्या संतुलन समीकरण है जो मैरियन स्मोलुचोव्स्की द्वारा 1916 के एक मौलिक प्रकाशन में पेश किया गया था, कणों की संख्या घनत्व के समय के विकास का वर्णन करते हुए वे जमावट करते हैं (इस संदर्भ में "एक साथ टकराते हुए") समय t पर x आकार में ।

बहुलकीकरण, एयरोसौल्ज़ का सहसंयोजन, पायसीकरण, फ्लोकुलेशन से जुड़ी प्रक्रियाओं में एक साथ जमावट (या एकत्रीकरण) का सामना करना पड़ता है

समीकरण
तंत्र के सभी कणों के परस्पर संबंध के अनुसार कण आकार का वितरण समय में बदलता है। इसलिए, स्मोलुचोव्स्की जमावट समीकरण कण-आकार के वितरण का एक पूर्णांक समीकरण है। कारको में जब जमा हुए कणों के आकार निरंतर चर होते हैं, तो समीकरण में एक अभिन्न अंग सम्मलित होता है:

यदि dy की व्याख्या असतत माप के रूप में की जाती है, अर्थात जब कण असतत आकारों में जुड़ते हैं, तो समीकरण का असतत रूप योग होता है:


 * $$\frac{\partial n(x_i,t)}{\partial t}=\frac{1}{2}\sum^{i-1}_{j=1}

K(x_i-x_j,x_j)n(x_i-x_j,t)n(x_j,t) - \sum^\infty_{j=1}K(x_i,x_j)n(x_i,t)n(x_j,t).$$ चुने गए कर्नेल प्रकार्य के लिए एक अद्वितीय उपाय मौजूद है।

जमावट गिरी
रैखिक ऑपरेटर, K, को कोएग्युलेशन अभिन्न कर्नेल के रूप में जाना जाता है और उस दर का वर्णन करता है जिस पर आकार के कण होते हैं $$x_1$$ आकार के कणों के साथ जमना $$x_2$$. समीकरण के लिए बंद रूप की अभिव्यक्ति तब मौजूद होती है जब कर्नेल तीन सरल रूपों में से एक लेता है:


 * $$K = 1,\quad K = x_1 + x_2, \quad K = x_1x_2,$$

स्थिरांक (गणित), योज्य मानचित्र, और गुणात्मक फलन गुठली के रूप में जाना जाता है। कारक के लिए $$K = 1$$ यह गणितीय रूप से सिद्ध किया जा सकता है कि Smoluchowski जमावट समीकरणों के समाधान में विषम रूप से गतिशील स्केलिंग संपत्ति है। यह स्व-समान व्यवहार स्केल इनवेरियन से निकटता से संबंधित है जो एक चरण संक्रमण की एक विशेषता विशेषता हो सकती है।

हालांकि, अधिकांश व्यावहारिक अनुप्रयोगों में कर्नेल काफी अधिक जटिल रूप धारण कर लेता है। उदाहरण के लिए, मुक्त-आण्विक कर्नेल जो तनु गैस-चरण (पदार्थ) प्रणाली में टकराव का वर्णन करता है,


 * $$K = \sqrt{\frac{\pi k_B T}{2}}\left(\frac{1}{m(x_1)}+\frac{1}{m(x_2)}\right)^{1/2}\left(d(x_1)+d(x_2)\right)^2.$$

कुछ जमावट गुठली समूहों के एक विशिष्ट भग्न आयाम के लिए खाते हैं, जैसा कि प्रसार-सीमित एकत्रीकरण में है:
 * $$K = \frac{2}{3} \frac{ k_B T} {\eta} \left(x_1^{1/y_1} +x_2^{1/y_2}\right)\left(x_1^{-1/y_1} +x_2^{-1/y_2}\right),$$

या रिएक्शन-सीमित एकत्रीकरण:
 * $$K = \frac{2}{3} \frac{ k_B T} {\eta} \frac{(x_1x_2)^\gamma}{W}\left(x_1^{1/y_1} +x_2^{1/y_2}\right)\left(x_1^{-1/y_1} +x_2^{-1/y_2}\right),$$

कहाँ $$y_1,y_2$$ समूहों के भग्न आयाम हैं, $$k_B$$ बोल्ट्जमैन स्थिरांक है, $$T$$ तापमान है, $$W$$ फुच्स स्थिरता अनुपात है, $$\eta$$ निरंतर चरण चिपचिपाहट है, और $$\gamma$$ उत्पाद कर्नेल का प्रतिपादक है, जिसे आमतौर पर एक फिटिंग पैरामीटर माना जाता है। बादल के लिए, बादल कणों के जमाव के लिए कर्नेल को आमतौर पर इस प्रकार व्यक्त किया जाता है:
 * $$K = \pi [r(x_1)+r(x_2)]^2 |v(x_1)-v(x_2)| E_{coll}(x_1,x_2), $$ कहाँ $$ r(x)$$ और $$ v(x)$$ बादल के कणों की त्रिज्या और गिरने की गति को आमतौर पर शक्ति कानून का उपयोग करके व्यक्त किया जाता है।

समान्यता इस तरह के भौतिक यथार्थवादी गुठली से उत्पन्न जमावट समीकरण हल करने योग्य नहीं होते हैं, और इस तरह, संख्यात्मक विश्लेषण के लिए अपील करना आवश्यक है। अधिकांश नियतात्मक विधियों का उपयोग तब किया जा सकता है जब ब्याज की केवल एक कण संपत्ति (x) हो, दो प्रमुख हैं क्षणों की विधि (संभाव्यता सिद्धांत)    और अनुभागीय तरीके। बहुभिन्नरूपी कैलकुलस | बहुभिन्नरूपी कारक में, हालांकि, जब दो या दो से अधिक गुण (जैसे आकार, आकार, संरचना, आदि) पेश किए जाते हैं, तो किसी को विशेष सन्निकटन विधियों की तलाश करनी होती है जो आयामीता के अभिशाप से कम पीड़ित होती हैं। गॉसियन रेडियल आधार कार्यों के आधार पर सन्निकटन को एक से अधिक आयामों में जमावट समीकरण पर सफलतापूर्वक लागू किया गया है। जब समाधान की सटीकता प्राथमिक महत्व की नहीं है, मोंटे कार्लो विधि | स्टोकेस्टिक कण (मोंटे कार्लो) विधियां एक आकर्षक विकल्प हैं।

संक्षेपण-संचालित एकत्रीकरण
एकत्रीकरण के अलावा, कण संघनन, निक्षेपण या द्वारा आकार में भी बढ़ सकते हैं अभिवृद्धि। हसन और हसन ने हाल ही में एक संक्षेपण-संचालित एकत्रीकरण (सीडीए) मॉडल का प्रस्ताव किया है जिसमें एकत्रीकरण कण टकराव पर विलय के बीच लगातार बढ़ते रहते हैं। सीडीए मॉडल को निम्नलिखित प्रतिक्रिया योजना द्वारा समझा जा सकता है


 * $$A_x(t) + A_y(t) \stackrel{

v(x,t)}{\longrightarrow} A_{(\alpha + 1)(x + y)}(t + \tau),$$ कहाँ $$A_x(t)$$ आकार के योग को दर्शाता है $$x$$ समय पर $$t$$ और $$\tau$$ बीता हुआ समय है। इस प्रतिक्रिया योजना को निम्नलिखित सामान्यीकृत स्मोलुचोव्स्की समीकरण द्वारा वर्णित किया जा सकता है


 * $$\Big[{{\partial }\over{\partial t}} + {{\partial}\over{\partial x}} v(x,t) \Big]n(x,t)

=-n(x,t)\int_0^\infty K(x,y)n(y,t)dy + {{1}\over{2}}\int_0^x dy K(y,x-y) n(y,t)n(x-y,t).$$ आकार का एक कण मानते हुए $$x$$ टक्कर के समय के बीच संघनन के कारण बढ़ता है $$\tau(x)$$ के व्युत्क्रम के बराबर $$\int_0^\infty K(x,y)n(y,t)dy$$ एक राशि से $$\alpha x$$ अर्थात।


 * $$v(x,t)={{\alpha x}\over{\tau(x)}}=\alpha x\int_0^\infty dyK(x,y)n(y,t).$$

निरंतर कर्नेल देने के लिए कोई भी सामान्यीकृत स्मोलुचोव्स्की समीकरण को हल कर सकता है


 * $$n(x,t)\sim t^{-(2+2\alpha)}e^{-{{x}\over{t^{1+2\alpha}}}},$$

जो गतिशील स्केलिंग प्रदर्शित करता है। एक साधारण भग्न विश्लेषण से पता चलता है कि संक्षेपण-संचालित एकत्रीकरण को आयाम के भग्न का सबसे अच्छा वर्णन किया जा सकता है


 * $$d_f={{1}\over{1+2\alpha}}.$$

$$d_f$$वें> वें क्षण $$n(x,t)$$ हमेशा एक संरक्षित मात्रा होती है जो डायनेमिक स्केलिंग के सभी घातांकों को ठीक करने के लिए जिम्मेदार होती है। ऐसा संरक्षण नियम कैंटर सेट में भी पाया गया है।

यह भी देखें

 * आइंस्टीन संबंध (गतिज सिद्धांत)|आइंस्टीन-स्मोलुचोव्स्की संबंध
 * फ्लोकुलेशन
 * स्मोलुचोव्स्की कारक
 * विलियम्स स्प्रे समीकरण