Bir kristallaşma süreci laboratuvardan üretim ölçeğine taşındığında qarışdırma, sıcaklık ve turbulentlik artık tankın her noktasında aynı kalmaz. Duvar yakınındaki daha soğuk bölgeler, qarışdırıcının çevresindeki yüksek turbulentlik ve düşük akışlı durgun bölgeler farklı nüvələnmə ve kristal böyümə koşulları oluşturabilir. Bunun sonucunda aynı kimyasal madde, tankın farklı bölgelerinde farklı ölçü ve formalerde kristalleşebilir.

Araştırmacılar bu sorunu, hesablama maye dinamikasınden elde edilen bilgileri ikiölçülü populyasiya balansı modeliyle birleştiren fizik esinli çox bölməli bir çerçeveyle ele aldı. On litrelik, çift qarışdırıcılı kristalizatör yedi hidrodinamik bölgeye ayrıldı. Her bölge için sıcaklık, akış değişimi, turbulentliklı enerji dissipasiyaı, həddən artıq doymuşluq, nüvələnmə ve kristalların uzunluq-en yönlerindeki böyüməsi ayrı ayrı hesaplandı.

Model, standart momentlər üsulunin hesaplama hızını yüksek həllolma qabiliyyətilü sonlu həcm yönteminin ayrıntılı kristal ölçüu dağılımıyla birleştiriyor. Tek bölməli 400 × 400 həllolma qabiliyyətilü karşılaştırmada hesaplama süresi 262 saniyeden 36,2 saniyeye düştü. Yedi bölməli 100 × 100 həllolma qabiliyyətilü karşılaştırmada ise tam sonlu həcm çözümü 1.245,6 saniye sürerken hibrit model aynı işlemi 53,1 saniyede tamamladı; bu yaklaşık 23,5 kat sürətlənmə anlamına geliyor.

Sonuçlar, çözeltinin ortalama konsentrasiyai tank genelinde neredeyse eşit kalsa bile divar kaynaklı sıcaklık farklarının yerel aşırı doygunluğu ve böyümə hızını değiştirebildiğini gösterdi. Türbülanslı enerji dissipasiyaı nüvələnmə denklemine dahil edildiğinde qarışdırıcıya yakın bölgeler baskın nüvələnmə merkezlerine dönüştü. Bu farklılıklar nihai kristal ölçüu dağılımını ve tərəflər nisbətinı belirgin biçimde değiştirdi.

Kristalizasyon ölçek büyütme neden zordur?

Kristalizasyon; ilaç, tarım kimyasalları ve özel kimyasalların ayrıştırılması ve saflaştırılmasında kullanılan temel işlemlerden biridir. Ürün yalnızca kimyasal saflığıyla değerlendirilmez. Kristallerin büyüklüğü, ölçü dağılımı, şekli, tərəflər nisbəti ve kristal formu da sonraki filtrasyon, yıkama, kurutma, taşıma ve formülasyon aşamalarını etkiler.

İğne veya levha biçimli anizotropik kristaller özellikle önemlidir. Uzun ve ince kristaller filtre yataklarında yüksek direnç oluşturabilir, mayelık sorunlarına yol açabilir veya mekanik işlemler sırasında kırılabilir. Kristalin yalnızca tek bir karakteristik çapla temsil edildiği bir ölçülu modeller, uzunluq ve genişliğin farklı hızlarda büyüdüğü bu tür sistemleri yeterince açıklayamaz.

Laboratuvar ölçeğinde karışım ideal duruma yaklaşabilir. Tank büyüdükçe ise:

• Duvar ile merkez arasında sıcaklık farkları oluşabilir.
• Karıştırıcı çevresinde turbulentliklı enerji dissipasiyaı yükselirken uzak bölgelerde azalabilir.
• Karıştırma süresi ile kristal böyümə süresi birbirine yaklaşabilir.
• Çekirdeklenme yalnızca ortalama aşırı doygunluğa değil, yerel çarpışma ve aşınma koşullarına bağlı hale gelebilir.
• Tank genelindeki ortalama değerler, ürünün oluştuğu yerel koşulları gizleyebilir.

