2025-11-22T10:49:15.885853

Transient rod-climbing in an Oldroyd-B fluid

Ruangkriengsin, Brandão, Wu et al.
The Weissenberg effect, or rod-climbing phenomenon, occurs in non-Newtonian fluids where the fluid interface ascends along a rotating rod. Despite its prominence, theoretical insights into this phenomenon remain limited. In earlier work, Joseph \& Fosdick (\emph{Arch. Rat. Mech. Anal.}, vol. 49, 1973, pp. 321--380) employed domain perturbation methods for second-order fluids to determine the equilibrium interface height by expanding solutions based on the rotation speed. In this work, we investigate the time-dependent interface height through asymptotic analysis with dimensionless variables and equations using the Oldroyd-B model. We begin by neglecting surface tension and inertia to focus on the interaction between gravity and viscoelasticity. In the small-deformation scenario, the governing equations indicate the presence of a boundary layer in time, where the interface rises rapidly over a short time scale before gradually approaching a steady state. By employing a stretched time variable, we derive the transient velocity field and corresponding interface profile on this short time scale and recover the steady-state profile on a longer time scale. Subsequently, we reintroduce small but finite inertial effects to investigate their interplay with viscoelasticity and propose a criterion for determining the conditions under which rod-climbing occurs.
academic

Oldroyd-B তরলে ক্ষণস্থায়ী রড-আরোহণ

মৌলিক তথ্য

  • পেপার আইডি: 2501.01016
  • শিরোনাম: Oldroyd-B তরলে ক্ষণস্থায়ী রড-আরোহণ
  • লেখক: Tachin Ruangkriengsin, Rodolfo Brandão, Katie Wu, Jonghyun Hwang, Evgeniy Boyko, Howard A. Stone
  • শ্রেণীবিভাগ: physics.flu-dyn (তরল গতিবিদ্যা)
  • প্রকাশনার সময়: ২০২৫ সালের ২ জানুয়ারি (arXiv-এ জমা দেওয়া)
  • জার্নাল: Journal of Fluid Mechanics (জমা দেওয়া হয়েছে)
  • পেপার লিঙ্ক: https://arxiv.org/abs/2501.01016

সারসংক্ষেপ

এই পেপারটি অ-নিউটোনীয় তরলে Weissenberg প্রভাব (রড-আরোহণ ঘটনা) অধ্যয়ন করে, যেখানে তরলের ইন্টারফেস ঘূর্ণায়মান রডের উপর উঠে যায়। লেখকরা Oldroyd-B মডেল ব্যবহার করে, অ্যাসিম্পটোটিক বিশ্লেষণ এবং অ-মাত্রিক চলক মাধ্যমে সময়-নির্ভর ইন্টারফেস উচ্চতা অধ্যয়ন করেছেন। গবেষণায় দেখা গেছে যে ছোট বিকৃতির ক্ষেত্রে, নিয়ন্ত্রক সমীকরণ একটি সময় সীমানা স্তরের অস্তিত্ব নির্দেশ করে, যেখানে ইন্টারফেস স্বল্প সময়ের স্কেলে দ্রুত উঠে এবং তারপর ধীরে ধীরে স্থিতিশীল অবস্থার দিকে এগিয়ে যায়। প্রসারিত সময় চলক প্রবর্তন করে, লেখকরা ক্ষণস্থায়ী বেগ ক্ষেত্র এবং সংশ্লিষ্ট ইন্টারফেস প্রোফাইল উদ্ভাবন করেছেন এবং ছোট কিন্তু সীমিত জড়তা প্রভাব পুনরায় প্রবর্তন করে তাদের সাথে viscoelasticity-এর মিথস্ক্রিয়া অধ্যয়ন করেছেন।

গবেষণা পটভূমি এবং প্রেরণা

সমস্যা সংজ্ঞা

