Bu çalışma, düşük sıcaklıkta jel oluşturan parafinli ham petrolün elastik davranıştan akışkan davranışa geçişini farklı kesme yüklemeleri altında tek bir reolojik modelle açıklamayı amaçlamıştır. Araştırmacılar klasik Maxwell modelini, mum kristali ağının mevcut durumunu 0 ile 1 arasında temsil eden λ yapısal parametresiyle genişletmiş; kayma modülünü ve etkin viskoziteyi bu parametreye bağlamıştır. Model, araştırma kapsamında yapılan salınımlı kesme deneylerinin yanı sıra sabit kesme hızı, kesme hızı histerezis döngüsü, sabit kesme gerilmesi ve basamaklı kesme gerilmesi deneylerinden alınan verilerle doğrulanmıştır. Önerilen model salınımlı kesmede ortalama %13,70; diğer dört yükleme biçiminde %2,10 ile %3,20 arasında ortalama mutlak sapmalar üretmiştir. Bununla birlikte çalışma hakem değerlendirmesinden geçmemiş bir preprinttir; model gerçek bir boru hattında, dinamik soğuma veya yeniden başlatma basıncı ölçümleriyle doğrulanmamıştır.

Modelin temel yeniliği, parafinli ham petrolü akma noktasından önce tamamen katı, akma noktasından sonra tamamen akışkan kabul etmemesidir. Bunun yerine malzemenin mum kristali ağı sağlamken elastik enerji depolayabildiği, kesme arttıkça ağın kademeli olarak parçalandığı ve kesme kaldırıldığında yeniden yapılanabildiği sürekli bir geçiş tanımlanmıştır. Böylece akma öncesi viskoelastik davranış, akma geçişi, kaymayla incelme ve zamana bağlı tiksotropi aynı matematiksel çerçevede ele alınmıştır.

Önerilen model karşılaştırılan bütün deneylerde mutlak olarak en iyi sonucu vermemiştir. Kesme hızı histerezis döngüsünde Teng Model III’ün ortalama mutlak sapması %1,87 iken yeni modelin sapması %3,02 olmuştur. Buna karşılık Teng Model III, kesme gerilmesinin kontrol edildiği problemlerde matematiksel çözüm güçlüğü taşıdığı için sabit gerilme ve gerilme basamağı karşılaştırmalarına alınmamıştır. Çalışmanın ana katkısı bu nedenle tek bir deneyde en düşük hatayı üretmekten çok, farklı yükleme biçimlerini çözülebilir ve sürekli bir modelde birleştirmesidir.

Araştırmanın temel sorusu nedir?

Araştırmanın temel sorusu, parafinli ham petrolün farklı kesme geçmişlerine verdiği mekanik tepkilerin tek bir bünye denklemiyle açıklanıp açıklanamayacağıdır. Bu tepki yalnızca viskozite değişiminden oluşmamaktadır. Malzeme düşük deformasyonda elastik bir katıya benzer davranabilmekte, daha yüksek deformasyonda hem elastik hem viskoz özellik gösterebilmekte, akma sonrasında ise viskoz akışın baskın olduğu bir duruma geçmektedir.

Araştırmacılar özellikle şu sorunları aynı model içinde çözmeye çalışmıştır:

• Akma öncesindeki geri döndürülebilir elastik ve viskoelastik deformasyonun tanımlanması.
• Akma noktasında gerilme veya modül eğrilerinde yapay süreksizlik oluşmaması.
• Mum kristali ağının kesmeyle parçalanması ve dinlenmeyle yeniden kurulmasının temsil edilmesi.
• Salınımlı ve dönel kesme deneylerinin aynı yapısal parametre yaklaşımıyla açıklanması.
• Kesme hızı kontrollü ve kesme gerilmesi kontrollü problemlerde matematiksel olarak çözülebilir bir denklem elde edilmesi.

Parafinli ham petrol neden karmaşık davranmaktadır?

Çalışmada parafinli ham petrolün soğumasıyla, sıvı fazda çözünmüş mum moleküllerinin kristalleşmeye başladığı açıklanmaktadır. Sıcaklık akma noktasına yaklaştığında bu kristaller birbirine bağlanarak üç boyutlu, sünger benzeri bir ağ meydana getirmektedir. Sıvı petrol bu ağın içinde tutulduğu için malzeme jel hâline gelmekte ve katı ile sıvı arasında karmaşık mekanik özellikler göstermektedir.

Oluşan mum kristali ağı dört temel makroskopik davranışla ilişkilendirilmiştir:

• Viskoelastisite: Malzeme hem elastik enerji depolar hem de akış sırasında enerji kaybeder.
• Akma gerilmesi: Kalıcı akışın başlaması için belirli bir gerilme düzeyinin aşılması gerekir.
• Tiksotropi: İç yapı kesme altında zamanla bozulur, kesme azaltıldığında veya kaldırıldığında yeniden yapılanabilir.
• Kaymayla incelme: Kesme hızı yükseldikçe görünen viskozite azalır.