En ayrıntılı çözüm, hesablama maye dinamikası ile populyasiya balansı denklemlerinin tankın her hücresinde eş zamanlı çözülmesidir. Ancak çok saatlik bir kristallaşma sürecinde iki iç ölçüa sahip kristal dağılımının milyonlarca akış hücresiyle birlikte hesaplanması, tek bir işletme koşulu için bile saatler veya günler sürebilir.

Araştırmanın temel yaklaşımı

Önerilen PIMC-PBM yaklaşımı, tamamen karışmış tek bölgeli model ile tam CFD-PBE çözümü arasında bir ara düzey oluşturuyor. CFD yalnızca reaktorün hızlı hidrodinamik ve kısa süreli ısıl davranışını belirlemek için kullanılıyor. Daha sonra bu bilgiler, uzun kristallaşma sürecini çok daha düşük maliyetle hesaplayan bölməli modele aktarılıyor.

Çerçeve dört temel bileşenden oluşuyor:

1. Reaktörün hız ve turbulentlik alanını hesaplayan CFD modeli,
2. Reaktörü fiziksel olarak anlamlı bölgelere ayıran optimizasyon yöntemi,
3. Her bölgedeki sıcaklığı hesaplayan düşük mertebeli ısıl vekil model,
4. Kristal uzunluğu ve genişliğini birlikte izleyen ikiölçülü populyasiya balansı modeli.

On litrelik reaktor yedi bölgeye ayrıldı

Çalışmada çift kanat düzeneğine sahip, 10 litre çalışma hacimli bir qarışdırmalı kesikli kristalizatör modellendi. Yaklaşık 1,5 milyon hücreden oluşan O-grid ağı kullanıldı. Sayfa 26’daki reaktor geometrisi ve ağ görseli, alt ve üst qarışdırıcı kanatlarını, sıvı-hava arayüzünü ve divar yakınındaki sıklaştırılmış hesaplama ağını gösteriyor.

CFD sonuçlarına göre reaktor şu yedi bölgeye ayrıldı:

Bölme sınırları elle veya yalnızca görsel değerlendirmeyle belirlenmedi. Araştırmacılar NSGA-II çok amaçlı optimizasyon algoritmasını kullanarak her bölmənin içindeki hız dağılımını ve turbulentliklı enerji dissipasiyaını mümkün olduğunca homojen hale getirmeye çalıştı.

İki amaç fonksiyonu, hız ile enerji dissipasiyaının çeyrekler açıklığına dayanıyordu:

\[ f_1(\mathbf{d})=\operatorname{IQR}(u) \]

\[ f_2(\mathbf{d})=\operatorname{IQR}(\varepsilon) \]

Burada u hız büyüklüğünü, ε turbulentliklı kinetik enerji yayılım hızını, d ise bölmə geometrisini belirleyen ölçüları gösteriyor. Sayfa 27’deki yedi Pareto grafiğinde, her bölmə için hız homojenliği ile enerji dissipasiyaı homojenliği arasındaki ödünleşim ve seçilen çözüm kırmızı yıldızlarla gösteriliyor.

En yüksek hassasiyet alt qarışdırıcının bulunduğu dördüncü bölməde görüldü. Bölme sınırları qarışdırıcının oluşturduğu yeniden dolaşım yapılarını daha iyi kapsayacak biçimde değiştirildiğinde, bölge içindeki enerji dissipasiyaı dağılımı belirgin biçimde daraldı.

CFD modeli nasıl kuruldu?

Akış hesaplamaları OpenFOAM ortamında gerçekleştirildi. Akışkan sıkıştırılamaz ve Newtonyen kabul edildi. Karıştırıcı dönüşü, kanatların zamana bağlı hareketini doğrudan çözmek yerine Çoklu Referans Çerçevesi yöntemiyle temsil edildi.

