2025-11-14T22:25:11.017657

Lattice Boltzmann model for non-ideal compressible fluid dynamics

Hosseini, Feinberg, Karlin
We present a lattice Boltzmann formulation for the simulation of compressible, non-ideal fluid flows. The method employs first-neighbor lattices and introduces a consistent set of correction terms through quasi-equilibrium attractors, ensuring positive-definite and Galilean-invariant Navier-Stokes dissipation rates. This construction circumvents the need for extended stencils or ad hoc regularization, while maintaining numerical stability and thermodynamic consistency across a broad range of flow regimes. The resulting model accurately reproduces both Euler- and Navier-Stokes-level hydrodynamics. As a stringent validation, we demonstrate, for the first time within a lattice Boltzmann framework, quantitatively accurate simulations of drop-shock interactions at Mach numbers up to 1.47. The proposed approach thus extends the applicability of lattice Boltzmann methods to high-speed, non-ideal compressible flows with a minimal kinetic stencil.
academic

অ-আদর্শ সংকোচনযোগ্য তরল গতিশীলতার জন্য জালি বোল্টজম্যান মডেল

মৌলিক তথ্য

  • পেপার আইডি: 2510.14712
  • শিরোনাম: অ-আদর্শ সংকোচনযোগ্য তরল গতিশীলতার জন্য জালি বোল্টজম্যান মডেল
  • লেখক: এস. এ. হোসেইনি, এম. ফেইনবার্গ, আই. ভি. কার্লিন (ইটিএইচ জুরিখ)
  • শ্রেণীবিভাগ: physics.flu-dyn
  • প্রকাশিত জার্নাল: জার্নাল অফ ফ্লুইড মেকানিক্স
  • পেপার লিঙ্ক: https://arxiv.org/abs/2510.14712

সারসংক্ষেপ

এই পেপারটি সংকোচনযোগ্য অ-আদর্শ তরল প্রবাহ অনুকরণের জন্য একটি জালি বোল্টজম্যান পদ্ধতি প্রস্তাব করে। এই পদ্ধতিটি প্রথম নিকটতম প্রতিবেশী জালি ব্যবহার করে, প্রস্তুত-সাম্যাবস্থা আকর্ষক মাধ্যমে সামঞ্জস্যপূর্ণ সংশোধন পদের একটি সেট প্রবর্তন করে, যা ইতিবাচক-নির্দিষ্টতা এবং গ্যালিলিয়ান অপরিবর্তনীয় নাভিয়ার-স্টোকস অপচয় হার নিশ্চিত করে। এই নির্মাণটি সম্প্রসারিত টেমপ্লেট বা অস্থায়ী নিয়মিতকরণের প্রয়োজনীয়তা এড়ায়, একই সাথে প্রবাহ অবস্থার বিস্তৃত পরিসরে সংখ্যাগত স্থিতিশীলতা এবং তাপগতিগত সামঞ্জস্য বজায় রাখে। ফলস্বরূপ মডেলটি অয়লার এবং নাভিয়ার-স্টোকস স্তরের তরল গতিশীলতা সঠিকভাবে পুনরুৎপাদন করে। কঠোর যাচাইকরণ হিসাবে, লেখকরা প্রথমবারের মতো জালি বোল্টজম্যান কাঠামোর মধ্যে ১.৪৭ মাখ সংখ্যা পর্যন্ত ফোঁটা-শক ইন্টারঅ্যাকশনের পরিমাণগতভাবে নির্ভুল অনুকরণ প্রদর্শন করেছেন।

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

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

অ-আদর্শ সংকোচনযোগ্য তরল গতিশীলতা একটি উদীয়মান এবং দ্রুত বিকশিত তরল গতিশীলতার শাখা, প্রধানত নিকট-সংকটীয়, অতিক্রম-সংকটীয় এবং অতি-সংকটীয় অবস্থায় কাজ করা পদ্ধতি এবং প্রযুক্তির উদ্ভবের কারণে। এই ধরনের তরলগুলির মধ্যে রয়েছে:

  1. অতি-সংকটীয় তরল: সংকটীয় বিন্দুর উপরে চাপ এবং তাপমাত্রার অবস্থায় তরল
  2. ঘন বাষ্প: সংকোচন ফ্যাক্টর যা আদর্শ গ্যাস আচরণ থেকে বিচ্যুত হয়
  3. বিজেডটি তরল: নেতিবাচক মৌলিক গ্যাস গতিশীলতা ডেরিভেটিভ সহ উচ্চ-পলিমার জটিলতার তরল