রড-আরোহণ ঘটনা (Weissenberg প্রভাব) অ-নিউটোনীয় তরলের সবচেয়ে বৈশিষ্ট্যপূর্ণ আচরণগুলির মধ্যে একটি। যখন একটি রড viscoelastic তরলে ঘোরে, তরলের পৃষ্ঠ রডের সাথে উপরের দিকে আরোহণ করে। এই ঘটনাটি তাত্ত্বিকভাবে এখনও গভীর বোঝাপড়ার অভাব রয়েছে।

গুরুত্ব

  1. তাত্ত্বিক তাৎপর্য: এই ঘটনাটি জটিল তরলের আচরণ বোঝার একটি মৌলিক ধারণা এবং শিক্ষা ও গবেষণায় গুরুত্বপূর্ণ ভূমিকা পালন করে
  2. ব্যবহারিক মূল্য: রিওলজি পরিমাপে গুরুত্বপূর্ণ প্রয়োগ রয়েছে, যা স্বাভাবিক চাপ সহগ এবং শিথিলকরণ সময়ের মতো উপাদান পরামিতি নির্ধারণ করতে ব্যবহৃত হয়
  3. প্রকৌশল প্রয়োগ: পলিমার দ্রবণ এবং অন্যান্য জটিল তরলের শিল্প প্রয়োগের জন্য নির্দেশনা প্রদান করে

বিদ্যমান পদ্ধতির সীমাবদ্ধতা

  1. তাত্ত্বিক সীমাবদ্ধতা: পূর্ববর্তী গবেষণা প্রধানত স্থিতিশীল অবস্থার আচরণে কেন্দ্রীভূত, দ্বিতীয় ক্রমের তরল মডেল ব্যবহার করে
  2. মডেল সমস্যা: দ্বিতীয় ক্রমের তরল মডেল শুধুমাত্র ছোট Deborah সংখ্যা বা Weissenberg সংখ্যায় কার্যকর এবং উচ্চ অস্থিরতা প্রদর্শন করে
  3. সময় নির্ভরতা: ক্ষণস্থায়ী গতিশীলতার তাত্ত্বিক বিশ্লেষণের অভাব রয়েছে

গবেষণা প্রেরণা

লেখকরা দ্বিতীয় ক্রমের তরল মডেলের অ-স্থিতিশীল প্রবাহে সমস্যা এড়াতে Oldroyd-B মডেল নির্বাচন করেছেন এবং প্রথমবারের মতো রড-আরোহণ ঘটনার সময় নির্ভরতা সিস্টেমেটিকভাবে অধ্যয়ন করেছেন।

মূল অবদান

  1. প্রথম তাত্ত্বিক বিশ্লেষণ: Oldroyd-B তরলে রড-আরোহণ ঘটনার প্রথম সম্পূর্ণ ক্ষণস্থায়ী তাত্ত্বিক বিশ্লেষণ প্রদান করে
  2. সময় সীমানা স্তর তত্ত্ব: সময় সীমানা স্তরের কাঠামো আবিষ্কার এবং বিশ্লেষণ করে, দ্রুত উত্থান এবং স্থিতিশীল অবস্থার দিকে ধীর অগ্রগতির দুই-পর্যায়ের প্রক্রিয়া প্রকাশ করে
  3. অ-মাত্রিক বিশ্লেষণ কাঠামো: Weissenberg সংখ্যার উপর ভিত্তি করে একটি সিস্টেমেটিক অ-মাত্রিক বিশ্লেষণ পদ্ধতি প্রতিষ্ঠা করে
  4. জড়তা প্রভাব গবেষণা: প্রথমবারের মতো ছোট কিন্তু সীমিত জড়তা প্রভাব এবং viscoelasticity-এর মিথস্ক্রিয়া অধ্যয়ন করে
  5. আরোহণ মানদণ্ড: রড-আরোহণ ঘটনার শর্ত নির্ধারণের জন্য একটি তাত্ত্বিক মানদণ্ড প্রস্তাব করে: Re<4βpWiRe < 4\beta_p Wi