Dönen referans çerçevesindeki momentum denklemi çalışmada şu biçimde verildi:

\[ \nabla\cdot(\rho\mathbf{u}_r\mathbf{u}) +\rho(\boldsymbol{\omega}\times\mathbf{u}) = -\nabla p+\rho\mathbf{g}+\mu\nabla^2\mathbf{u} \]

Burada u mutlak hızı, u_r dönen çerçeveye göre bağıl hızı, ω açısal hızı, p basıncı ve μ dinamik viskoziteyi temsil ediyor.

Türbülans standart k-ε RANS modeliyle hesaplandı. Duvarlarda sürüşməma koşulu ve divar fonksiyonları, sıvı-hava arayüzünde ise sürüşmə koşulu kullanıldı. Bütün değişkenler için ikinci dereceden uzaysal ayrıklaştırma uygulanarak sayısal yayınımın azaltılması amaçlandı.

Akış alanı hesaplamaları 64 işlemci çekirdeğinde yürütüldü ve her CFD koşulu yaklaşık 20 dakika sürdü. Sıcaklığın zamana bağlı CFD örneklemeleri ise yaklaşık 24 saatlik divar süresi gerektirdi.

Uzun süreli istilik ötürülməsi nasıl hızlandırıldı?

Dört saatlik veya daha uzun bir kristallaşma sürecinin tamamını CFD ile çözmek yerine, üç dakikalık kısa ısıl CFD örnekleri kullanıldı. Duvar sıcaklığının 273 K ve 373 K olduğu soyutma ve ısıtma koşulları örneklenerek düşük mertebeli bir ısıl model oluşturuldu.

Her bölmənin enerji dengesi şu denklemle ifade edildi:

\[ \rho C_pV_i\frac{dT_i}{dt} = \sum_{j\neq i}F_{j\rightarrow i}C_p(T_j-T_i) + hA_i(T_w-T_i) \]

V_i bölmə hacmini, F_j→i bölmələr arasındaki hacimsel akışı, T_i bölmə sıcaklığını, T_w divar sıcaklığını ve h etkin ısı transfer katsayısını gösteriyor.

Yedi hidrodinamik bölmə, ısıl hesaplamayı daha da azaltmak için iki ana grupta toplandı:

• Duvara komşu 1-3 numaralı bölmələr,
• Merkez ve qarışdırıcı etkisindeki 4-7 numaralı bölmələr.

Etkin ısı transfer katsayısı sıcaklığın birinci dereceden fonksiyonu olarak modellendi:

\[ h=aT_i+b \]

Parametreler önce genetik algoritmayla küresel olarak arandı, ardından doğrusal olmayan en küçük kareler yöntemiyle iyileştirildi. Güven aralıkları Levenberg-Marquardt yaklaşımıyla hesaplandı.

Kristaller iki ölçüta izlendi

Model kristalleri tek bir eşdeğer çapla değil, uzunluq l ve en w koordinatlarıyla temsil ediyor. Her bölmədeki sayı yoğunluğu n_z(l,w,t) ile gösteriliyor.

Basitleştirilmiş ikiölçülü populyasiya balansı denklemi şu yapıdadır:

\[ \frac{\partial n_z}{\partial t} + \frac{\partial(n_zG_{l,z})}{\partial l} + \frac{\partial(n_zG_{w,z})}{\partial w} = \delta(l-l_n)\delta(w-w_n)B_z + \sum_{k\neq z} \frac{F_{k\rightarrow z}}{V_z} (n_k-n_z) \]

Denklemin ilk terimi zamana bağlı değişimi, ikinci ve üçüncü terimler kristalların uzunluq ve en yönündeki böyüməsini, Dirac delta terimi yeni çekirdeklerin belirli bir başlangıç ölçüunda oluşmasını, son terim ise kristalların bölmələr arasındaki taşınmasını ifade ediyor.

Başlangıç tohumu ikiölçülü Gauss dağılımıyla diaqnostikamlandı. Ortalama kristal uzunluğu 200 µm, ortalama genişliği 100 µm; iki yöndeki dağılım parametreleri ise 150 µm olarak seçildi.

