Bir kristalizasyon süreci laboratuvardan üretim ölçeğine taşındığında karıştırma, sıcaklık ve türbülans artık tankın her noktasında aynı kalmaz. Duvar yakınındaki daha soğuk bölgeler, karıştırıcının çevresindeki yüksek türbülans ve düşük akışlı durgun bölgeler farklı çekirdeklenme ve kristal büyüme koşulları oluşturabilir. Bunun sonucunda aynı kimyasal madde, tankın farklı bölgelerinde farklı boyut ve şekillerde kristalleşebilir.

Araştırmacılar bu sorunu, hesaplamalı akışkanlar dinamiğinden elde edilen bilgileri iki boyutlu popülasyon denge modeliyle birleştiren fizik esinli çok bölmeli bir çerçeveyle ele aldı. On litrelik, çift karıştırıcılı kristalizatör yedi hidrodinamik bölgeye ayrıldı. Her bölge için sıcaklık, akış değişimi, türbülanslı enerji yayılımı, aşırı doygunluk, çekirdeklenme ve kristallerin uzunluk-genişlik yönlerindeki büyümesi ayrı ayrı hesaplandı.

Model, standart momentler yönteminin hesaplama hızını yüksek çözünürlüklü sonlu hacim yönteminin ayrıntılı kristal boyutu dağılımıyla birleştiriyor. Tek bölmeli 400 × 400 çözünürlüklü karşılaştırmada hesaplama süresi 262 saniyeden 36,2 saniyeye düştü. Yedi bölmeli 100 × 100 çözünürlüklü karşılaştırmada ise tam sonlu hacim çözümü 1.245,6 saniye sürerken hibrit model aynı işlemi 53,1 saniyede tamamladı; bu yaklaşık 23,5 kat hızlanma anlamına geliyor.

Sonuçlar, çözeltinin ortalama derişimi tank genelinde neredeyse eşit kalsa bile duvar kaynaklı sıcaklık farklarının yerel aşırı doygunluğu ve büyüme hızını değiştirebildiğini gösterdi. Türbülanslı enerji yayılımı çekirdeklenme denklemine dahil edildiğinde karıştırıcıya yakın bölgeler baskın çekirdeklenme merkezlerine dönüştü. Bu farklılıklar nihai kristal boyutu dağılımını ve en-boy oranını 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üğü, boyut dağılımı, şekli, en-boy oranı 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, akışkanlı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 boyutlu modeller, uzunluk 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 türbülanslı enerji yayılımı yükselirken uzak bölgelerde azalabilir.
• Karıştırma süresi ile kristal büyüme 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, hesaplamalı akışkanlar dinamiği ile popülasyon denge denklemlerinin tankın her hücresinde eş zamanlı çözülmesidir. Ancak çok saatlik bir kristalizasyon sürecinde iki iç boyuta 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 reaktörün hızlı hidrodinamik ve kısa süreli ısıl davranışını belirlemek için kullanılıyor. Daha sonra bu bilgiler, uzun kristalizasyon sürecini çok daha düşük maliyetle hesaplayan bölmeli modele aktarılıyor.

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

1. Reaktörün hız ve türbülans 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 boyutlu popülasyon denge modeli.

On litrelik reaktör yedi bölgeye ayrıldı

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

CFD sonuçlarına göre reaktör ş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ölmenin içindeki hız dağılımını ve türbülanslı enerji yayılımını mümkün olduğunca homojen hale getirmeye çalıştı.

İki amaç fonksiyonu, hız ile enerji yayılımı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ü, ε türbülanslı kinetik enerji yayılım hızını, d ise bölme geometrisini belirleyen boyutları gösteriyor. Sayfa 27’deki yedi Pareto grafiğinde, her bölme için hız homojenliği ile enerji yayılımı 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 karıştırıcının bulunduğu dördüncü bölmede görüldü. Bölme sınırları karıştı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 yayılımı 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 kaymama koşulu ve duvar fonksiyonları, sıvı-hava arayüzünde ise kayma 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 duvar süresi gerektirdi.

Uzun süreli ısı transferi nasıl hızlandırıldı?

Dört saatlik veya daha uzun bir kristalizasyon 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 soğutma ve ısıtma koşulları örneklenerek düşük mertebeli bir ısıl model oluşturuldu.

Her bölmenin 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ölme hacmini, F_j→i bölmeler arasındaki hacimsel akışı, T_i bölme sıcaklığını, T_w duvar sıcaklığını ve h etkin ısı transfer katsayısını gösteriyor.

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

• Duvara komşu 1-3 numaralı bölmeler,
• Merkez ve karıştırıcı etkisindeki 4-7 numaralı bölmeler.

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 boyutta izlendi

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