পদ্ধতির বিস্তারিত বর্ণনা

কাজের সংজ্ঞা

অসীম দীর্ঘ বৃত্তাকার রড (ব্যাসার্ধ aa) যা কৌণিক বেগ Ω\Omega সহ viscoelastic পলিমার দ্রবণে ঘোরে, তরলের মুক্ত পৃষ্ঠ উচ্চতা h(r,t)h(r,t) এর সময় বিবর্তন নিয়ম অধ্যয়ন করা।

গাণিতিক মডেল

নিয়ন্ত্রক সমীকরণ

Oldroyd-B গঠনমূলক মডেল ব্যবহার করা হয়েছে:

  • কনফিগারেশন টেনসর বিবর্তন সমীকরণ: At+uAAu(u)TA=1λ(AI)\frac{\partial\mathbf{A}}{\partial t} + \mathbf{u} \cdot \nabla\mathbf{A} - \mathbf{A} \cdot \nabla\mathbf{u} - (\nabla\mathbf{u})^T \cdot \mathbf{A} = -\frac{1}{\lambda}(\mathbf{A} - \mathbf{I})
  • চাপ টেনসর: σ=pI+2μsE+G(AI)\boldsymbol{\sigma} = -p\mathbf{I} + 2\mu_s\mathbf{E} + G(\mathbf{A} - \mathbf{I})
  • গতিবেগ সমীকরণ: σ=ρgez\nabla \cdot \boldsymbol{\sigma} = \rho g\mathbf{e}_z

অ-মাত্রিকীকরণ

মূল অ-মাত্রিক পরামিতি প্রবর্তন করা হয়েছে:

  • Weissenberg সংখ্যা: Wi=λΩWi = \lambda\Omega
  • মাধ্যাকর্ষণ পরামিতি: G=ρga(μs+μp)ΩG = \frac{\rho ga}{(\mu_s + \mu_p)\Omega}
  • Reynolds সংখ্যা: Re=ρa2Ωμs+μpRe = \frac{\rho a^2\Omega}{\mu_s + \mu_p}

বৈশিষ্ট্যপূর্ণ সীমাবদ্ধতা

লেখকরা একটি মূল বৈশিষ্ট্যপূর্ণ সীমাবদ্ধতা চিহ্নিত করেছেন: G=WiGG = Wi \cdot \mathcal{G}, যেখানে G=O(1)\mathcal{G} = O(1), এটি দ্বি-পরামিতি অ্যাসিম্পটোটিক সম্প্রসারণকে একক-পরামিতি সম্প্রসারণে সরল করে।

প্রযুক্তিগত উদ্ভাবন পয়েন্ট

1. সময় স্কেল বিচ্ছেদ

দুটি সময় স্কেলের অস্তিত্ব আবিষ্কৃত হয়েছে:

  • বাহ্যিক স্তরের সময়: T=O(1)T = O(1), স্থিতিশীল অবস্থার আচরণের সাথে সংশ্লিষ্ট
  • অভ্যন্তরীণ স্তরের সময়: τ=T/Wi\tau = T/Wi, ক্ষণস্থায়ী সীমানা স্তরের সাথে সংশ্লিষ্ট

2. ডোমেইন বিঘ্ন পদ্ধতি

ছোট ইন্টারফেস বিকৃতির জন্য, সীমানা শর্তকে Z=H(R,T)Z = H(R,T) থেকে Z=0Z = 0 এ Taylor সম্প্রসারিত করা হয়েছে, মুক্ত পৃষ্ঠ সীমানা শর্তের পরিচালনা সরল করে।

3. Weber রূপান্তর

জড়তা প্রভাব পরিচালনা করার সময়, অর্ধ-অসীম ডোমেইনে Bessel ডিফারেনশিয়াল অপারেটর পরিচালনা করতে Weber রূপান্তর ব্যবহার করা হয়েছে, আংশিক ডিফারেনশিয়াল সমীকরণকে সাধারণ ডিফারেনশিয়াল সমীকরণে রূপান্তরিত করে।