Model qırılma, aglomerasyon ve çözünmeyi ihmal ediyor. Bu nedenle sonuçlar böyümə ve ikincil nüvələnmənin baskın olduğu kristallaşma koşullarıyla sınırlıdır.

Momentler yöntemi hesaplamayı nasıl azalttı?

Her bölməde bütün uzunluq-en hücrelerini zaman boyunca çözmek yerine önce sayı dağılımının çapraz momentleri hesaplandı:

\[ \mu_{ij,z}(t) = \int_0^\infty\int_0^\infty w^il^j n_z(l,w,t)\,dl\,dw \]

Düşük mertebeli momentler kristal sayısını, ortalama ölçüu ve toplam kristal hacmini düşük sayıda diferansiyel denklemle temsil ediyor. Böylece temel kristallaşma dinamikleri standart momentlər üsuluyle hızla çözülebiliyor.

Ancak yalnızca momentler kullanıldığında dağılımın tamamı, çok tepeli yapı veya kuyruk davranışı kaybedilebilir. Araştırmacılar bu nedenle süreç dinamiklerini SMOM ile çözüp son ikiölçülü kristal ölçüu dağılımını yüksek həllolma qabiliyyətilü sonlu həcm yöntemiyle yeniden kurdu.

Yüksek həllolma qabiliyyətilü sonlu həcm yöntemi

HRFVM yönteminde uzunluq ve en alanı düzenli bir hesaplama ağına ayrıldı. Kristaller büyüdükçe sayı yoğunluğu küçük ölçü hücrelerinden daha büyük ölçü hücrelerine taşındı.

Birinci dereceden yukarı akış şeması kararlılığı sağlarken ikinci dereceden düzeltme sayısal yayınımı azalttı. Keskin dağılım cephelerinde fiziksel olmayan salınımları önlemek için akı sınırlayıcı kullanıldı:

\[ \phi(\theta)=\frac{|\theta|+\theta}{\theta+1} \]

Düzgün bölgelerde yöntem ikinci dereceden doğruluğa yaklaşırken keskin gradyanlarda daha kararlı birinci dereceden davranışa dönüyor.

Aşırı doygunluk, böyümə ve nüvələnmə

Her bölmənin aşırı doygunluğu yerel konsentrasiya ile sıcaklığa bağlı həllolma qabiliyyətiten hesaplandı:

\[ \sigma_z=\frac{C_z-C_s(T_z)}{C_s(T_z)} \]

Çözünürlük ikinci dereceden sıcaklık polinomuyla temsil edildi:

\[ C_s=a+bT_z+cT_z^2 \]

Kristalin uzunluq ve en yönlerindeki böyümə hızları ayrı bağıntılarla verildi:

\[ G_{l,z}=k_{gl}\sigma_z^{g_l} \]

\[ G_{w,z}=k_{gw}\sigma_z^{g_w} \]

Bu ayrım, kristalin iki yönde aynı hızla böyüməmesine ve tərəflər nisbətinın süreç boyunca değişmesine izin veriyor.

İkincil nüvələnmə ise həddən artıq doymuşluq, mevcut kristal hacmi ve yerel enerji dissipasiyaıyla ilişkilendirildi:

\[ B_z=k_s\sigma_z^sM_{t,z}^{\,j}D_{\varepsilon,z} \]

Enerji yayılımı etkisi etkinleştirildiğinde qarışdırıcı çevresindeki yüksek turbulentliklı bölgelerde kristal-kristal çarpışmaları ve yüzey aşınmasıyla ilişkili nüvələnmə kuvvetlenebiliyor.

Hesaplama performansı

İlk karşılaştırma, tek bölməli ve tamamen karışmış bir kristallaşma modeli üzerinde yapıldı. Tam HRFVM ile hibrit SMOM-HRFVM aynı ağ həllolma qabiliyyətilerinde karşılaştırıldı.