গুরুত্ব

  1. শিল্প প্রয়োগের চাহিদা: জৈব র‍্যাঙ্কাইন চক্র, অতি-সংকটীয় CO₂ টারবাইন এবং অন্যান্য শক্তি প্রযুক্তি
  2. পরীক্ষামূলক ডেটার স্বল্পতা: পরীক্ষামূলক সেটআপ জটিল, ডেটা অর্জন কঠিন
  3. সংখ্যাগত সরঞ্জামের চাহিদা: জটিল ভৌত ঘটনা বোঝার জন্য সামঞ্জস্যপূর্ণ এবং দক্ষ সংখ্যাগত সরঞ্জামের প্রয়োজন

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

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

মূল অবদান

  1. নতুন দ্বি-বিতরণ ফাংশন এলবিএম মডেল প্রস্তাব: ক্লাসিক্যাল প্রথম নিকটতম প্রতিবেশী জালি (D3Q27) ব্যবহার করে অ-আদর্শ সংকোচনযোগ্য প্রবাহ অনুকরণ বাস্তবায়ন
  2. প্রস্তুত-সাম্যাবস্থা আকর্ষক প্রক্রিয়া প্রবর্তন: fᵢ এবং gᵢ সংশোধন পদের মাধ্যমে ইতিবাচক-নির্দিষ্টতা এবং গ্যালিলিয়ান অপরিবর্তনীয় অপচয় হার নিশ্চিত করা
  3. স্বাধীন আয়তন সান্দ্রতা নিয়ন্ত্রণ বাস্তবায়ন: বিজিকে কাঠামোতে আয়তন সান্দ্রতা সম্ভাব্যত নেতিবাচক হওয়ার সমস্যা সমাধান
  4. প্রথমবারের মতো উচ্চ মাখ সংখ্যা ফোঁটা-শক ইন্টারঅ্যাকশন অনুকরণ বাস্তবায়ন: এলবিএম কাঠামোর মধ্যে Ma=1.47 এর পরিমাণগতভাবে নির্ভুল অনুকরণ অর্জন
  5. তাপগতিগত সামঞ্জস্য বজায় রাখা: প্রবাহ অবস্থার বিস্তৃত পরিসরে সংখ্যাগত স্থিতিশীলতা

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

কাজের সংজ্ঞা

নাভিয়ার-স্টোকস-কর্টেওয়েগ (এনএসকে) সমীকরণ সেট সন্তুষ্ট করে এমন সংকোচনযোগ্য অ-আদর্শ তরল অনুকরণ:

  • ভর সংরক্ষণ: ∂ₜρ + ∇·(ρu) = 0
  • গতিবেগ সংরক্ষণ: ∂ₜ(ρu) + ∇·(ρu⊗u) + ∇P + ∇·T^NS + ∇·T^K = 0
  • শক্তি সংরক্ষণ: ∂ₜ(ρE) + ∇·((ρE+P)u) + u·(∇·T^K) + ∇·(u·T^NS) + ∇·q = 0

যেখানে T^K হল কর্টেওয়েগ পৃষ্ঠ টেনশন টেনসর, T^NS হল সান্দ্র চাপ টেনসর।

মডেল স্থাপত্য

১. দ্বি-বিতরণ ফাংশন ডিজাইন

দুটি বিতরণ ফাংশন সেট {fᵢ, gᵢ}, i=1,...,27 ব্যবহার করা হয়:

  • fᵢ: ভর এবং গতিবেগ সংরক্ষণ পরিচালনা করে
  • gᵢ: শক্তি সংরক্ষণ পরিচালনা করে

২. জালি বোল্টজম্যান সমীকরণ

fᵢ(x+cᵢδt, t+δt) = fᵢ + 2β(fᵢᵉᵠ - fᵢ) + (1-β)(fᵢ* - fᵢᵉᵠ)
gᵢ(x+cᵢδt, t+δt) = gᵢ + 2β(gᵢᵉᵠ - gᵢ) + (1-β)(gᵢ* - gᵢᵉᵠ)

যেখানে β = Pδt/(2μ + Pδt) শিথিলকরণ প্যারামিটার।

৩. সাম্যাবস্থা বিতরণ ফাংশন

পণ্য ফর্ম গ্রহণ করা হয়:

fᵢᵉᵠ = ρ ∏_α Ψᵢα(uα, P/ρ + uα²)

যেখানে Ψᵢα হল মৌলিক ফাংশন:

Ψᵢα(ξα, ζαα) = 1 - cᵢα² + ½[(3cᵢα² - 2)ζαα + cᵢαξα]

৪. প্রস্তুত-সাম্যাবস্থা সংশোধন পদ

f বিতরণ ফাংশন সংশোধন:

  • বল পদ: Fα = -κρ∂α∇²ρ (কর্টেওয়েগ বল)
  • চাপ সংশোধন: Φαα কর্ণ চাপ টেনসর সংশোধন এবং স্বাধীন আয়তন সান্দ্রতা নিয়ন্ত্রণ অন্তর্ভুক্ত করে

