Bir кристаллдашуу süreci laboratuvardan üretim ölçeğine taşındığında аралаштыруу, sıcaklık ve турбуленттүүлүк artık tankın her noktasında aynı kalmaz. Duvar yakınındaki daha soğuk bölgeler, аралаштыргычnın çevresindeki yüksek турбуленттүүлүк ve düşük akışlı durgun bölgeler farklı ядролонуу ve кристалл өсүү koşulları oluşturabilir. Bunun sonucunda aynı kimyasal madde, tankın farklı bölgelerinde farklı өлчөм ve формаlerde кристаллleşebilir.

Araştırmacılar bu sorunu, эсептик суюктук динамикасыnden elde edilen bilgileri эки өлчөмдүү популяция балансы modeliyle birleştiren fizik esinli көп бөлүмдүү bir çerçeveyle ele aldı. On litrelik, çift аралаштыргычlı кристаллizatör yedi hidrodinamik bölgeye ayrıldı. Her bölge için sıcaklık, akış değişimi, турбуленттүүлүкlı энергия диссипациясы, ашыкча каныгуу, ядролонуу ve кристаллдардын узундук-кеңдик yönlerindeki өсүүsi ayrı ayrı hesaplandı.

Model, standart моменттер ыкмасыnin hesaplama hızını yüksek эригичтикlü чектүү көлөм yönteminin ayrıntılı кристалл өлчөмu dağılımıyla birleştiriyor. Tek бөлүмli 400 × 400 эригичтикlü karşılaştırmada hesaplama süresi 262 saniyeden 36,2 saniyeye düştü. Yedi бөлүмli 100 × 100 эригичтикlü karşılaştırmada ise tam чектүү көлөм çözümü 1.245,6 saniye sürerken hibrit model aynı işlemi 53,1 saniyede tamamladı; bu yaklaşık 23,5 kat ылдамдоо anlamına geliyor.

Sonuçlar, çözeltinin ortalama концентрацияi tank genelinde neredeyse eşit kalsa bile дубал kaynaklı sıcaklık farklarının yerel aşırı doygunluğu ve өсүү hızını değiştirebildiğini gösterdi. Türbülanslı энергия диссипациясы ядролонуу denklemine dahil edildiğinde аралаштыргычya yakın bölgeler baskın ядролонуу merkezlerine dönüştü. Bu farklılıklar nihai кристалл өлчөмu dağılımını ve тараптар катышы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üğü, өлчөм dağılımı, şekli, тараптар катышы ve кристалл formu da sonraki filtrasyon, yıkama, kurutma, taşıma ve formülasyon aşamalarını etkiler.

İğne veya levha biçimli anizotropik кристаллler özellikle önemlidir. Uzun ve ince кристаллler filtre yataklarında yüksek direnç oluşturabilir, суюктукlı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 өлчөмlu modeller, узундук 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 турбуленттүүлүкlı энергия диссипациясы yükselirken uzak bölgelerde azalabilir.
• Karıştırma süresi ile кристалл өсүү 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, эсептик суюктук динамикасы ile популяция балансы denklemlerinin tankın her hücresinde eş zamanlı çözülmesidir. Ancak çok saatlik bir кристаллдашуу sürecinde iki iç өлчөмa sahip кристалл 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 реакторün hızlı hidrodinamik ve kısa süreli ısıl davranışını belirlemek için kullanılıyor. Daha sonra bu bilgiler, uzun кристаллдашуу sürecini çok daha düşük maliyetle hesaplayan бөлүмli modele aktarılıyor.

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

1. Reaktörün hız ve турбуленттүүлүк 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 эки өлчөмдүү популяция балансы modeli.

On litrelik реактор yedi bölgeye ayrıldı

Çalışmada çift kanat düzeneğine sahip, 10 litre çalışma hacimli bir аралаштырууlı kesikli кристаллizatör modellendi. Yaklaşık 1,5 milyon hücreden oluşan O-grid ağı kullanıldı. Sayfa 26’daki реактор geometrisi ve ağ görseli, alt ve üst аралаштыргыч kanatlarını, sıvı-hava интерфейсünü ve дубал yakınındaki sıklaştırılmış hesaplama ağını gösteriyor.

CFD sonuçlarına göre реактор ş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 бөлүмnin içindeki hız dağılımını ve турбуленттүүлүкlı энергия диссипациясыnı mümkün olduğunca homojen hale getirmeye çalıştı.