İlk yükleme aşamasında kristal ağı büyük ölçüde sağlamdır ve sistemin tepkisi elastiktir. Kesme yoğunlaştıkça ağ önce daha büyük kümelere, ardından daha küçük yapılara ayrılmaktadır. Akma bölgesine ulaşıldığında elastik etkinin yerini giderek viskoz akış almaktadır. Bu geçiş ani bir fiziksel anahtar değişimi değil, ağın gelişen hasarıyla bağlantılı sürekli bir süreç olarak modellenmiştir.

Bu davranış boru hattı taşımacılığı açısından neden önemlidir?

Çalışmada modellemenin temel mühendislik bağlamı, soğuyan parafinli ham petrolün boru hattında jelleşmesi ve bir duruş sonrasında yeniden harekete geçirilmesidir. Boru hattı kapatıldığında kristal ağı zamanla yeniden kurulabilir. Akış tekrar başlatıldığında pompa basıncı bu yapının elastik deformasyonunu, kademeli parçalanmasını ve akışa dönüşmesini yenmek zorundadır.

Yalnızca sabit bir viskozite veya tek bir akma gerilmesi kullanan yaklaşım, yüklemenin geçmişini ve yapının zaman içinde değişmesini temsil edemeyebilir. Araştırmacılar bu nedenle farklı yükleme biçimlerini kapsayan bir modelin, yeniden başlatma basıncı ve güvenli taşıma hesapları için daha uygun bir kuramsal temel oluşturabileceğini savunmuştur.

Bununla birlikte çalışma doğrudan boru hattı basıncı hesaplamamış, pompa seçimi yapmamış ve saha yeniden başlatma deneyi gerçekleştirmemiştir. Mühendislik katkısı, bu tür hesaplarda gelecekte kullanılabilecek bir malzeme modeli geliştirilmesiyle sınırlıdır.

Önceki modellerde hangi eksiklikler belirlenmiştir?

Çalışma, tiksotropik malzeme modellerini mikroyapısal ve fenomenolojik olmak üzere iki ana grupta değerlendirmiştir. Mikroyapısal modeller parçacıklar arasındaki fiziksel bağlantıları doğrudan tanımlamaya çalışmaktadır. Ancak parafinli ham petrolün kimyasal bileşimi ve kristal biçimleri karmaşık olduğundan, kristaller arasındaki bağ parametrelerinin güvenilir biçimde belirlenmesi pratik değildir.

Fenomenolojik yapısal kinetik modeller ise mikroskobik ayrıntıların tamamını çözmek yerine, yapının genel durumunu λ gibi boyutsuz bir değişkenle temsil etmektedir. Çalışmada λ = 1 tamamen kurulmuş, λ = 0 ise tamamen parçalanmış yapıya karşılık gelmektedir.

Modelin temel fikri nedir?

Model, seri bağlanmış elastik yay ile viskoz sönümleyiciden oluşan klasik Maxwell yaklaşımını temel almaktadır. Klasik Maxwell modelinde kayma modülü G ve viskozite η sabittir. Yeni yaklaşımda ise her iki özellik mum kristali ağının mevcut durumuna göre değişmektedir.

Yapı sağlamken λ değeri 1’e yaklaşmakta, kayma modülü ve yapısal viskozite yüksek kalmaktadır. Kesme yapıyı bozdukça λ azalmakta, malzemenin elastik gerilme taşıma ve viskoz enerji dağıtma kapasitesi değişmektedir. Kesme azaldığında yapısal geri kazanım terimi λ’yı yeniden yükseltebilmektedir.

Klasik Maxwell denklemi neyi ifade etmektedir?

Çalışmada başlangıç noktası olarak kullanılan klasik Maxwell denklemi şöyledir:

\[ \frac{\eta}{G}\frac{d\tau}{dt}+\tau=\eta\frac{d\gamma}{dt} \]

Burada τ kayma gerilmesini ve birimi Pa’yı; γ kayma şekil değiştirmesini ve boyutsuz deformasyonu; dγ/dt veya γ̇ kayma hızını ve birimi s^-1’i; G kayma modülünü ve birimi Pa’yı; η ise dinamik viskoziteyi ve birimi Pa·s’yi göstermektedir.

η/G oranı bir zaman ölçeği oluşturur. Bu süre, malzemenin elastik gerilmeyi ne kadar hızlı gevşettiğini ifade etmektedir. Klasik biçimde G ve η sabit olduğu için kristal ağının kesme sırasında bozulması veya dinlenme sırasında yeniden kurulması tanımlanamamaktadır.

Kayma modülü yapısal parametreye nasıl bağlanmıştır?

Yapıya bağlı kayma modülü şu şekilde tanımlanmıştır:

\[ G_s(\lambda)=G_0\lambda^p \]

G_s(λ), mevcut yapısal durumdaki kayma modülüdür ve birimi Pa’dır. G₀, ağ tamamen kurulmuşken teorik kayma modülüdür. λ, 0 ile 1 arasında değişen boyutsuz yapısal parametredir. p ise modülün yapı bozuldukça ne kadar doğrusal olmayan biçimde azalacağını belirleyen boyutsuz üstür.