g বিতরণ ফাংশন সংশোধন:

  • স্থানান্তরিত বেগ: u* = u + δt F/ρ
  • স্থানান্তরিত তাপমাত্রা: T* = T - δt²F·F/(2ρ²cᵥ)
  • তাপ প্রবাহ সংশোধন: অ-আদর্শ গ্যাসের এনথ্যালপি গ্রেডিয়েন্ট প্রভাব বিবেচনা করা

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

১. প্রস্তুত-সাম্যাবস্থা আকর্ষক প্রক্রিয়া

fᵢ এবং gᵢ প্রবর্তনের মাধ্যমে, সিস্টেম সহজভাবে স্থানীয় সাম্যাবস্থার দিকে শিথিল হয় না, বরং সংশোধিত প্রস্তুত-সাম্যাবস্থা অবস্থার দিকে শিথিল হয়, যা নিশ্চিত করে:

  • সঠিক ম্যাক্রোস্কোপিক গতিশীলতা
  • সংখ্যাগত স্থিতিশীলতা
  • গ্যালিলিয়ান অপরিবর্তনীয়তা

২. স্বাধীন আয়তন সান্দ্রতা নিয়ন্ত্রণ

ঐতিহ্যবাহী বিজিকে মডেলে আয়তন সান্দ্রতা η' = ((D+2)/D - ρcs²/P)μ নেতিবাচক হতে পারে। সংশোধন পদ Φ' এর মাধ্যমে স্বাধীন নিয়ন্ত্রণ বাস্তবায়ন:

Φ' = P((D+2)/D - ρcs²/P - η/μ)(∇·u)

৩. তাপগতিগত সামঞ্জস্য

সাবধানে ডিজাইন করা স্থানান্তরিত তাপমাত্রা T* এবং শক্তি সংশোধনের মাধ্যমে, অ-রৈখিক অবস্থা সমীকরণের অধীনে তাপগতিগত সামঞ্জস্য নিশ্চিত করা হয়।

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

যাচাইকরণ কেস

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

পদার্থ প্যারামিটার

প্রধানত নাইট্রোজেন (N₂) ব্যবহার করে যাচাইকরণ:

  • সংকটীয় চাপ: Pc = 3.4×10⁶ Pa
  • সংকটীয় ঘনত্ব: ρc = 241.96 kg/m³
  • সংকটীয় তাপমাত্রা: Tc = 126.2 K
  • R/cv = 0.4

অবস্থা সমীকরণ

ভ্যান ডার ওয়ালস অবস্থা সমীকরণ গ্রহণ করা হয়:

P(ρ,T) = ρRT/(1-bρ) - aρ²

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

প্রধান ফলাফল

১. মৌলিক ভৌত পরিমাণ যাচাইকরণ

  • শব্দের গতি: সম্পৃক্ত তরল এবং বাষ্প শাখায়, Tr∈0.7,1 পরিসরে বিশ্লেষণাত্মক সমাধানের সাথে নিখুঁত সামঞ্জস্য
  • শিয়ার সান্দ্রতা: বিভিন্ন মাখ সংখ্যায় গ্যালিলিয়ান অপরিবর্তনীয়তা বজায় রাখে, ত্রুটি <1%
  • তাপীয় কুয়েট প্রবাহ: Pr∈{0.6,1.2,4.9} এবং Ma∈{0.8,1.2,1.6} এ বিশ্লেষণাত্মক সমাধানের সাথে সামঞ্জস্যপূর্ণ

২. বহু-পর্যায়ের প্রবাহ বৈশিষ্ট্য

  • সহাবস্থান ঘনত্ব: Tr∈0.3,0.99 পরিসরে ম্যাক্সওয়েল সমান-ক্ষেত্র নির্মাণের সাথে সম্পূর্ণ সামঞ্জস্য
  • ইন্টারফেস সংমিশ্রণ: বিভিন্ন জালি রেজোলিউশনে ভাল সংমিশ্রণ বৈশিষ্ট্য প্রদর্শন

३. অ-আদর্শ শক ট্যুব

তিনটি ভিন্ন অ-আদর্শ শক ট্যুব কনফিগারেশন সফলভাবে অনুকরণ:

  • কেস I: Γ<0 অঞ্চলে বিরল শক
  • কেস II: ক্লাসিক্যাল সংকোচন তরঙ্গ আচরণ
  • কেস III: গতিশীল Γ=0 লাইন অতিক্রম করা মিশ্র বিরল তরঙ্গ

সমস্ত কেস রেফারেন্স ডেটার সাথে চমৎকার সামঞ্জস্য।