Basitleştirilmiş iki boyutlu popülasyon denge 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 kristallerin uzunluk ve genişlik yönündeki büyümesini, Dirac delta terimi yeni çekirdeklerin belirli bir başlangıç boyutunda oluşmasını, son terim ise kristallerin bölmeler arasındaki taşınmasını ifade ediyor.

Başlangıç tohumu iki boyutlu Gauss dağılımıyla tanımlandı. 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 kırılma, aglomerasyon ve çözünmeyi ihmal ediyor. Bu nedenle sonuçlar büyüme ve ikincil çekirdeklenmenin baskın olduğu kristalizasyon koşullarıyla sınırlıdır.

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

Her bölmede bütün uzunluk-genişlik 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 boyutu ve toplam kristal hacmini düşük sayıda diferansiyel denklemle temsil ediyor. Böylece temel kristalizasyon dinamikleri standart momentler yöntemiyle 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 boyutlu kristal boyutu dağılımını yüksek çözünürlüklü sonlu hacim yöntemiyle yeniden kurdu.

Yüksek çözünürlüklü sonlu hacim yöntemi

HRFVM yönteminde uzunluk ve genişlik alanı düzenli bir hesaplama ağına ayrıldı. Kristaller büyüdükçe sayı yoğunluğu küçük boyut hücrelerinden daha büyük boyut 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üme ve çekirdeklenme

Her bölmenin aşırı doygunluğu yerel derişim ile sıcaklığa bağlı çözünürlükten 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 uzunluk ve genişlik yönlerindeki büyüme 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ümemesine ve en-boy oranının süreç boyunca değişmesine izin veriyor.

İkincil çekirdeklenme ise aşırı doygunluk, mevcut kristal hacmi ve yerel enerji yayılımıyla ilişkilendirildi:

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

Enerji yayılımı etkisi etkinleştirildiğinde karıştırıcı çevresindeki yüksek türbülanslı bölgelerde kristal-kristal çarpışmaları ve yüzey aşınmasıyla ilişkili çekirdeklenme kuvvetlenebiliyor.

Hesaplama performansı

İlk karşılaştırma, tek bölmeli ve tamamen karışmış bir kristalizasyon modeli üzerinde yapıldı. Tam HRFVM ile hibrit SMOM-HRFVM aynı ağ çözünürlüklerinde 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ölmeli hesaplamada fark daha belirgin hale geldi:

Araştırmacılar genel hız kazanımını problem boyutuna bağlı olarak bir ile iki mertebe arasında tanımlıyor. Tablodaki doğrudan karşılaştırılabilen en güçlü sonuç yaklaşık 23,5 kat hızlanmadır.

Sıcaklık farkı derişim eşitken bile büyümeyi değiştirdi

Uzaysal heterojenlik kademeli biçimde modele eklendi. İlk durumda bölmeli yapı ideal karışım davranışına yakındı. İkinci durumda duvar kaynaklı sıcaklık farkları, sonraki durumlarda ise enerji yayılımına bağlı çekirdeklenme etkileri eklendi.

Sayfa 30’daki grafiklerde bütün bölmelerin çözünen madde derişimi eğrileri neredeyse üst üste geliyor. Bunun nedeni, incelenen temel koşullarda bölmeler arası karışımın kristalizasyon kinetiğinden daha hızlı olmasıdır.

Buna karşılık aynı grafikteki aşırı doygunluk ve genişlik yönündeki büyüme hızı eğrileri iki gruba ayrılıyor. Duvara yakın 1-3 numaralı bölmeler daha düşük sıcaklık nedeniyle daha yüksek aşırı doygunluğa ve daha yüksek büyüme hızına ulaşıyor. Merkezdeki 4-7 numaralı bölmeler 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 derişim ölçümü homojen görünse bile yerel sıcaklık farkları kristal büyümesini heterojen hale getirebilir.

Türbülans çekirdeklenmenin yerini belirledi

Enerji yayılımının çekirdeklenme üzerinde etkili olmadığı ilk durumlarda bölmelerin çekirdeklenme 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 çekirdeklenme grafikleri, karıştırıcıya yakın ve yüksek enerji yayılımına sahip bölmelerin baskın çekirdeklenme merkezleri haline geldiğini gösteriyor. Daha düşük türbülanslı bölgelerde ise çekirdeklenme büyük ölçüde bastırıldı.

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

• Yüksek türbülanslı bölgelerde daha fazla yeni çekirdek oluşarak küçük kristal nüfusu artıyor.
• Düşük çekirdeklenmeli bölgelerde mevcut kristaller büyümeye devam ederek daha büyük parçacıklar oluşturuyor.

Kristal boyutu ve en-boy oranı nasıl değişti?

Nominal kinetik ve 0,0111 °C/s soğutma hızında, sıcaklık heterojenliği tek başına nihai en-boy oranını çok az değiştirdi. Enerji yayılımına bağlı çekirdeklenme kuvvetlendikçe ortalama kristal boyutları belirgin biçimde büyüdü.

