While recent approaches, such as the macroscopic forcing method (MFM) or Green's function-based approaches, can be used to compute Reynolds-averaged Navier--Stokes closure operators using forced direct numerical simulations, MFM can also be used to directly compute moments of the effective nonlocal and anisotropic eddy diffusivities. The low-order spatial and temporal moments contain limited information about the eddy diffusivity but are often sufficient for quantification and modeling of nonlocal and anisotropic effects. However, when using MFM to compute eddy diffusivity moments, the statistical convergence can be slow for higher-order moments. In this work, we demonstrate that using the same direct numerical simulation (DNS) for all forced MFM simulations improves statistical convergence of the eddy diffusivity moments. We present its implementation in conjunction with a decomposition method that handles the MFM forcing semi-analytically and allows for consistent boundary condition treatment, which we develop for both scalar and momentum transport. We demonstrate that for a two-dimensional Rayleigh--Taylor instability case study, using the same DNS for all forced MFM simulations results in convergence with O(100) simulations rather than O(1000) simulations. We then demonstrate the impacts of improved convergence on the quantification of the eddy diffusivity.
- পেপার আইডি: 2503.06418
- শিরোনাম: এডি ডিফিউসিভিটি মোমেন্টস পরিমাপে উন্নত পরিসংখ্যানগত সংমিশ্রণের কৌশল
- লেখক: ডানা এল. ও.-এল. লাভাকট (স্ট্যানফোর্ড বিশ্ববিদ্যালয়), জেসি লিউ (স্ট্যানফোর্ড বিশ্ববিদ্যালয়), ব্র্যান্ডন ই. মরগান (লরেন্স লিভারমোর ন্যাশনাল ল্যাবরেটরি), আলি মানি (স্ট্যানফোর্ড বিশ্ববিদ্যালয়)
- শ্রেণীবিভাগ: physics.flu-dyn
- প্রকাশনার সময়: অক্টোবর ১৫, ২০২৫ (arXiv প্রি-প্রিন্ট)
- পেপার লিঙ্ক: https://arxiv.org/abs/2503.06418v2
এই পেপারটি ম্যাক্রোস্কোপিক ফোর্সিং মেথড (MFM) এ এডি ডিফিউসিভিটি কোএফিশিয়েন্টের মোমেন্টস গণনার পরিসংখ্যানগত সংমিশ্রণ সমস্যার সমাধানের জন্য উন্নত কৌশল প্রস্তাব করে। গবেষণায় দেখা যায় যে সমস্ত ফোর্সিং MFM সিমুলেশনের জন্য একই ডাইরেক্ট নিউমেরিক্যাল সিমুলেশন (DNS) ব্যবহার করা এডি ডিফিউসিভিটি মোমেন্টসের পরিসংখ্যানগত সংমিশ্রণকে উল্লেখযোগ্যভাবে উন্নত করতে পারে। লেখকরা ডিকম্পোজিশন মেথডের সাথে সংমিশ্রিত একটি বাস্তবায়ন স্কিম তৈরি করেছেন যা MFM ফোর্সিং অর্ধ-বিশ্লেষণাত্মকভাবে পরিচালনা করতে এবং সীমানা শর্তের সামঞ্জস্যপূর্ণ পরিচালনা নিশ্চিত করতে পারে। দ্বিমাত্রিক Rayleigh-Taylor অস্থিতিশীলতা কেস স্টাডিতে, এই পদ্ধতি সংমিশ্রণের জন্য প্রয়োজনীয় সিমুলেশনের সংখ্যা O(1000) থেকে O(100) এ হ্রাস করে।
টার্বুলেন্ট মডেলিংয়ে, রেনল্ডস গড় Navier-Stokes (RANS) সমীকরণগুলি অবিবৃত পদগুলি যেমন স্কেলার ফ্লাক্স বা রেনল্ডস স্ট্রেস পরিচালনা করার জন্য ক্লোজার মডেলের প্রয়োজন। ঐতিহ্যবাহী Boussinesq অনুমান স্থানীয় এবং আইসোট্রপিক এডি ডিফিউসিভিটি কোএফিশিয়েন্ট অনুমান করে, কিন্তু জটিল প্রবাহের জন্য প্রায়শই অকার্যকর।
- অ-স্থানীয় এবং অ্যানিসোট্রপিক প্রভাব: প্রকৃত এডি ডিফিউসিভিটি কোএফিশিয়েন্ট সাধারণত অ-স্থানীয় এবং অ্যানিসোট্রপিক, যার জন্য আরও নির্ভুল বর্ণনা প্রয়োজন
- ধীর পরিসংখ্যানগত সংমিশ্রণ: MFM ব্যবহার করে উচ্চ-ক্রম এডি ডিফিউসিভিটি মোমেন্টস গণনা করার সময়, পরিসংখ্যানগত সংমিশ্রণ দুর্বল
- উচ্চ গণনামূলক খরচ: ঐতিহ্যবাহী পদ্ধতিগুলি পরিসংখ্যানগত সংমিশ্রণ অর্জনের জন্য প্রচুর DNS সিমুলেশনের প্রয়োজন
বিদ্যমান ম্যাক্রোস্কোপিক ফোর্সিং মেথড (MFM) এবং গ্রীন ফাংশন পদ্ধতি অ-স্থানীয় এবং অ্যানিসোট্রপিক এডি ডিফিউসিভিটি কোএফিশিয়েন্ট গণনা করতে পারে, কিন্তু উচ্চ-ক্রম মোমেন্টস গণনায় ধীর পরিসংখ্যানগত সংমিশ্রণ সমস্যা রয়েছে। এই পেপারটি পরিসংখ্যানগত সংমিশ্রণ ত্বরান্বিত করতে এবং গণনামূলক খরচ হ্রাস করার জন্য উন্নত কৌশল বিকাশের লক্ষ্য রাখে।
- পরিসংখ্যানগত ত্রুটি বৃদ্ধি সমস্যা চিহ্নিত এবং সমাধান: স্বাধীন ডোনর সিমুলেশন ব্যবহার করলে উচ্চ-ক্রম মোমেন্টসের পরিসংখ্যানগত ত্রুটি বৃদ্ধি পায় এবং একক ডোনর ব্যবহারের সমাধান প্রস্তাব করে
- ডিকম্পোজিশন MFM পদ্ধতি বিকাশ: সাধারণ অ-স্থির এবং বিশৃঙ্খল প্রবাহ এবং মোমেন্টাম পরিবহনের জন্য ডিকম্পোজিশন পদ্ধতি সম্প্রসারিত করা
- গণনামূলক দক্ষতা উল্লেখযোগ্যভাবে উন্নত করা: 2D Rayleigh-Taylor অস্থিতিশীলতায় প্রায় 10 গুণ গণনামূলক ত্বরণ অর্জন করা
- সম্পূর্ণ তাত্ত্বিক এবং বাস্তবায়ন ফ্রেমওয়ার্ক প্রদান: স্কেলার এবং মোমেন্টাম পরিবহনের জন্য সম্পূর্ণ গাণিতিক ফ্রেমওয়ার্ক এবং বাস্তবায়ন বিবরণ প্রদান করা
এই পেপারটি অ-স্থানীয় এবং অ্যানিসোট্রপিক এডি ডিফিউসিভিটি কোএফিশিয়েন্টের স্থানকালীন মোমেন্টস দক্ষতার সাথে কীভাবে গণনা করতে হয় তা অধ্যয়ন করে, যেখানে এই মোমেন্টসগুলি এডি ডিফিউসিভিটি কোএফিশিয়েন্ট কার্নেলের সীমিত কিন্তু গুরুত্বপূর্ণ তথ্য ধারণ করে, যা অ-স্থানীয় এবং অ্যানিসোট্রপিক প্রভাব পরিমাপ এবং মডেল করার জন্য যথেষ্ট।
এডি ডিফিউসিভিটি কোএফিশিয়েন্ট Kramers-Moyal সম্প্রসারণের মাধ্যমে মোমেন্টসের আকারে প্রকাশ করা যায়:
−⟨ui′c′⟩(x,t)=[Dij00(x,t)+Dijk10(x,t)∂xk∂+⋯+Dij01(x,t)∂t∂+⋯]∂xj∂⟨c⟩
যেখানে বিভিন্ন ক্রমের মোমেন্টস নিম্নরূপ সংজ্ঞায়িত:
- Dij00(x,t)=∫∫Dij(x,x′,t,t′)dx′dt′ (শূন্য-ক্রম স্থানকালীন মোমেন্ট)
- Dijk10(x,t)=∫∫(xk′−xk)Dij(x,x′,t,t′)dx′dt′ (প্রথম-ক্রম স্থানীয় মোমেন্ট)
- Dij01(x,t)=∫∫(t′−t)Dij(x,x′,t,t′)dx′dt′ (প্রথম-ক্রম সময়ীয় মোমেন্ট)
MFM স্কেলার পরিবহন সমীকরণে ম্যাক্রোস্কোপিক ফোর্সিং যোগ করে এডি ডিফিউসিভিটি কোএফিশিয়েন্টের মোমেন্টস গণনা করে:
∂t∂c+∂xi∂(uic)=DM∂xi∂xi∂2c+s
বিভিন্ন গড় স্কেলার গ্রেডিয়েন্ট নির্দিষ্ট করে (যেমন বহুপদী আকার), সংশ্লিষ্ট এডি ডিফিউসিভিটি কোএফিশিয়েন্টের মোমেন্টস সরাসরি পাওয়া যায়।
ডিকম্পোজিশন পদ্ধতি স্কেলার ওঠানামা c′ কে সম্প্রসারিত করে:
c′(x,t)=[cj00(x,t)+cjk10(x,t)∂xk∂+⋯+cj01(x,t)∂t∂+⋯]∂xj∂⟨c⟩
এটি MFM ফোর্সিং অর্ধ-বিশ্লেষণাত্মকভাবে পরিচালনা করতে এবং পর্যায়ক্রমিক সীমানা শর্তের সামঞ্জস্যতা সমস্যা সমাধান করতে অনুমতি দেয়।
- একক ডোনর কৌশল: সমস্ত রিসিভার সমীকরণে প্রবাহ ক্ষেত্র সরবরাহ করার জন্য একই DNS সিমুলেশন ব্যবহার করে, স্বাধীন ডোনর সিমুলেশনের মধ্যে পার্থক্য থেকে উদ্ভূত ত্রুটি বৃদ্ধি এড়ানো
- ত্রুটি প্রচার বিশ্লেষণ: তাত্ত্বিক বিশ্লেষণ দেখায় যে স্বাধীন ডোনর ব্যবহার করলে, উচ্চ-ক্রম মোমেন্টসের পরিসংখ্যানগত ত্রুটি স্থানীয় স্থানাঙ্কের শক্তি অনুযায়ী বৃদ্ধি পায়:
D1=F1−x1F0+x1(ε01−ε00)+ε11
- ডিকম্পোজিশন পদ্ধতির সম্প্রসারণ: মূলত স্থির স্তরীয় প্রবাহের জন্য ব্যবহৃত ডিকম্পোজিশন পদ্ধতি সাধারণ অ-স্থির বিশৃঙ্খল প্রবাহ এবং মোমেন্টাম পরিবহনে সম্প্রসারিত করা
- গণনামূলক ডোমেইন: 2049×2049 গ্রিডের বর্গাকার ডোমেইন
- সীমানা শর্ত: x₁ দিকে পর্যায়ক্রমিক, x₂ দিকে নো-স্লিপ এবং নো-পেনিট্রেশন
- ভৌত পরামিতি:
- Atwood সংখ্যা: A = 0.05
- সর্বোচ্চ Mach সংখ্যা: Ma_max = 0.05
- Grashof সংখ্যা: Gr = 1
- Schmidt সংখ্যা: Sc = 1
- রেনল্ডস সংখ্যা: Re_T = 54, Re_L = 8,000
- পরিসংখ্যানগত সংমিশ্রণ: বিভিন্ন বাস্তবায়ন সংখ্যার অধীনে মূল-গড়-বর্গ ত্রুটি (RMSE) দ্বারা মূল্যায়ন করা
- গণনামূলক দক্ষতা: একই সংমিশ্রণ নির্ভুলতা অর্জনের জন্য প্রয়োজনীয় সিমুলেশনের সংখ্যা তুলনা করা
- মডেলিং কর্মক্ষমতা: মোমেন্ট ম্যাচিং ইনভার্সন (MMI) পদ্ধতির মাধ্যমে নির্মিত ক্লোজার অপারেটরের কর্মক্ষমতা মূল্যায়ন করা
- মান MFM (স্বাধীন ডোনর): প্রতিটি মোমেন্টের জন্য স্বাধীন ডোনর-রিসিভার জোড়া ব্যবহার করা
- ডিকম্পোজিশন MFM (একক ডোনর): সমস্ত রিসিভার সমীকরণ একই ডোনর সিমুলেশন ব্যবহার করা
একক ডোনর ব্যবহার করে ডিকম্পোজিশন MFM স্বাধীন ডোনরের মান MFM এর তুলনায় পরিসংখ্যানগত সংমিশ্রণে উল্লেখযোগ্য উন্নতি প্রদর্শন করে:
- শূন্য-ক্রম মোমেন্ট D⁰⁰: উভয় পদ্ধতি একই রকম কর্মক্ষমতা প্রদর্শন করে, কারণ গণনা পদ্ধতি গাণিতিকভাবে সমতুল্য
- প্রথম-ক্রম সময়ীয় মোমেন্ট D⁰¹: ডিকম্পোজিশন MFM 100টি বাস্তবায়নে ভাল সংমিশ্রণ অর্জন করে, যখন মান MFM 1000টি বাস্তবায়নের প্রয়োজন
- প্রথম-ক্রম স্থানীয় মোমেন্ট D¹⁰: ডিকম্পোজিশন MFM 100টি বাস্তবায়নে মান MFM এর তুলনায় স্পষ্টভাবে কম ত্রুটি প্রদর্শন করে
- ত্বরণ অনুপাত: প্রায় 10 গুণ গণনামূলক ত্বরণ (O(1000) থেকে O(100) সিমুলেশনে হ্রাস)
- ত্রুটি স্কেলিং: তাত্ত্বিক বিশ্লেষণ এবং সংখ্যাগত ফলাফল উভয়ই সময় এবং স্থানীয় স্থানাঙ্কের সাথে ত্রুটি বৃদ্ধির ঘটনা নিশ্চিত করে
Re_τ = 180 এর টার্বুলেন্ট চ্যানেল প্রবাহে পদ্ধতির কার্যকারিতা যাচাই করা:
- মান MFM (একক ডোনর) এবং ডিকম্পোজিশন MFM ফলাফল প্রায় সম্পূর্ণভাবে সামঞ্জস্যপূর্ণ
- উন্নতি প্রধানত ডোনর একীকরণ থেকে আসে, ডিকম্পোজিশন পদ্ধতি থেকে নয়
উন্নত সংমিশ্রিত এডি ডিফিউসিভিটি কোএফিশিয়েন্টের মোমেন্টস ব্যবহার করে নির্মিত MMI ক্লোজার অপারেটর:
- ডিকম্পোজিশন MFM: 200টি বাস্তবায়নে নির্মিত মডেল DNS ফলাফলের সাথে ভাল সামঞ্জস্য প্রদর্শন করে
- মান MFM: একই বাস্তবায়ন সংখ্যায় নির্মিত মডেল উল্লেখযোগ্য বিচ্যুতি প্রদর্শন করে, যা ভুল মডেলিং সিদ্ধান্তের দিকে পরিচালিত করতে পারে
- গ্রীন ফাংশন পদ্ধতি: Hamba ইত্যাদি দ্বারা বিকশিত নির্ভুল কিন্তু গণনামূলকভাবে ব্যয়বহুল পদ্ধতি
- দ্রুত MFM: Bryngelson ইত্যাদি দ্বারা বিকশিত লুকানো বিরলতা ব্যবহার করে ত্বরণ পদ্ধতি
- সহযোগী MFM: Liu ইত্যাদি দ্বারা প্রস্তাবিত সহযোগী পদ্ধতির উপর ভিত্তি করে লক্ষ্য গণনা কৌশল
- Kraichnan এর নির্ভুল অ-স্থানীয় অভিব্যক্তি
- বিস্তার মডেলিংয়ে মোমেন্ট-গ্রেডিয়েন্ট পদ্ধতি
- টার্বুলেন্ট মডেলিংয়ে Kramers-Moyal সম্প্রসারণ
- পরিসংখ্যানগত ত্রুটির উৎস চিহ্নিতকরণ: স্বাধীন ডোনর সিমুলেশনের মধ্যে ক্ষুদ্র পার্থক্য উচ্চ-ক্রম মোমেন্টস গণনায় O(1) এর ত্রুটিতে বৃদ্ধি পায় তা প্রমাণ করা
- কার্যকর সমাধান পদ্ধতি: একক ডোনর ব্যবহার করা ত্রুটি বৃদ্ধি দূর করতে পারে এবং পরিসংখ্যানগত সংমিশ্রণ উল্লেখযোগ্যভাবে উন্নত করতে পারে
- পদ্ধতির সার্বজনীনতা: ডিকম্পোজিশন MFM স্কেলার এবং মোমেন্টাম পরিবহনে প্রযোজ্য, পর্যায়ক্রমিক সীমানা শর্ত সমস্যা পরিচালনা করতে পারে
- প্রবাহ প্রকার সীমাবদ্ধতা: প্রধানত 2D Rayleigh-Taylor অস্থিতিশীলতায় যাচাই করা, আরও প্রবাহ প্রকারের পরীক্ষার প্রয়োজন
- সংখ্যাগত বাস্তবায়ন নির্ভরতা: উন্নতির মাত্রা নির্দিষ্ট সংখ্যাগত সমাধানকারী এবং সমান্তরাল বাস্তবায়নের উপর নির্ভর করতে পারে
- তাত্ত্বিক বিশ্লেষণ সরলীকরণ: ত্রুটি বিশ্লেষণ সরলীকৃত অনুমানের উপর ভিত্তি করে, প্রকৃত পরিস্থিতি আরও জটিল হতে পারে
- প্রয়োগের সম্প্রসারণ: আরও জটিল প্রবাহে পদ্ধতির কার্যকারিতা যাচাই করা
- তাত্ত্বিক উন্নতি: আরও সম্পূর্ণ ত্রুটি প্রচার তত্ত্ব বিকাশ করা
- অ্যালগরিদম অপ্টিমাইজেশন: অন্যান্য ত্বরণ কৌশলের সাথে সংমিশ্রণ করে গণনামূলক দক্ষতা আরও উন্নত করা
- সমস্যা চিহ্নিতকরণ নির্ভুল: MFM এ পরিসংখ্যানগত ত্রুটি বৃদ্ধি সমস্যা নির্ভুলভাবে চিহ্নিত করা
- সমাধান পদ্ধতি কার্যকর: প্রস্তাবিত একক ডোনর কৌশল সহজ এবং কার্যকর
- তাত্ত্বিক বিশ্লেষণ গভীর: ত্রুটি প্রচারের গাণিতিক বিশ্লেষণ প্রদান করা
- পরীক্ষামূলক যাচাইকরণ যথেষ্ট: একাধিক কেস দ্বারা পদ্ধতির কার্যকারিতা যাচাই করা
- ব্যবহারিক মূল্য উচ্চ: গণনামূলক খরচ উল্লেখযোগ্যভাবে হ্রাস করা, পদ্ধতির ব্যবহারিকতা উন্নত করা
- সীমিত কেস: প্রধানত 2D RT অস্থিতিশীলতায় যাচাই করা, 3D এবং অন্যান্য প্রবাহ প্রকারের যাচাইকরণের অভাব
- তাত্ত্বিক সীমাবদ্ধতা: ত্রুটি বিশ্লেষণ এক-মাত্রিক সরলীকৃত পরিস্থিতির উপর ভিত্তি করে, প্রকৃত বহু-মাত্রিক পরিস্থিতি আরও জটিল
- নির্ভরতা সমস্যা: উন্নতির প্রভাব নির্দিষ্ট সংখ্যাগত বাস্তবায়ন এবং সমান্তরাল কৌশলের উপর নির্ভর করতে পারে
- একাডেমিক অবদান: MFM পদ্ধতির উন্নতির জন্য গুরুত্বপূর্ণ প্রযুক্তিগত পথ প্রদান করা
- ব্যবহারিক মূল্য: এডি ডিফিউসিভিটি কোএফিশিয়েন্টের মোমেন্টস গণনার দক্ষতা উল্লেখযোগ্যভাবে উন্নত করা
- পদ্ধতি প্রচার: ডিকম্পোজিশন MFM পদ্ধতির ভাল প্রচার সম্ভাবনা
- টার্বুলেন্ট মডেলিং: নির্ভুল এডি ডিফিউসিভিটি প্রয়োজনীয় RANS মডেলিং
- মিশ্রণ সমস্যা: Rayleigh-Taylor, Richtmyer-Meshkov ইত্যাদি অস্থিতিশীলতা গবেষণা
- অ-স্থানীয় প্রভাব গবেষণা: অ-স্থানীয় এবং অ্যানিসোট্রপিক প্রভাব পরিমাপ করার প্রয়োজনীয় প্রবাহ সমস্যা
এই পেপারটি 37টি গুরুত্বপূর্ণ সংদর্ভ উদ্ধৃত করে, যা টার্বুলেন্ট তত্ত্ব, সংখ্যাগত পদ্ধতি এবং MFM সম্পর্কিত গবেষণা অন্তর্ভুক্ত করে, গবেষণার জন্য একটি দৃঢ় তাত্ত্বিক ভিত্তি প্রদান করে।
সামগ্রিক মূল্যায়ন: এটি একটি উচ্চ-মানের গণনামূলক তরল গতিবিদ্যা পেপার যা ব্যবহারিক গণনামূলক সমস্যার জন্য কার্যকর সমাধান প্রস্তাব করে। পদ্ধতি সহজ এবং ব্যবহারিক, তাত্ত্বিক বিশ্লেষণ গভীর, পরীক্ষামূলক যাচাইকরণ যথেষ্ট, গুরুত্বপূর্ণ একাডেমিক মূল্য এবং ব্যবহারিক তাৎপর্য রয়েছে।