Raporlanan moment farkları çoğunlukla yüzde 0,514’ün altında kaldı. Böylece hız kazancı, temel momentlerin doğruluğunda büyük bir kayıp olmadan elde edildi.

Yedi bölməli hesaplamada fark daha belirgin hale geldi:

Araştırmacılar genel hız kazanımını problem ölçüuna bağlı olarak bir ile iki mertebe arasında diaqnostikamlıyor. Tablodaki doğrudan karşılaştırılabilen en güçlü sonuç yaklaşık 23,5 kat sürətlənmədır.

Sıcaklık farkı konsentrasiya eşitken bile böyüməyi değiştirdi

Uzaysal heterojenlik kademeli biçimde modele eklendi. İlk durumda bölməli yapı ideal karışım davranışına yakındı. İkinci durumda divar kaynaklı sıcaklık farkları, sonraki durumlarda ise enerji dissipasiyaına bağlı nüvələnmə etkileri eklendi.

Sayfa 30’daki grafiklerde bütün bölmələrin çözünen madde konsentrasiyai eğrileri neredeyse üst üste geliyor. Bunun nedeni, incelenen temel koşullarda bölmələr arası karışımın kristallaşma kinetiğinden daha hızlı olmasıdır.

Buna karşılık aynı grafikteki həddən artıq doymuşluq ve en yönündeki böyümə hızı eğrileri iki gruba ayrılıyor. Duvara yakın 1-3 numaralı bölmələr daha düşük sıcaklık nedeniyle daha yüksek aşırı doygunluğa ve daha yüksek böyümə hızına ulaşıyor. Merkezdeki 4-7 numaralı bölmələr ise daha sıcak kalıyor.

Bu sonuç önemli bir ölçek büyütme mekanizmasına işaret ediyor: tankın tamamından alınan ortalama konsentrasiya ölçümü homojen görünse bile yerel sıcaklık farkları kristal böyüməsini heterojen hale getirebilir.

Türbülans nüvələnmənin yerini belirledi

Enerji yayılımının nüvələnmə üzerinde etkili olmadığı ilk durumlarda bölmələrin nüvələnmə eğrileri birbirine yakındı. Enerji yayılımına bağımlılık kuvvetlendirildikçe eğriler belirgin biçimde ayrıldı.

Sayfa 31’deki nüvələnmə grafikleri, qarışdırıcıya yakın ve yüksek enerji dissipasiyaına sahip bölmələrin baskın nüvələnmə merkezleri haline geldiğini gösteriyor. Daha düşük turbulentliklı bölgelerde ise nüvələnmə büyük ölçüde bastırıldı.

Bu mekanizma nihai ürün dağılımını iki karşıt yolla etkiliyor:

• Yüksek turbulentliklı bölgelerde daha fazla yeni çekirdek oluşarak küçük kristal nüfusu artıyor.
• Düşük nüvələnməli bölgelerde mevcut kristaller böyüməye devam ederek daha büyük parçacıklar oluşturuyor.

Kristal ölçüu ve tərəflər nisbəti nasıl değişti?

Nominal kinetik ve 0,0111 °C/s soyutma hızında, sıcaklık heterojenliği tek başına nihai tərəflər nisbətinı çok az değiştirdi. Enerji yayılımına bağlı nüvələnmə kuvvetlendikçe ortalama kristal ölçüları belirgin biçimde büyüdü.

Birinci ve beşinci durum karşılaştırıldığında hacim ağırlıklı en yaklaşık yüzde 24,5, uzunluq ise yaklaşık yüzde 25,5 arttı. En-boy oranındaki değişim daha sınırlı kalarak 2,23’ten 2,25’e yükseldi.

Sayfa 32’deki ikiölçülü dağılımlar, ortalama tərəflər nisbəti benzer görünse bile parçacık yoğunluğunun uzunluq-en düzlemindeki konumunun ve dağılım biçiminin değiştiğini gösteriyor. Bu nedenle yalnızca ortalama tərəflər nisbətina bakmak, ürün dağılımındaki bütün farkları ortaya koymuyor.