Birinci ve beşinci durum karşılaştırıldığında hacim ağırlıklı genişlik yaklaşık yüzde 24,5, uzunluk 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 boyutlu dağılımlar, ortalama en-boy oranı benzer görünse bile parçacık yoğunluğunun uzunluk-genişlik düzlemindeki konumunun ve dağılım biçiminin değiştiğini gösteriyor. Bu nedenle yalnızca ortalama en-boy oranına bakmak, ürün dağılımındaki bütün farkları ortaya koymuyor.

Yavaş soğutma daha uzun kristaller oluşturdu

Soğutma hızı 0,0111 °C/s’den 0,0011 °C/s’ye düşürüldüğünde kristallerin büyüme için daha fazla zamanı oldu. Uzunluk yönündeki büyüme genişlik 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ş soğutmada hacim ağırlıklı genişlik 269,7 µm, uzunluk 785,1 µm ve en-boy oranı 2,91 olarak hesaplandı.

Bu bulgu, soğutma programının yalnızca ortalama kristal boyutunu değil, büyümenin iki eksen arasındaki dağılımını da değiştirebildiğini gösteriyor.

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

Araştırmacılar kristal büyüme sabitlerini 90 kat artırarak kristalizasyon 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ölmeler arası akış, çözünen madde derişimini 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 derişim eğrisi yerine yedi ayrı bölme eğrisi oluştu.

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

Hızlandırılmış büyümede çekirdeklenme hızları bölmeler arasında büyük farklar gösterse de süreç büyüme baskın kaldı. Bu nedenle hacim ağırlıklı nihai kristal dağılımları birbirine yakın kaldı. Sayfa 35’teki dağılımlarda en-boy oranı 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, kristalizasyon ölçek büyütmede yalnızca tankın ortalama sıcaklığına, karıştırıcı devrine veya ortalama aşırı doygunluğa bakmanın yeterli olmayabileceğini gösteriyor. Aynı ortalama derişimin arkasında farklı sıcaklık, büyüme ve çekirdeklenme geçmişleri bulunabilir.

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

• Laboratuvar ve üretim reaktörleri arasındaki morfoloji değişiminin öngörülmesi,
• Karıştırıcı konumu ve hızının kristal boyutuna etkisinin araştırılması,
• Duvar soğutma programlarının ve termal çevrimlerin tasarlanması,
• İğnemsi veya levhamsı kristallerin en-boy oranının kontrol edilmesi,
• Çok sayıda işletme koşulunun dijital ortamda hızlı biçimde taranması,
• Hesaplama maliyeti düşük kristalizasyon 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 çözünürlük, büyüme, çekirdeklenme ve enerji yayılımı 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ş çok bölmeli iki boyutlu popülasyon denge 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ölmeli modellerde çözünürlüğe bağlı olarak yaklaşık 3,2-7,2 kat hızlanma sağladı.
• Yedi bölmeli 100 × 100 ağda hesaplama süresi 1.245,6 saniyeden 53,1 saniyeye indi.
• 200 × 200 ve 300 × 300 çok bölmeli 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ölmeler arası derişim farkları ihmal edilebilir düzeyde kaldı.
• Duvar sıcaklığı farkları derişim aynı kalırken aşırı doygunluk ve büyüme hızını bölgesel olarak değiştirdi.
• Enerji yayılımına bağlı çekirdeklenme, karıştırıcı çevresindeki bölgeleri baskın çekirdeklenme merkezlerine dönüştürdü.
• Nominal hızlı soğutmada en-boy oranı 2,23-2,25 aralığında kaldı.
• Yavaş soğutmada en-boy oranı yaklaşık 2,91-2,95 aralığına yükseldi.
• Büyüme sabitleri 90 kat artırıldığında derişim alanı yedi ayrı bölgesel eğriye ayrıldı.
• Hızlandırılmış büyümede güçlü bölgesel çekirdeklenme farklarına rağmen nihai hacim dağılımı büyüme 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ğ çözünürlüklerindeki 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 reaktörden alınmış hız, sıcaklık, kristal boyutu dağılımı veya en-boy oranı ö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 boyutu ve en-boy oranı sonuçları fiziksel deneyle doğrulanmamıştır.
• CFD hız ve türbülans 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ölmenin kendi içinde sıcaklık, derişim ve enerji yayılımı bakımından homojen olduğu kabul edilmiştir.
• Akış alanı kristalizasyon 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ölmeyi iki sıcaklık bölgesine indirgemektedir.
• Çekirdeklenmenin enerji yayılımına bağımlılığı ampirik bir terimle temsil edilmiştir.
• Kinetik parametreler başka maddelere veya reaktörlere 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 derişim, aşırı doygunluk ve büyüme grafikleriyle uyumlu görünmemektedir; bu durum preprintte bir şekil 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 boyutu sonuçları, şekil 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.