λ = 1 olduğunda G_s = G₀ olur. λ sıfıra yaklaştıkça kristaller arasındaki etkin bağlantılar azalmakta ve modelin elastik taşıma kapasitesi zayıflamaktadır. p büyüdükçe küçük bir yapısal kayıp modülde daha keskin bir azalma oluşturabilir.

Etkin viskozite neden değiştirilmiştir?

Klasik güç yasası viskozitesi şu genel biçimdedir:

\[ \mu=K|\dot{\gamma}|^{n-1} \]

Kaymayla incelen bir akışta n < 1 olduğunda, γ̇ sıfıra yaklaşırken bu ilişki sonsuz viskoziteye yaklaşabilir. Araştırmacılar, salınımlı deneylerin düşük kesme hızlı doğrusal bölgesinde bu matematiksel davranışın fiziksel olmayan ve uyumu bozan sonuçlar doğurduğunu belirtmiştir.

Bu sorunu önlemek için etkin viskozite şu şekilde tanımlanmıştır:

\[ \eta_{\mathrm{eff}}(\lambda,\dot{\gamma})=(\lambda\Delta k+k)\left(\dot{\gamma}_c+|\dot{\gamma}|\right)^{n-1} \]

η_eff, mevcut yapı ve kayma hızındaki etkin viskozitedir. k, yapı tamamen bozulduğunda kalan kıvam katsayısıdır ve birimi Pa·sⁿ’dir. Δk, yapının viskoz katkısının büyüklüğünü belirleyen katsayıdır; λΔk terimi yapı sağlamlaştıkça artmaktadır. n akış davranış indisi, γ̇_c ise doğrusal viskoelastik bölgeden doğrusal olmayan akış bölgesine geçişi temsil eden kritik kayma hızıdır ve birimi s^-1’dir.

γ̇ sıfıra yaklaşırken γ̇_c terimi sıfır olmayan bir taban oluşturduğu için viskozite sonsuza gitmemekte ve sonlu bir değere yaklaşmaktadır:

\[ \eta_{\mathrm{eff}}\rightarrow(\lambda\Delta k+k)\dot{\gamma}_c^{\,n-1} \]

Genişletilmiş Maxwell denklemi nasıl kurulmuştur?

Yapıya bağlı modül ve etkin viskozite klasik Maxwell denklemine yerleştirildiğinde çalışma şu kapalı biçimi vermektedir:

\[ \frac{\eta_{\mathrm{eff}}(\lambda,\dot{\gamma})}{G_s(\lambda)}\frac{d\tau}{dt}+\tau=\eta_{\mathrm{eff}}(\lambda,\dot{\gamma})\frac{d\gamma}{dt} \]

Bu denklemde gerilmenin zamana bağlı gelişimi, hem mevcut kristal ağına hem de uygulanan kayma hızına göre değişmektedir. Yapı sağlam olduğunda elastik depolama ve viskoz direnç yüksek; yapı parçalandığında ise parametreler daha düşük bir yapısal duruma karşılık gelmektedir.

Çalışmanın açılmış Denklem 5 ve Denklem 7 biçimlerinde sağ tarafta bağımsız ve tanımsız görünen bir “(λ, γ̇)” ifadesi bulunmaktadır. Bu nedenle burada denklem sessizce yeniden yazılmamış; çalışmada açık ve boyutsal olarak yorumlanabilir olan kapalı Denklem 4 kullanılmıştır.

Mum kristali ağının bozulması ve iyileşmesi nasıl modellenmiştir?

Yapısal parametrenin zamana göre değişimi şu hız denklemiyle verilmiştir:

\[ \frac{d\lambda}{dt}=k_{\mathrm{rec}}(1-\lambda)-k_{\mathrm{de}}\lambda|\dot{\gamma}|^c \]

İlk terim yapısal geri kazanımı, ikinci terim kesmenin neden olduğu yapısal parçalanmayı temsil etmektedir.

• k_rec: Yapısal geri kazanım hız sabitidir. Kesme bulunmadığında λ’yı 1’e doğru taşımaktadır.
• k_de: Kesmenin kristal ağını parçalama etkinliğini belirleyen hız sabitidir.
• c: Parçalanma hızının kayma hızına duyarlılığını belirleyen üstür.
• λ: Mevcut yapının tamamen bozulmuş ve tamamen kurulmuş durumlar arasındaki konumudur.
• |γ̇|: Kaymanın yönünden bağımsız büyüklüğüdür; salınımlı yüklemede pozitif ve negatif yarı çevrimlerin her ikisi de yapıyı bozabilir.

Kesme sıfırken parçalanma terimi ortadan kalkar ve λ, k_rec tarafından 1’e doğru yükselir. Kayma hızı arttığında parçalanma terimi güçlenir ve λ azalır. Böylece model yalnızca uygulanan anlık kaymaya değil, malzemenin önceki yükleme geçmişine de duyarlı hâle gelmektedir.