İki amaç fonksiyonu, hız ile энергия диссипациясы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ü, ε турбуленттүүлүкlı kinetik enerji yayılım hızını, d ise бөлүм geometrisini belirleyen өлчөмları gösteriyor. Sayfa 27’deki yedi Pareto grafiğinde, her бөлүм için hız homojenliği ile энергия диссипациясы 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 аралаштыргычnın bulunduğu dördüncü бөлүмde görüldü. Bölme sınırları аралаштыргычnın oluşturduğu yeniden dolaşım yapılarını daha iyi kapsayacak biçimde değiştirildiğinde, bölge içindeki энергия диссипациясы 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 жылмаma koşulu ve дубал fonksiyonları, sıvı-hava интерфейсünde ise жылма 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 дубал süresi gerektirdi.

Uzun süreli жылуулук өткөрүү nasıl hızlandırıldı?

Dört saatlik veya daha uzun bir кристаллдашуу 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 муздатуу ve ısıtma koşulları örneklenerek düşük mertebeli bir ısıl model oluşturuldu.

Her бөлүм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 бөлүм hacmini, F_j→i бөлүмдөр arasındaki hacimsel akışı, T_i бөлүм sıcaklığını, T_w дубал sıcaklığını ve h etkin ısı transfer katsayısını gösteriyor.

Yedi hidrodinamik бөлүм, ısıl hesaplamayı daha da azaltmak için iki ana grupta toplandı:

• Duvara komşu 1-3 numaralı бөлүмдөр,
• Merkez ve аралаштыргыч etkisindeki 4-7 numaralı бөлүмдөр.

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 өлчөмta izlendi

Model кристаллleri tek bir eşdeğer çapla değil, узундук l ve кеңдик w koordinatlarıyla temsil ediyor. Her бөлүмdeki sayı yoğunluğu n_z(l,w,t) ile gösteriliyor.

Basitleştirilmiş эки өлчөмдүү популяция балансы 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 кристаллдардын узундук ve кеңдик yönündeki өсүүsini, Dirac delta terimi yeni çekirdeklerin belirli bir başlangıç өлчөмunda oluşmasını, son terim ise кристаллдардын бөлүмдөр arasındaki taşınmasını ifade ediyor.

Başlangıç tohumu эки өлчөмдүү Gauss dağılımıyla диагностикаmlandı. Ortalama кристалл uzunluğu 200 µm, ortalama genişliği 100 µm; iki yöndeki dağılım parametreleri ise 150 µm olarak seçildi.

Model сынуу, aglomerasyon ve çözünmeyi ihmal ediyor. Bu nedenle sonuçlar өсүү ve ikincil ядролонууnin baskın olduğu кристаллдашуу koşullarıyla sınırlıdır.

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

Her бөлүмde bütün узундук-кеңдик 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 кристалл sayısını, ortalama өлчөмu ve toplam кристалл hacmini düşük sayıda diferansiyel denklemle temsil ediyor. Böylece temel кристаллдашуу dinamikleri standart моменттер ыкмасыyle 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 эки өлчөмдүү кристалл өлчөмu dağılımını yüksek эригичтикlü чектүү көлөм yöntemiyle yeniden kurdu.

Yüksek эригичтикlü чектүү көлөм yöntemi

HRFVM yönteminde узундук ve кеңдик alanı düzenli bir hesaplama ağına ayrıldı. Kristaller büyüdükçe sayı yoğunluğu küçük өлчөм hücrelerinden daha büyük өлчөм 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, өсүү ve ядролонуу

Her бөлүмnin aşırı doygunluğu yerel концентрация ile sıcaklığa bağlı эригичтикten 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 узундук ve кеңдик yönlerindeki өсүү 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, кристаллin iki yönde aynı hızla өсүүmesine ve тараптар катышыnın süreç boyunca değişmesine izin veriyor.

İkincil ядролонуу ise ашыкча каныгуу, mevcut кристалл hacmi ve yerel энергия диссипациясыyla ilişkilendirildi:

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

Enerji yayılımı etkisi etkinleştirildiğinde аралаштыргыч çevresindeki yüksek турбуленттүүлүкlı bölgelerde кристалл-кристалл çarpışmaları ve yüzey aşınmasıyla ilişkili ядролонуу kuvvetlenebiliyor.

Hesaplama өндүрүмдүүлүкı