Yavaş soyutma daha uzun kristaller oluşturdu

Soğutma hızı 0,0111 °C/s’den 0,0011 °C/s’ye düşürüldüğünde kristalların böyümə için daha fazla zamanı oldu. Uzunluk yönündeki böyümə en yönünden daha güçlü olduğu için en-boy oranları yaklaşık 2,91-2,95 aralığına çıktı.

En güçlü hidrodinamik heterojenliğe sahip beşinci durumda, yavaş soyutmada hacim ağırlıklı en 269,7 µm, uzunluq 785,1 µm ve tərəflər nisbəti 2,91 olarak hesaplandı.

Bu bulgu, soyutma programının yalnızca ortalama kristal ölçüunu değil, böyümənin iki eksen arasındaki dağılımını da değiştirebildiğini gösteriyor.

Karışım ve kristallaşma süreleri yarıştığında ne oldu?

Araştırmacılar kristal böyümə sabitlerini 90 kat artırarak kristallaşma süresini yapay olarak kısalttı. Karıştırma süresi aynı bırakıldığı için iki süreç arasındaki rekabet doğrudan incelenebildi.

Temel kinetikte bölmələr arası akış, çözünen madde konsentrasiyaini hızla eşitliyordu. Hızlandırılmış kinetikte ise kristaller aşırı doygunluğu karışımın eşitleyebileceğinden daha hızlı tüketmeye başladı. Bunun sonucunda tek bir konsentrasiya eğrisi yerine yedi ayrı bölmə eğrisi oluştu.

Duvara yakın daha soğuk bölgelerde böyümə ve çözünen madde tüketimi hızlandı. Merkez bölgeler ise bölmələr arası alışverişle seyreldi. Böylece ölçek büyütmede belirleyici unsurun yalnızca qarışdırma hızı veya yalnızca kristal böyümə hızı olmadığı; bu iki sürenin birbirine oranı olduğu gösterildi.

Hızlandırılmış böyüməde nüvələnmə hızları bölmələr arasında büyük farklar gösterse de süreç böyümə baskın kaldı. Bu nedenle hacim ağırlıklı nihai kristal dağılımları birbirine yakın kaldı. Sayfa 35’teki dağılımlarda tərəflər nisbəti temel durumda 1,44, heterojen durumlarda yaklaşık 1,50 olarak görülüyor.

Mühendislik açısından ne anlama geliyor?

Çalışma, kristallaşma ölçek büyütmede yalnızca tankın ortalama sıcaklığına, qarışdırıcı devrine veya ortalama aşırı doygunluğa bakmanın yeterli olmayabileceğini gösteriyor. Aynı ortalama konsentrasiyain arkasında farklı sıcaklık, böyümə ve nüvələnmə geçmişleri bulunabilir.

Önerilen çerçeve özellikle şu amaçlarla değerlendirilebilir:

• Laboratuvar ve üretim reaktorleri arasındaki morfoloji değişiminin öngörülmesi,
• Karıştırıcı konumu ve hızının kristal ölçüuna etkisinin araştırılması,
• Duvar soyutma programlarının ve termal çevrimlerin tasarlanması,
• İğnemsi veya levhamsı kristalların tərəflər nisbətinın kontrol edilmesi,
• Çok sayıda işletme koşulunun dijital ortamda hızlı biçimde taranması,
• Hesaplama maliyeti düşük kristallaşma dijital ikizleri geliştirilmesi.

Bununla birlikte model, deneysel doğrulama yapılmadan doğrudan üretim kontrolüne bağlanmamalıdır. Her kimyasal bileşik için həllolma qabiliyyəti, böyümə, nüvələnmə ve enerji dissipasiyaı bağımlılıklarının deneysel olarak belirlenmesi gerekir.

Çalışmanın Yöntemi ve Bulguları

Araştırma tasarımı