পরীক্ষামূলক সেটআপ

তাত্ত্বিক বিশ্লেষণ কাঠামো

এই পেপারটি একটি বিশুদ্ধ তাত্ত্বিক গবেষণা, অ্যাসিম্পটোটিক বিশ্লেষণ পদ্ধতি ব্যবহার করে:

  1. ছোট পরামিতি: Wi1Wi \ll 1 (দুর্বল viscoelasticity)
  2. মাধ্যাকর্ষণ পরিসীমা: GWi1G \gg Wi^{-1} (শক্তিশালী মাধ্যাকর্ষণ প্রভাব)
  3. ছোট বিকৃতি অনুমান: ইন্টারফেস বিকৃতি তুলনামূলকভাবে ছোট

বিশ্লেষণ কৌশল

  1. চতুর্থ বিভাগ: বাহ্যিক স্তর বিশ্লেষণ - স্থিতিশীল ইন্টারফেস প্রোফাইল
  2. পঞ্চম বিভাগ: অভ্যন্তরীণ স্তর বিশ্লেষণ - ক্ষণস্থায়ী ইন্টারফেস প্রোফাইল
  3. ষষ্ঠ বিভাগ: ছোট জড়তা প্রভাব অন্তর্ভুক্ত করে বিশ্লেষণ

পরীক্ষামূলক ফলাফল

প্রধান ফলাফল

1. স্থিতিশীল ইন্টারফেস প্রোফাইল

Hsteady(R)=2βpWi2GR4H_{steady}(R) = \frac{2\beta_p Wi^2}{\mathcal{G}R^4}

এটি Joseph & Fosdick (1973) দ্বিতীয় ক্রমের তরল মডেল ব্যবহার করে প্রাপ্ত ফলাফলের সাথে সামঞ্জস্যপূর্ণ।

2. ক্ষণস্থায়ী ইন্টারফেস প্রোফাইল

Htransient(R,τ)=2βp(1(τ+1)eτ)Wi2GR4H_{transient}(R,\tau) = \frac{2\beta_p(1-(\tau+1)e^{-\tau})Wi^2}{\mathcal{G}R^4}

যেখানে τ=T/Wi\tau = T/Wi প্রসারিত সময় চলক।

3. আরোহণ মানদণ্ড

জড়তা প্রভাব অন্তর্ভুক্ত করার সময়, রড-আরোহণ ঘটনার শর্ত হল: 4Elβp>14El\beta_p > 1 যেখানে El=Wi/ReEl = Wi/Re হল স্থিতিস্থাপকতা সংখ্যা।

সময় বিবর্তন বৈশিষ্ট্য

  1. দ্রুত উত্থান পর্যায়: ইন্টারফেস τ6\tau \approx 6 সময়ের মধ্যে দ্রুত উঠে যায়
  2. স্থিতিশীল অবস্থার কাছাকাছি: ক্ষণস্থায়ী সমাধান দীর্ঘ সময় সীমায় স্থিতিশীল সমাধানে রূপান্তরিত হয়
  3. অ-একঘেয়ে আচরণ: জড়তা অন্তর্ভুক্ত করার সময়, ইন্টারফেস উচ্চতা প্রথমে উঠতে পারে তারপর নেমে আসতে পারে

মূল আবিষ্কার

  1. পরিধীয় চাপ: পলিমার প্রসারণ দ্বারা উৎপাদিত পরিধীয় চাপ আরোহণ চালনার প্রধান প্রক্রিয়া
  2. জড়তা প্রতিযোগিতা: জড়তা দ্বারা উৎপাদিত কেন্দ্রাতিগ শক্তি viscoelasticity-এর পরিধীয় চাপের সাথে প্রতিযোগিতা করে
  3. সীমানা স্তর কাঠামো: সময় সীমানা স্তরের অস্তিত্ব দ্রুত ক্ষণস্থায়ী এবং ধীর স্থিতিশীল অবস্থার দুই-পর্যায়ের আচরণ ব্যাখ্যা করে