İlk karşılaştırma, tek бөлүмli ve tamamen karışmış bir кристаллдашуу modeli üzerinde yapıldı. Tam HRFVM ile hibrit SMOM-HRFVM aynı ağ эригичтикlerinde 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 бөлүмli hesaplamada fark daha belirgin hale geldi:

Araştırmacılar genel hız kazanımını problem өлчөмuna bağlı olarak bir ile iki mertebe arasında диагностикаmlıyor. Tablodaki doğrudan karşılaştırılabilen en güçlü sonuç yaklaşık 23,5 kat ылдамдооdır.

Sıcaklık farkı концентрация eşitken bile өсүүyi değiştirdi

Uzaysal heterojenlik kademeli biçimde modele eklendi. İlk durumda бөлүмli yapı ideal karışım davranışına yakındı. İkinci durumda дубал kaynaklı sıcaklık farkları, sonraki durumlarda ise энергия диссипациясыna bağlı ядролонуу etkileri eklendi.

Sayfa 30’daki grafiklerde bütün бөлүмдөрin çözünen madde концентрацияi eğrileri neredeyse üst üste geliyor. Bunun nedeni, incelenen temel koşullarda бөлүмдөр arası karışımın кристаллдашуу kinetiğinden daha hızlı olmasıdır.

Buna karşılık aynı grafikteki ашыкча каныгуу ve кеңдик yönündeki өсүү hızı eğrileri iki gruba ayrılıyor. Duvara yakın 1-3 numaralı бөлүмдөр daha düşük sıcaklık nedeniyle daha yüksek aşırı doygunluğa ve daha yüksek өсүү hızına ulaşıyor. Merkezdeki 4-7 numaralı бөлүмдөр 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 концентрация ölçümü homojen görünse bile yerel sıcaklık farkları кристалл өсүүsini heterojen hale getirebilir.

Türbülans ядролонууnin yerini belirledi

Enerji yayılımının ядролонуу üzerinde etkili olmadığı ilk durumlarda бөлүмдөрin ядролонуу 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 ядролонуу grafikleri, аралаштыргычya yakın ve yüksek энергия диссипациясыna sahip бөлүмдөрin baskın ядролонуу merkezleri haline geldiğini gösteriyor. Daha düşük турбуленттүүлүкlı bölgelerde ise ядролонуу büyük ölçüde bastırıldı.

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

• Yüksek турбуленттүүлүкlı bölgelerde daha fazla yeni çekirdek oluşarak küçük кристалл nüfusu artıyor.
• Düşük ядролонууli bölgelerde mevcut кристаллler өсүүye devam ederek daha büyük parçacıklar oluşturuyor.

Kristal өлчөмu ve тараптар катышы nasıl değişti?

Nominal kinetik ve 0,0111 °C/s муздатуу hızında, sıcaklık heterojenliği tek başına nihai тараптар катышыnı çok az değiştirdi. Enerji yayılımına bağlı ядролонуу kuvvetlendikçe ortalama кристалл өлчөмları belirgin biçimde büyüdü.

Birinci ve beşinci durum karşılaştırıldığında hacim ağırlıklı кеңдик yaklaşık yüzde 24,5, узундук 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 эки өлчөмдүү dağılımlar, ortalama тараптар катышы benzer görünse bile parçacık yoğunluğunun узундук-кеңдик düzlemindeki konumunun ve dağılım biçiminin değiştiğini gösteriyor. Bu nedenle yalnızca ortalama тараптар катышыna bakmak, ürün dağılımındaki bütün farkları ortaya koymuyor.

Yavaş муздатуу daha uzun кристаллler oluşturdu

Soğutma hızı 0,0111 °C/s’den 0,0011 °C/s’ye düşürüldüğünde кристаллдардын өсүү için daha fazla zamanı oldu. Uzunluk yönündeki өсүү кеңдик 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ş муздатууda hacim ağırlıklı кеңдик 269,7 µm, узундук 785,1 µm ve тараптар катышы 2,91 olarak hesaplandı.

Bu bulgu, муздатуу programının yalnızca ortalama кристалл өлчөмunu değil, өсүүnin iki eksen arasındaki dağılımını da değiştirebildiğini gösteriyor.

Karışım ve кристаллдашуу süreleri yarıştığında ne oldu?