Çalışma, CFD ile bilgilendirilmiş çox bölməli ikiölçülü populyasiya balansı modelinin geliştirilmesi ve sayısal karşılaştırmalarla değerlendirilmesine dayanan hesaplamalı bir kimya mühendisliği araştırmasıdır.

Temel işletme koşulları

Yazılım ve hesaplama kaynakları

• HRFVM kodu C dilinde geliştirildi ve MATLAB’a MEX olarak bağlandı.
• SMOM denklemleri MATLAB betikleriyle çözüldü.
• PBM karşılaştırmaları Intel Core i9-12900KF, 32 GB DDR5 belleğe sahip masaüstünde tek iş parçacığıyla yürütüldü.
• CFD akış hesaplamaları 64 CPU çekirdeğinde yaklaşık 20 dakika sürdü.
• Geçici sıcaklık CFD hesaplamaları yaklaşık 24 saat sürdü.
• Bölme optimizasyonu tek CPU çekirdeğinde yaklaşık 30 dakika sürdü.
• Bölmeler arası akış çıkarımı 8 CPU çekirdeğinde yaklaşık 5 dakikada tamamlandı.

Başlıca sayısal bulgular

• Hibrit yöntem tek bölməli modellerde çözünürlüğe bağlı olarak yaklaşık 3,2-7,2 kat sürətlənmə sağladı.
• Yedi bölməli 100 × 100 ağda hesaplama süresi 1.245,6 saniyeden 53,1 saniyeye indi.
• 200 × 200 ve 300 × 300 çox bölməli tam HRFVM çözümleri makul sürede tamamlanamazken hibrit yöntem sırasıyla 841,5 ve 1.667,3 saniyede çalıştı.
• Temel kinetik koşullarında bölmələr arası konsentrasiya farkları ihmal edilebilir düzeyde kaldı.
• Duvar sıcaklığı farkları konsentrasiya aynı kalırken həddən artıq doymuşluq ve böyümə hızını bölgesel olarak değiştirdi.
• Enerji yayılımına bağlı nüvələnmə, qarışdırıcı çevresindeki bölgeleri baskın nüvələnmə merkezlerine dönüştürdü.
• Nominal hızlı soyutmada tərəflər nisbəti 2,23-2,25 aralığında kaldı.
• Yavaş soyutmada tərəflər nisbəti yaklaşık 2,91-2,95 aralığına yükseldi.
• Büyüme sabitleri 90 kat artırıldığında konsentrasiya alanı yedi ayrı bölgesel eğriye ayrıldı.
• Hızlandırılmış böyüməde güçlü bölgesel nüvələnmə farklarına rağmen nihai hacim dağılımı böyümə baskınlığı nedeniyle sınırlı değişti.

Doğrulama ve karşılaştırma sınırı

Modelin sayısal doğruluğu, SMOM momentleri ile HRFVM sonuçlarının karşılaştırılması ve farklı ağ həllolma qabiliyyətilerindeki hesaplama süreleri üzerinden değerlendirildi. Reaktör CFD ağı bir ağ yakınsama çalışmasından sonra yaklaşık 1,5 milyon hücre olarak seçildiği belirtilmekle birlikte, ağ yakınsama değerlerinin ayrıntılı tablosu sunulmadı.

Çalışmada fiziksel reaktorden alınmış hız, sıcaklık, kristal ölçüu dağılımı veya tərəflər nisbəti ölçümleriyle doğrudan doğrulama yapılmadı. Bu nedenle sonuçlar model içi sayısal karşılaştırma ve mekanistik tutarlılık düzeyinde değerlendirilmelidir.

Kritik sınırlılıklar