সম্পর্কিত কাজ

ঐতিহাসিক উন্নয়ন

  1. Serrin (1959): প্রথম তাত্ত্বিক প্রচেষ্টা, দ্বিতীয় ক্রমের তরল মডেল ব্যবহার করে
  2. Giesekus (1961): জড়তা প্রভাব উপেক্ষা করে আরও ব্যাপক গঠনমূলক সমীকরণ
  3. Joseph & Fosdick (1973): সিস্টেমেটিক ডোমেইন বিঘ্ন পদ্ধতি কাঠামো
  4. Oldroyd (1950, 1951): প্রাথমিক Oldroyd-B মডেল প্রয়োগ

এই পেপারের সুবিধা

  1. সময় নির্ভরতা: প্রথম সম্পূর্ণ ক্ষণস্থায়ী বিশ্লেষণ
  2. মডেল স্থিতিশীলতা: দ্বিতীয় ক্রমের তরল মডেলের অস্থিরতা সমস্যা এড়ায়
  3. অ-মাত্রিক কাঠামো: সিস্টেমেটিক অ-মাত্রিক পরামিতি বিশ্লেষণ
  4. জড়তা প্রভাব: প্রথমবারের মতো জড়তা এবং viscoelasticity-এর মিথস্ক্রিয়া অধ্যয়ন

সিদ্ধান্ত এবং আলোচনা

প্রধান সিদ্ধান্ত

  1. সীমানা স্তর কাঠামো: রড-আরোহণ ঘটনায় স্পষ্ট সময় সীমানা স্তর রয়েছে, ইন্টারফেস দ্রুত উঠে তারপর ধীরে ধীরে স্থিতিশীল অবস্থার দিকে এগিয়ে যায়
  2. স্কেল বিচ্ছেদ: ক্ষণস্থায়ী সময় স্কেল O(Wi)O(Wi), স্থিতিশীল সময় স্কেল O(1)O(1)
  3. আরোহণ প্রক্রিয়া: পলিমার প্রসারণ দ্বারা উৎপাদিত পরিধীয় চাপ আরোহণ চালনার মূল প্রক্রিয়া
  4. জড়তা প্রভাব: ছোট জড়তা প্রভাব আরোহণ আচরণকে উল্লেখযোগ্যভাবে প্রভাবিত করতে পারে, একটি সমালোচনামূলক স্থিতিস্থাপকতা সংখ্যা রয়েছে

সীমাবদ্ধতা

  1. ছোট বিকৃতি অনুমান: শুধুমাত্র ছোট ইন্টারফেস বিকৃতির ক্ষেত্রে প্রযোজ্য
  2. দুর্বল viscoelasticity: Wi1Wi \ll 1 পরিসীমায় সীমাবদ্ধ
  3. শক্তিশালী মাধ্যাকর্ষণ: GWi1G \gg Wi^{-1} শর্ত প্রয়োজন
  4. সরলীকৃত মডেল: পৃষ্ঠ টেনশন এবং যোগাযোগ কোণ প্রভাব উপেক্ষা করে

ভবিষ্যত দিকনির্দেশনা

  1. উচ্চ Weissenberg সংখ্যা: Wi1Wi \gg 1 ক্ষেত্রে আচরণ অধ্যয়ন করা
  2. জটিল গঠনমূলক মডেল: FENE-P, Giesekus ইত্যাদি মডেলে সম্প্রসারণ
  3. সীমিত বিকৃতি: বড় ইন্টারফেস বিকৃতি ক্ষেত্রে পরিচালনা করা
  4. পরীক্ষামূলক যাচাইকরণ: পরীক্ষামূলক ডেটার সাথে তুলনা যাচাইকরণ