Depolama ve kayıp modülleri nasıl hesaplanmıştır?

Salınımlı kesmede depolama modülü G′, bir çevrim sırasında elastik olarak depolanan enerjiyi; kayıp modülü G′′ ise viskoz kayıpları temsil etmektedir. Parametreleştirilmiş Maxwell modeli için verilen ifadeler şöyledir:

\[ G'=G_s(\lambda)\frac{\left(\omega\frac{\eta_{\mathrm{eff}}(\lambda,\dot{\gamma})}{G_s(\lambda)}\right)^2}{1+\left(\omega\frac{\eta_{\mathrm{eff}}(\lambda,\dot{\gamma})}{G_s(\lambda)}\right)^2} \]

\[ G''=\frac{\omega\eta_{\mathrm{eff}}(\lambda,\dot{\gamma})}{1+\left(\omega\frac{\eta_{\mathrm{eff}}(\lambda,\dot{\gamma})}{G_s(\lambda)}\right)^2} \]

Burada ω = 2πf açısal frekanstır ve birimi rad/s’dir; f ise Hz cinsinden salınım frekansıdır. G′ ve G′′ değerlerinin birimi Pa’dır. λ hız denkleminden, η_eff etkin viskozite denkleminden ve G_s yapıya bağlı modül denkleminden elde edilmektedir.

Düşük deformasyonda ağ büyük ölçüde sağlam kaldığı için G′ yüksek ve yaklaşık sabittir. Deformasyon arttığında ağ bozulmakta, G′ düşmekte ve G′′ önce enerji kaybının artması nedeniyle yükselip daha sonra yapının ileri derecede parçalanmasıyla azalmaktadır.

Deneyde hangi ham petrol kullanılmıştır?

Araştırmacıların gerçekleştirdiği salınımlı kesme deneylerinde Daqing ham petrolü kullanılmıştır. Numunenin bildirilen temel özellikleri şöyledir:

34 °C’lik temel deney sıcaklığı, bildirilen 36 °C akma noktasının ve 49,85 °C mum görünüm sıcaklığının altındadır. Bu nedenle test sırasında mum kristali ağının gelişmesine izin veren düşük sıcaklıklı jel durumu incelenmiştir.

Numune hazırlama işlemi nasıl yapılmıştır?

Numunenin önceki sıcaklık ve kesme geçmişinden kaynaklanan “hafıza” etkisini azaltmak amacıyla petrol testten önce 80 °C’de işleme alınmıştır. Deney sırası şu şekilde verilmiştir:

1. Ön işleme tabi tutulan numune 50 °C’ye ısıtılmış ve 5 dakika bu sıcaklıkta tutulmuştur.
2. Numune reometreye yüklendikten sonra 5 dakika daha bekletilmiştir.
3. Sıcaklık dakikada 0,5 °C hızla test sıcaklığına düşürülmüştür.
4. Jel yapısının gelişmesi için numune test sıcaklığında 40 dakika hareketsiz tutulmuştur.
5. Bu işlemin ardından salınımlı kesme deneyi uygulanmıştır.

Bu hazırlama protokolü, her ölçümün benzer başlangıç yapısından başlamasını amaçlamaktadır. Çalışmada bağımsız deney tekrarlarının sayısı veya ölçümlerin tekrarlar arası değişkenliği verilmemiştir.

Salınımlı kesme deneyi nasıl gerçekleştirilmiştir?

Deneylerde Almanya üretimi HAAKE MARS III reometre ve 60 mm çapında parlatılmış koni-plaka rotor sistemi kullanılmıştır. HAAKE A25 programlanabilir su banyosunun sıcaklık kontrol doğruluğu ±0,01 °C olarak belirtilmiştir.

Ana deney sıcaklığı 34 °C, salınım frekansları ise 1, 2, 3, 4 ve 5 Hz’dir. Şekil değiştirme genliği artırılarak malzemenin doğrusal viskoelastik bölgeden doğrusal olmayan bölgeye ve akma sonrası bölgeye geçişi izlenmiştir.

İlk salınımlı kesme grafikleri ne göstermektedir?

Depolama modülü G′ grafiğinde yatay eksen yüzde cinsinden şekil değiştirme genliğini, dikey eksen Pa cinsinden G′ değerini göstermektedir. Her renk 1 ile 5 Hz arasındaki farklı bir frekansı temsil etmektedir. Düşük şekil değiştirmede G′ yaklaşık 10³-10⁴ Pa ölçeğinde bir plato göstermekte, akma bölgesine yaklaşıldığında kesintisiz biçimde düşmektedir.

Kayıp modülü G′′ grafiği doğrusal olmayan bir tepe davranışı göstermektedir. G′′ düşük deformasyonda görece sabit kalmakta, mum kristali ağı çözülmeye ve enerji kaybı artmaya başladığında yükselmekte, ileri deformasyonda yapı büyük ölçüde parçalandığında yeniden azalmaktadır.

Bu iki eğri birlikte değerlendirildiğinde malzemenin düşük deformasyonda elastik ağırlıklı, akma çevresinde güçlü enerji dağıtan viskoelastik ve yüksek deformasyonda viskoz akış ağırlıklı bir tepki verdiği görülmektedir. Grafikler ani bir mekanik sıçrama yerine sürekli bir geçiş göstermektedir.

Model nasıl eğitilmiş ve doğrulanmıştır?

Salınımlı kesme doğrulaması 32, 33, 34 ve 35 °C olmak üzere dört sıcaklıktaki toplam 14.000 veri noktası üzerinden yapılmıştır. Her sıcaklıkta deney verilerinin bir bölümü en küçük kareler yöntemiyle model parametrelerini belirlemek için kullanılmıştır. Elde edilen parametrelerle aynı sıcaklıktaki diğer frekansların modül değerleri hesaplanmış ve ölçümlerle karşılaştırılmıştır.

Çalışma, eğitim ve doğrulama için kullanılan veri noktalarının kesin oranını veya hangi frekansların her sıcaklıkta eğitim amacıyla seçildiğini toplu olarak açıklamamıştır. 34 °C örneğinde modelin 2 ve 4 Hz sonuçları ölçümlerle karşılaştırılmıştır.

Şekil 2’de noktalı eğriler ölçülen G′ ve G′′ değerlerini, düz kırmızı ve mavi eğriler model sonuçlarını göstermektedir. Model G′ eğrisindeki düzgün düşüşü ve G′′ eğrisindeki önce yükselip sonra azalan davranışı genel olarak yeniden üretmiştir.

Ortalama mutlak sapma sonuçları nasıl yorumlanmalıdır?

Çalışma model uyumunu ortalama mutlak sapma (average absolute deviation, AAD) yüzdeleriyle değerlendirmiştir. AAD’nin matematiksel tanımı metinde ayrıca verilmemiştir; bu nedenle burada yeni bir hata formülü eklenmemiştir. Daha küçük AAD, model eğrisinin ölçülen verilere genel olarak daha yakın olduğunu göstermektedir.

Dört sıcaklıkta ve toplam 14.000 salınımlı veri noktasında elde edilen genel ortalama AAD %13,70’tir. Sıcaklık ve frekanslara göre G′ ile G′′ sapmaları yaklaşık %10,83 ile %16,64 arasında değişmiştir.

Tablo 4’te verilen %11,98 değeri genel %13,70 sonucuyla çelişmemektedir. %11,98, yalnızca 34 °C’de 2 ve 4 Hz deneylerinden oluşan 3.500 veri noktalı model karşılaştırmasının ortalamasıdır. %13,70 ise dört sıcaklıktaki daha geniş salınımlı kesme doğrulamasının genel ortalamasıdır.

Salınımlı kesmede diğer modellerle karşılaştırma nasıldır?

Houska modeli bağımsız bir elastik eleman içermediği ve akma öncesini rijit kabul ettiği için G′ ve G′′ değerlerini tanımlayamamaktadır. Bu nedenle salınımlı kesme karşılaştırmasına alınmamıştır.

34 °C’de 2 ve 4 Hz için sonuçlar şöyledir:

Bu karşılaştırma, yeni modelin özellikle akma öncesi ve sonrası bölgeleri tek bir sürekli eğriyle birleştirme avantajını göstermektedir. Ancak salınımlı kesmedeki yaklaşık %12-%14 düzeyindeki hata, modelin bütün frekans ve deformasyon ayrıntılarını eksiksiz yeniden üretmediğini de göstermektedir.

Sabit kesme hızı deneyi ne göstermiştir?

Sabit kesme hızı doğrulamasında 34 °C’de 1, 2, 4 ve 8 s^-1 kesme hızlarına ait toplam 600 veri noktası kullanılmıştır. Bu veriler araştırmacıların bu çalışma kapsamında gerçekleştirdiği yeni deneylerden değil, belirtilen önceki çalışmadan alınmıştır.

Gerilme-zaman eğrileri başlangıçta yükselerek bir tepeye ulaşmış, ardından azalmış ve kararlı bir düzeye yaklaşmıştır. Başlangıçtaki yükseliş kristal ağının elastik gerilme taşımasına, tepe sonrasındaki düşüş ise yapısal parçalanma ve tiksotropik yumuşamaya karşılık gelmektedir.

Önerilen model bu yüklemede karşılaştırılan modeller arasında en düşük ortalama sapmayı üretmiştir.

Kesme hızı histerezis döngüsü ne anlatmaktadır?

Histerezis deneyinde kesme hızı önce artırılmış, ardından aynı rampayla azaltılmıştır. Yapı artan kesme sırasında parçalandığı için azalan kesme kolundaki gerilme aynı kayma hızındaki yükselen koldan farklı kalmış ve kapalı bir döngü oluşmuştur. Bu alan, yükleme geçmişinin ve tiksotropik yapısal değişimin etkisini göstermektedir.

Deneyler 34 ve 36 °C’de, 0,05 s^-1/s kesme hızı rampa oranında ve toplam 2.000 veri noktasıyla değerlendirilmiştir.

Bu sonuç, yeni modelin her veri kümesinde mutlaka en iyi olmadığını göstermektedir. Teng Model III histerezis eğrisini daha düşük sapmayla açıklamış, ancak aynı model gerilme kontrollü problemlerde çözüm güçlüğü nedeniyle sonraki karşılaştırmalara dahil edilememiştir.

Sabit kesme gerilmesi altında model nasıl davranmıştır?

Sabit gerilme deneyinde 34 °C’de 80, 90, 100 ve 110 Pa gerilmeler uygulanmış ve toplam 600 veri noktası kullanılmıştır. Ölçülen kesme hızı başlangıçta artmış, ardından kademeli biçimde kararlı düzeye yaklaşmıştır.

Önerilen model ile Teng Model IV bu eğri biçimini genel olarak izlemiştir. Houska modeli ise yaklaşık doğrusal yükseliş ve doğrusal kararlılaşma üreterek ölçümlerden daha fazla sapmıştır.

Teng Model III’ün bünye denklemi gerilme kontrollü yüklemede matematiksel çözüm güçlüğü oluşturduğu için bu karşılaştırmada sonuç verilmemiştir.

Basamaklı kesme gerilmesi deneyi ne göstermiştir?

Basamaklı deneyde gerilme belirli zamanlarda art arda daha yüksek düzeylere çıkarılmıştır. Modelin her yeni gerilme basamağında kesme hızındaki ani ve zamana bağlı değişimi izlemesi beklenmiştir.

34 °C’de 80, 90, 100 ve 110 Pa; 36 °C’de ise 30, 40, 50 ve 60 Pa basamakları kullanılmıştır. Toplam 1.920 veri noktası üzerinden elde edilen sonuçlar şöyledir:

Önerilen model en düşük ortalama hatayı üretmiş, ancak Teng Model IV ile arasındaki fark yalnızca 0,18 yüzde puan olmuştur. Çalışma bu fark için istatistiksel anlamlılık testi veya güven aralığı sunmamıştır.

Beş yükleme biçimindeki genel sonuç nedir?

Neden salınımlı kesmedeki hata daha yüksektir?

Çalışmaya göre salınımlı kesmede malzeme tek bir akış rejiminde kalmamaktadır. G′ ve G′′ değerleri doğrusal viskoelastik bölgeden doğrusal olmayan bölgeye, ardından akma sonrası davranışa geçerken birkaç aşamada değişmektedir. Ayrıca sonuçlar salınım frekansından etkilenmektedir.

Dönel kesme deneylerinde ise karşılaştırılan çıktı çoğunlukla tek bir gerilme veya kesme hızı eğrisidir. Bu nedenle salınımlı deneyde aynı modelin hem elastik enerji depolamasını hem viskoz kaybı hem frekans etkisini hem de yapısal bozulmayı eş zamanlı izlemesi gerekmektedir.

Modelin güçlü yönleri nelerdir?

• Akma öncesindeki viskoelastik davranışı model kapsamına almaktadır.
• Akma öncesi ve sonrası için ayrı parçalı denklemler kullanmadan sürekli geçiş sağlamaktadır.
• Kayma modülü ile etkin viskoziteyi aynı yapısal parametreye bağlamaktadır.
• Düşük kesme hızında güç yasasının sonsuz viskozite sorununu kritik kesme hızı terimiyle önlemektedir.
• Salınımlı, hız kontrollü ve gerilme kontrollü yüklemelerde aynı temel yapı kullanılabilmektedir.
• Gerilme kontrollü problemlerde cebirsel çözüm engeli bulunan Teng Model III’e göre daha geniş bir uygulama kapsamı sunmaktadır.
• Beş farklı yükleme biçimindeki eğilimleri genel olarak düşük hata oranlarıyla açıklamaktadır.

Çalışmanın sınırlılıkları nelerdir?

• Çalışma hakem değerlendirmesinden geçmemiş bir preprinttir.
• Yeni gerçekleştirilen deneysel bölüm salınımlı kesmeyle sınırlıdır; diğer dört yükleme biçimindeki veriler önceki çalışmalardan alınmıştır.
• Farklı veri kümeleri farklı ham petrol numunelerine, sıcaklıklara ve deney protokollerine ait olabilir; ortak model performansı bu heterojenlik içinde değerlendirilmiştir.
• Parametrelerin sıcaklıkla değişimini açıklayan tek bir evrensel bağıntı sunulmamıştır; her sıcaklık için parametre uyumu yapılmıştır.
• Model dinamik soğuma, boru boyunca sıcaklık gradyanı, basınç değişimi ve gerçek boru akışıyla doğrulanmamıştır.
• Bağımsız deney tekrar sayısı, ölçüm standart sapmaları ve parametre güven aralıkları verilmemiştir.
• Model parametrelerinin tekillik, korelasyon veya tanımlanabilirlik analizi sunulmamıştır.
• AAD karşılaştırmaları için istatistiksel anlamlılık testi yapılmamıştır.
• Akış içindeki kayma bantlaşması, duvar kayması veya uzamsal yapı değişimleri doğrudan ölçülmemiştir.
• Model mum kristali ağını λ adlı tek bir toplu değişkenle temsil etmektedir; kristal boyutu, biçimi ve bağ sayısı ayrı ayrı çözülmemektedir.

Çalışmanın söylediği sonuçlar

• Yapıya bağlı Maxwell modeli, parafinli ham petrolün akma öncesi ve sonrası davranışını sürekli biçimde temsil edebilmektedir.
• Yapısal parametreye bağlanan modül ve viskozite, kristal ağının bozulma ve yeniden kurulma etkilerini tek çerçevede birleştirmektedir.
• Model salınımlı kesmede karşılaştırılan Teng modellerinden daha düşük hata üretmiştir.
• Sabit kesme hızı, sabit gerilme ve gerilme basamağı deneylerinde karşılaştırılan modeller içinde en düşük ortalama AAD değerlerini vermiştir.
• Kesme hızı histerezisinde Teng Model III daha düşük hata üretmiştir.
• Model farklı yükleme biçimlerinde uygulanabilir bir kuramsal temel sunmaktadır.

Çalışmanın kanıtlamadığı sonuçlar

• Modelin gerçek bir petrol boru hattındaki yeniden başlatma basıncını doğru tahmin ettiği kanıtlanmamıştır.
• Pompa gücü, enerji tüketimi veya işletme maliyetinde azalma ölçülmemiştir.
• Model Türkiye’deki belirli bir ham petrol veya boru hattı üzerinde doğrulanmamıştır.
• Tek bir parametre takımının bütün sıcaklıklara, ham petrol türlerine ve kesme geçmişlerine uygulanabileceği gösterilmemiştir.
• Önerilen model her yükleme biçiminde diğer bütün modellerden daha iyi değildir.
• Düşük AAD değerleri, modelin mikroskobik mum kristali yapısını fiziksel olarak eksiksiz temsil ettiğini kanıtlamamaktadır.
• Endüstriyel ölçekli yazılım entegrasyonu veya gerçek zamanlı boru hattı kontrolü uygulanmamıştır.

Türkiye açısından ne ifade etmektedir?

Çalışmanın Türkiye açısından aktarılabilir yönü, soğuk koşullarda jelleşebilen parafinli ham petrolün boru hattındaki davranışını tek bir sabit viskozite veya keskin akma sınırıyla açıklamanın yetersiz kalabileceğini göstermesidir. Model, kesme geçmişini, elastik deformasyonu, yapısal bozulmayı ve yeniden yapılanmayı birlikte dikkate alan bir malzeme tanımı sunmaktadır.

Bu yaklaşım, gelecekte boru hattı duruş ve yeniden başlatma simülasyonlarında kullanılabilecek bir bünye denklemi adayıdır. Ancak Türkiye’de uygulama için yerel ham petrol numuneleriyle deney, boru çapı ve sıcaklık profilleri, duruş süresi, gerçek basınç kayıpları ve saha yeniden başlatma verileriyle ayrı doğrulama gerekmektedir.

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

Araştırma tasarımının teknik özeti

Deney ve veri kaynaklarının ayrımı

Karşılaştırmalı performans tablosu

Salınımlı kesmedeki %11,98 değeri 34 °C ve iki frekanstaki karşılaştırmalı tablodan gelmektedir. Dört sıcaklığı kapsayan geniş doğrulamanın ortalama değeri %13,70’tir.

Grafiklerin teknik mesajı

• Şekil 1: G′ düşük genlikte plato göstermekte ve akma sonrasında sürekli düşmektedir. G′′ önce yükselip daha sonra azalmaktadır.
• Şekil 2: Önerilen model 34 °C’de 2 ve 4 Hz için ölçülen G′ ve G′′ eğilimlerini genel olarak izlemektedir.
• Şekil 3: Teng Model IV akma noktasında süreksizlik, Teng Model III ise modül büyüklüklerinde zayıf uyum göstermektedir.
• Şekil 4: Sabit kesme hızında gerilme önce tepe yapmakta, yapı bozuldukça düşmekte ve kararlı değere yaklaşmaktadır.
• Şekil 5: Artan ve azalan kesme hızı kolları arasında tiksotropik histerezis döngüsü oluşmaktadır.
• Şekil 6: Sabit gerilmede kesme hızı zamana bağlı artıp kararlılaşmakta; önerilen model bu eğriyi Houska modelinden daha iyi izlemektedir.
• Şekil 7: Gerilme basamaklarında kesme hızı her basamağa karşılık yükselmekte; önerilen model ile Teng IV ölçümlere yakın sonuç vermektedir.

İstatistiksel ve yöntemsel değerlendirme

Model parametreleri en küçük kareler yöntemiyle belirlenmiş ve başarım AAD yüzdeleriyle raporlanmıştır. Çalışmada p değerleri, güven aralıkları, çapraz doğrulama dağılımları veya bağımsız test numuneleri için performans aralıkları verilmemiştir. Veri noktası sayısının yüksek olması, bağımsız numune veya deney tekrar sayısının aynı ölçüde yüksek olduğu anlamına gelmemektedir.

Model karşılaştırmalarında aynı veri üzerindeki AAD değerleri kullanılmıştır; ancak modellerin parametre sayısı, bilgi ölçütleri, hesaplama süresi veya aşırı uyum riski ortak bir ölçütle karşılaştırılmamıştır. Bu nedenle sonuçlar öncelikle eğri uyumu ve matematiksel uygulanabilirlik açısından değerlendirilmelidir.

Metinde belirlenen iç tutarsızlıklar

• Özette “structural dynamics theory”, ana metinde ise “structural kinetics theory” ifadesi kullanılmıştır.
• Denklem 5 ve 7’nin açılmış sağ tarafında tanımlanmayan bağımsız bir “(λ, γ̇)” çarpanı görünmektedir.
• Houska modeline verilen [17] ve [18] numaralı metin içi atıflar, kaynakçada aynı numaralarda yer alan eserlerle açıkça eşleşmemektedir.
• SSRN metadata kaydı Lei Wang’ı yazar listesinde göstermemekte, yüklenen çalışma ise Lei Wang’ı üçüncü yazar olarak vermektedir.

Kaynak ve Yöntem Notu

Çalışmanın tam özgün adı: Study on the rheological model of waxy crude oil under multi-mode shear loading

Yüklenen çalışmadaki yazarlar: Liping Guo; Jiajun Zhang; Lei Wang; Chenzi Shi; Baojun Liu.

Yazar sıralaması: Yukarıdaki liste yüklenen çalışmanın ilk sayfasındaki özgün sırayı korumaktadır.

Metadata uyuşmazlığı: SSRN’nin 6944013 numaralı kayıt sayfası Lei Wang’ı yazar listesinde göstermemektedir. Yüklenen çalışmada Lei Wang üçüncü yazardır. Bu nedenle yayımdan önce SSRN metadata kaydının veya yazar listesinin araştırmacılarla doğrulanması uygun olacaktır.

Eş birinci yazar: Çalışmada eş birinci yazarlık veya eş katkı beyanı bulunmamaktadır.

Sorumlu yazar: Liping Guo.

Kurum: Key Laboratory for Enhanced Oil & Gas Recovery of the Ministry of Education, Northeast Petroleum University, Daqing 163318, China.

DOI:10.2139/ssrn.6944013

Dergi: Hakemli bir dergi veya nihai dergi yayını bu sürüm üzerinden doğrulanmamıştır.

Özgün yayınevi: Nihai bir dergi yayınevi bulunmamaktadır. SSRN bir yayınevi veya hakemli dergi değil, Elsevier tarafından işletilen preprint ve erken araştırma platformudur.

Yayın platformu: SSRN.

Yayın yılı: 2026.

SSRN sürüm tarihi: 6944013 numaralı mevcut sürümün kayıt tarihi 15 Haziran 2026 olarak görünmektedir.

Kaynak türü: Yeni bir reolojik bünye modeli geliştiren, salınımlı kesme deneyi ve önceki deney verileriyle model doğrulaması içeren preprint araştırma çalışması.

Hakemlik durumu: Bu çalışma hakem değerlendirmesinden geçmemiş bir preprinttir; sonuçları bu sınırlılık dikkate alınarak okunmalıdır.

Resmî kaynak bağlantısı:https://ssrn.com/abstract=6944013

Finansman: Çalışma, National Natural Science Foundation of China tarafından 52574082 numaralı destek kapsamında finanse edilmiştir.

Çıkar çatışması: Yazarlar, çalışmayı etkileyebilecek bilinen mali çıkar veya kişisel ilişki bulunmadığını beyan etmiştir.

Yazar katkıları: Ayrıntılı CRediT yazar katkı beyanı çalışmada yer almamaktadır.

Bu Türkçe bilimsel açıklama, yüklenen çalışmanın metni, denklemleri, tabloları ve grafiklerinin baştan sona incelenmesiyle hazırlanmıştır. Bilimsel içerik yalnızca çalışmaya dayanmaktadır; çalışmada bulunmayan bilimsel bulgu veya dış literatür sonucu eklenmemiştir. Dış kaynaklar yalnızca DOI, SSRN kaydı, sürüm tarihi, platform niteliği ve bibliyografik kimlik kontrolü amacıyla kullanılmıştır.

Çalışmanın temel sınırlılıkları; hakem değerlendirmesinden geçmemesi, yeni deneylerin salınımlı kesmeyle sınırlı olması, diğer yükleme verilerinin önceki çalışmalardan alınması, bağımsız tekrar ve parametre belirsizliği bilgilerinin verilmemesi, sıcaklığın modele evrensel bir değişken olarak dahil edilmemesi ve gerçek boru hattı koşullarında doğrulama yapılmamasıdır. Önerilen model, boru hattı güvenliği veya yeniden başlatma basıncı için doğrulanmış saha garantisi olarak değerlendirilmemelidir.