• Çalışma hakem değerlendirmesinden geçmemiş bir preprinttir.
• Nihai kristal ölçüu ve tərəflər nisbəti sonuçları fiziksel deneyle doğrulanmamıştır.
• CFD hız ve turbulentlik alanlarının deneysel ölçümlerle karşılaştırması sunulmamıştır.
• Kırılma, aglomerasyon ve çözünme mekanizmaları modele dahil edilmemiştir.
• Her bölmənin kendi içinde sıcaklık, konsentrasiya ve enerji dissipasiyaı bakımından homojen olduğu kabul edilmiştir.
• Akış alanı kristallaşma boyunca yarı kararlı kabul edilmiştir.
• Sıcaklık değişimlerinin yoğunluk farkı ve doğal taşınım yoluyla akış alanını değiştirmediği varsayılmıştır.
• Isıl model yedi hidrodinamik bölməyi iki sıcaklık bölgesine indirgemektedir.
• Çekirdeklenmenin enerji dissipasiyaına bağımlılığı ampirik bir terimle temsil edilmiştir.
• Kinetik parametreler başka maddelere veya reaktorlere doğrudan aktarılamaz.
• Çalışmada açık veri seti veya kaynak kod deposu belirtilmemiştir.
• Tablo 6 ve Tablo 7’de hata sütun başlıkları aynı moment sembolüyle tekrarlanmıştır; hangi farklı momentlerin karşılaştırıldığı tablo başlıklarından açıkça anlaşılamamaktadır.
• Sayfa 33’teki Şekil 11 açıklaması, görüntülenen konsentrasiya, həddən artıq doymuşluq ve böyümə grafikleriyle uyumlu görünmemektedir; bu durum preprintte bir forma açıklaması veya numaralandırma hatasına işaret etmektedir.

Kaynak ve Yöntem Notu

Özgün çalışma: “A Physics-Inspired Multi-compartmental (PIMC) Population Balance Framework for Size and Aspect Ratio Control in Crystallization Scale-Up.”

Yazarlar: Yung-Shun Kang, Antonello Raponi ve Zoltan K. Nagy.

Sorumlu yazar: Zoltan K. Nagy. Eş katkı beyanı bulunmamaktadır.

Kurumlar: Davidson School of Chemical Engineering, Purdue University, West Lafayette, Indiana, ABD; Mines Saint-Étienne, Univ Lyon, CNRS, UMR 5307 LGF, Centre SPIN, Saint-Étienne, Fransa.

Yazar katkıları: Yung-Shun Kang kavramsallaştırma, yazılım, araştırma, yöntem ve ilk taslağı; Antonello Raponi yazılım, yöntem ve metin değerlendirmesini; Zoltan K. Nagy kavramsallaştırma, danışmanlık, proje yönetimi ve finansman edinimini üstlenmiştir.

Kaynak türü: Özgün hesaplamalı kimya mühendisliği araştırması ve preprint.

Yıl: 2026.

Hakemlik durumu: Çalışma hakem değerlendirmesinden geçmemiştir. Bulgular, hakemli ve nihai dergi sonuçları olarak değerlendirilmemelidir.

Platform: SSRN.

Hakemli dergi: Herhangi bir hakemli dergide kabul veya yayımlanma bilgisi bulunmamaktadır.

DOI:10.2139/ssrn.6948349

Resmî çalışma bağlantısı:SSRN çalışma kaydı

Finansman: Çalışma, ABD Ulusal Bilim Vakfının 2132142 numaralı hibesiyle kısmen desteklenmiştir. Corteva Agriscience’ın kısmi mali desteği ve araştırma iş birliği de belirtilmiştir.

Çıkar çatışması: Çıkar çatışması beyanı metinde yer almamaktadır.

Veri ve kod erişimi: Açık veri seti, model dosyaları veya kaynak kod deposu bağlantısı verilmemiştir.

Bu içerikteki yöntem, denklemler, CFD koşulları, hesaplama süreleri, kristal ölçüu sonuçları, forma ve tablo yorumları, bulgular ve sınırlılıklar incelenen çalışmaya dayanmaktadır. Dış kaynaklar yalnızca başlık, yazar kimliği, kurum, DOI, platform ve resmî çalışma bağlantısının bibliyografik doğrulanması amacıyla kullanılmış; çalışmada bulunmayan bilimsel sonuç eklenmemiştir.