গভীর মূল্যায়ন

সুবিধা

  1. তাত্ত্বিক কঠোরতা: গাণিতিক অনুমান কঠোর, অ্যাসিম্পটোটিক বিশ্লেষণ পদ্ধতি যথাযথভাবে প্রয়োগ করা হয়েছে
  2. ভৌত অন্তর্দৃষ্টি: সময় সীমানা স্তর কাঠামো এবং আরোহণ প্রক্রিয়া প্রকাশ করে
  3. পদ্ধতি উদ্ভাবন: Weber রূপান্তর এবং ডোমেইন বিঘ্ন পদ্ধতির চতুর সমন্বয়
  4. ব্যবহারিক মূল্য: রিওলজি পরিমাপের জন্য তাত্ত্বিক ভিত্তি প্রদান করে

অপূর্ণতা

  1. পরামিতি পরিসীমা: প্রযোজ্য পরামিতি পরিসীমা তুলনামূলকভাবে সীমিত
  2. পরীক্ষামূলক তুলনা: পরীক্ষামূলক ডেটার সাথে সরাসরি তুলনার অভাব
  3. সংখ্যাগত যাচাইকরণ: কিছু ফলাফল শুধুমাত্র সংখ্যাগত গণনা দ্বারা দেওয়া হয়, বিশ্লেষণাত্মক অভিব্যক্তির অভাব

প্রভাব

  1. একাডেমিক অবদান: ক্ষণস্থায়ী রড-আরোহণ তত্ত্বের শূন্যতা পূরণ করে
  2. পদ্ধতি মূল্য: প্রদত্ত বিশ্লেষণ পদ্ধতি অন্যান্য জটিল তরল সমস্যায় সম্প্রসারণযোগ্য
  3. প্রকৌশল প্রয়োগ: রিওমিটার ডিজাইন এবং উপাদান বৈশিষ্ট্যকরণের জন্য তাত্ত্বিক সমর্থন প্রদান করে

প্রযোজ্য পরিস্থিতি

  1. পাতলা পলিমার দ্রবণ: Boger তরল ইত্যাদি ধ্রুবক শিয়ার ভিসকোসিটি viscoelastic তরল
  2. কম ঘূর্ণন গতি: ছোট Weissenberg সংখ্যার পরীক্ষামূলক শর্ত
  3. রিওলজি পরিমাপ: স্বাভাবিক চাপ সহগ এবং শিথিলকরণ সময়ের পরিমাপ

তথ্যসূত্র

  1. Joseph, D.D. & Fosdick, R.L. (1973). The free surface on a liquid between cylinders rotating at different speeds. Arch. Rat. Mech. Anal., 49, 321-380.
  2. Oldroyd, J.G. (1950). On the formulation of rheological equations of state. Proc. R. Soc. A, 200, 523-541.
  3. Bird, R.B., Armstrong, R.C. & Hassager, O. (1987). Dynamics of Polymeric Liquids, Vol. 1, 2nd edn. John Wiley and Sons.

সামগ্রিক মূল্যায়ন: এটি একটি উচ্চ মানের তাত্ত্বিক তরল গতিবিদ্যা পেপার যা প্রথমবারের মতো Oldroyd-B তরলে রড-আরোহণ ঘটনার ক্ষণস্থায়ী আচরণ সিস্টেমেটিকভাবে বিশ্লেষণ করে। পেপারটির গাণিতিক অনুমান কঠোর, ভৌত অন্তর্দৃষ্টি গভীর এবং পদ্ধতি উদ্ভাবনী, জটিল তরল তত্ত্ব এবং রিওলজি প্রয়োগে গুরুত্বপূর্ণ অবদান রাখে। যদিও পরামিতি পরিসীমা সীমাবদ্ধতা রয়েছে, তবে এর তাত্ত্বিক কাঠামো এবং বিশ্লেষণ পদ্ধতির ভালো সম্প্রসারণ মূল্য রয়েছে।