४. শক-ফোঁটা ইন্টারঅ্যাকশন

  • এলবিএম কাঠামোর মধ্যে প্রথমবারের মতো Ma=1.47 এর উচ্চ মাখ সংখ্যা অনুকরণ বাস্তবায়ন
  • আপতিত শক, সংক্রমিত তরঙ্গ, প্রতিফলিত শক, মাখ রড ইত্যাদি জটিল তরঙ্গ কাঠামো সঠিকভাবে ক্যাপচার করা
  • ফোঁটা বিকৃতি বিবর্তন পরীক্ষামূলক ডেটার সাথে উচ্চ সামঞ্জস্য

বিলোপন পরীক্ষা

বিভিন্ন জালি রেজোলিউশনের ইন্টারফেস পরীক্ষার মাধ্যমে মডেলের সংমিশ্রণ যাচাই করা হয়েছে:

  • δx থেকে 5μm থেকে 0.1μm এর চারটি রেজোলিউশন স্তর
  • ভাল জালি সংমিশ্রণ বৈশিষ্ট্য প্রদর্শন করে

মূল আবিষ্কার

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

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

ঐতিহ্যবাহী এলবিএম বহু-পর্যায়ের প্রবাহ মডেল

  1. রঙ গ্রেডিয়েন্ট মডেল (গানস্টেনসেন এট আল. 1991)
  2. সিউডো-সম্ভাব্যতা মডেল (শান এবং চেন 1993)
  3. মুক্ত শক্তি মডেল (সুইফট এট আল. 1996)

সংকোচনযোগ্য এলবিএম উন্নয়ন

  • হে এবং ডোলেন (2002): সংকোচনযোগ্য অ-আদর্শ প্রবাহের প্রথম প্রচেষ্টা
  • বিদ্যমান মডেল প্রধানত ফুটন্ত অ্যাপ্লিকেশনের জন্য
  • সংকোচনযোগ্য অ-আদর্শ প্রবাহ মডেলের সিস্টেমেটিক অভাব

এই পেপারের আপেক্ষিক সুবিধা

  1. সিস্টেমেটিক পদ্ধতি: সম্পূর্ণ দ্বি-বিতরণ ফাংশন কাঠামো
  2. তাত্ত্বিক ভিত্তি: কঠোর বহু-স্কেল বিশ্লেষণ
  3. যাচাইকরণ সম্পূর্ণতা: মৌলিক ভৌত পরিমাণ থেকে জটিল প্রবাহ ঘটনার সম্পূর্ণ যাচাইকরণ

উপসংহার এবং আলোচনা

প্রধান উপসংহার

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

সীমাবদ্ধতা

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

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

  1. আরও শক্তিশালী সংঘর্ষ মডেলে সম্প্রসারণ: অশান্ত কনফিগারেশন অনুকরণের জন্য
  2. আরও জটিল অবস্থা সমীকরণ: যেমন পেং-রবিনসন ইত্যাদি
  3. বড় আকারের ত্রিমাত্রিক প্রয়োগ: প্রকৃত প্রকৌশল সমস্যার প্রয়োগ

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

সুবিধা

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

অপূর্ণতা

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

প্রভাব

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

প্রযোজ্য দৃশ্যকল্প

  1. অতি-সংকটীয় তরল প্রযুক্তি: CO₂ টারবাইন, জৈব র‍্যাঙ্কাইন চক্র
  2. বহু-পর্যায়ের প্রবাহ: ফোঁটা ভাঙ্গন, পর্যায় পরিবর্তন প্রক্রিয়া
  3. শক গতিশীলতা: উচ্চ-গতির প্রবাহে জটিল তরঙ্গ ইন্টারঅ্যাকশন
  4. সরাসরি সংখ্যাগত অনুকরণ: উচ্চ নির্ভুলতা প্রয়োজন এমন মৌলিক ভৌত গবেষণা

তথ্যসূত্র

  1. গার্ডোন, এ. এট আল. (2024)। অ-আদর্শ সংকোচনযোগ্য তরল গতিশীলতা ঘন বাষ্প এবং অতি-সংকটীয় তরলের। বার্ষিক পর্যালোচনা তরল মেকানিক্স, 56, 241-269।
  2. হোসেইনি, এস.এ. এবং কার্লিন, আই.ভি. (2023)। অ-আদর্শ তরলের জন্য জালি বোল্টজম্যান: মৌলিক বিষয় এবং অনুশীলন। পদার্থবিজ্ঞান প্রতিবেদন, 1030, 1-137।
  3. শান, এক্স. এবং চেন, এইচ. (1993)। একাধিক পর্যায় এবং উপাদান সহ প্রবাহ অনুকরণের জন্য জালি বোল্টজম্যান মডেল। শারীরিক পর্যালোচনা ই, 47(3), 1815।