Araştırmacılar кристалл өсүү sabitlerini 90 kat artırarak кристаллдашуу 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 бөлүмдөр arası akış, çözünen madde концентрацияini hızla eşitliyordu. Hızlandırılmış kinetikte ise кристаллler aşırı doygunluğu karışımın eşitleyebileceğinden daha hızlı tüketmeye başladı. Bunun sonucunda tek bir концентрация eğrisi yerine yedi ayrı бөлүм eğrisi oluştu.

Duvara yakın daha soğuk bölgelerde өсүү ve çözünen madde tüketimi hızlandı. Merkez bölgeler ise бөлүмдөр arası alışverişle seyreldi. Böylece ölçek büyütmede belirleyici unsurun yalnızca аралаштыруу hızı veya yalnızca кристалл өсүү hızı olmadığı; bu iki sürenin birbirine oranı olduğu gösterildi.

Hızlandırılmış өсүүde ядролонуу hızları бөлүмдөр arasında büyük farklar gösterse de süreç өсүү baskın kaldı. Bu nedenle hacim ağırlıklı nihai кристалл dağılımları birbirine yakın kaldı. Sayfa 35’teki dağılımlarda тараптар катышы 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, кристаллдашуу ölçek büyütmede yalnızca tankın ortalama sıcaklığına, аралаштыргыч devrine veya ortalama aşırı doygunluğa bakmanın yeterli olmayabileceğini gösteriyor. Aynı ortalama концентрацияin arkasında farklı sıcaklık, өсүү ve ядролонуу geçmişleri bulunabilir.

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

• Laboratuvar ve üretim реакторleri arasındaki morfoloji değişiminin öngörülmesi,
• Karıştırıcı konumu ve hızının кристалл өлчөмuna etkisinin araştırılması,
• Duvar муздатуу programlarının ve termal çevrimlerin tasarlanması,
• İğnemsi veya levhamsı кристаллдардын тараптар катышыnın kontrol edilmesi,
• Çok sayıda işletme koşulunun dijital ortamda hızlı biçimde taranması,
• Hesaplama maliyeti düşük кристаллдашуу 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 эригичтик, өсүү, ядролонуу ve энергия диссипациясы 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ş көп бөлүмдүү эки өлчөмдүү популяция балансы 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 бөлүмli modellerde çözünürlüğe bağlı olarak yaklaşık 3,2-7,2 kat ылдамдоо sağladı.
• Yedi бөлүмli 100 × 100 ağda hesaplama süresi 1.245,6 saniyeden 53,1 saniyeye indi.
• 200 × 200 ve 300 × 300 көп бөлүмдүү 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 бөлүмдөр arası концентрация farkları ihmal edilebilir düzeyde kaldı.
• Duvar sıcaklığı farkları концентрация aynı kalırken ашыкча каныгуу ve өсүү hızını bölgesel olarak değiştirdi.
• Enerji yayılımına bağlı ядролонуу, аралаштыргыч çevresindeki bölgeleri baskın ядролонуу merkezlerine dönüştürdü.
• Nominal hızlı муздатууda тараптар катышы 2,23-2,25 aralığında kaldı.
• Yavaş муздатууda тараптар катышы yaklaşık 2,91-2,95 aralığına yükseldi.
• Büyüme sabitleri 90 kat artırıldığında концентрация alanı yedi ayrı bölgesel eğriye ayrıldı.
• Hızlandırılmış өсүүde güçlü bölgesel ядролонуу farklarına rağmen nihai hacim dağılı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ğ эригичтикlerindeki 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 реакторden alınmış hız, sıcaklık, кристалл өлчөмu dağılımı veya тараптар катышы ö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 кристалл өлчөмu ve тараптар катышы sonuçları fiziksel deneyle doğrulanmamıştır.
• CFD hız ve турбуленттүүлүк 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 бөлүмnin kendi içinde sıcaklık, концентрация ve энергия диссипациясы bakımından homojen olduğu kabul edilmiştir.
• Akış alanı кристаллдашуу 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 бөлүмyi iki sıcaklık bölgesine indirgemektedir.
• Çekirdeklenmenin энергия диссипациясыna bağımlılığı ampirik bir terimle temsil edilmiştir.
• Kinetik parametreler başka maddelere veya реакторlere 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 концентрация, ашыкча каныгуу ve өсүү grafikleriyle uyumlu görünmemektedir; bu durum preprintte bir форма 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, кристалл өлчөмu sonuçları, форма 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.