আপনি গণ সংরক্ষণের সমীকরণটি দেখছেন:
dmdt=0
প্রতি ইউনিট ভলিউম ভর বিবর্তনের বিবেচনা করার সময়, এটি ফ্লাক্স আকারে ঘনত্ব অ্যাডভেকশন সমীকরণে ফোটে:
∂ρ∂t=−∇⋅(ρu)
এই সম্পর্কে ভাল জিনিস এটা একটি অবাধ স্কালে মাঠের মাত্র advection সমীকরণ যে (আমাদের ক্ষেত্রে, এই ঘনত্ব হতে হবে হয় পর্যাপ্ত সময় এবং স্থান ডিফারেন্সিং স্কিম, এবং প্রাথমিক এবং প্রদত্ত) এবং (অপেক্ষাকৃত) হল সহজ সমাধান, সীমানা শর্ত.ρ
সীমাবদ্ধ বিবিধ স্কিম ডিজাইন করার সময়, আমরা রূপান্তর, স্থিতিশীলতা এবং নির্ভুলতার বিষয়ে চিন্তা করি। একটি প্রকল্প যদি সমকেন্দ্রি হয় যখনΔt→0। স্কিম নিশ্চিত স্থায়িত্ব পরিমাণ যেএকটিসসীম যখন রয়ে যায়টন→∞। স্কিমের আনুষ্ঠানিক নির্ভুলতা জানায় যে আংশিক ডেরাইভেটিভ মিথ্যাচারের টেলর সম্প্রসারণ সিরিজে কাটা ত্রুটি। বিচ্ছিন্ন স্কিমের এই মৌলিক বৈশিষ্ট্যগুলির বিষয়ে আরও তথ্যের জন্য একটি সিএফডি পাঠ্যপুস্তকটি সন্ধান করুন।ΔAΔt→∂A∂tΔt→0At→∞
এখন, সহজ পদ্ধিতিটি হ'ল সরাসরি প্রবাহে পৃথক পৃথক 1 ম আদেশে যাওয়া। এই স্কিমটি ইতিবাচক-সুনির্দিষ্ট, রক্ষণশীল এবং গণনামূলকভাবে দক্ষ। প্রথম দুটি বৈশিষ্ট্য বিশেষত গুরুত্বপূর্ণ যখন আমরা একটি পরিমাণের বিবর্তনকে মডেল করি যা সর্বদা ধনাত্মক (যেমন ভর বা ঘনত্ব)।
সরলতার জন্য, আসুন 1-D কেসটি দেখুন:
∂ρ∂t=−∂(ρu)∂x
ফ্লাক্স সংজ্ঞায়িত করা এখন সুবিধাজনক , যাতে:Φ=ρu
∂(ρu)∂x=∂Φ∂x≈ΔΦΔx≈Φi+1/2−Φi−1/2Δx
আমরা যা অনুকরণ করছি তার একটি পরিকল্পনাকারী এখানে:
u u
| --> --> |
| rho | rho | rho |
x-----o-----x-----o-----x-----o-----x
i-1 i-1/2 i i+1/2 i+1
আমরা বিবর্তন মূল্যায়নের হয় সেল এ আমি । নেট লাভ বা লস কি, আসে পার্থক্য থেকে আসে Φ আমি - 1 / 2 এবং কি খুঁজে যায়, Φ আমি + + 1 / 2ρiΦi−1/2Φi+1/2। এখানেই আমরা পলের উত্তর থেকে সরে যেতে শুরু করি। সত্য রক্ষণশীল উজান প্রবাহের ভিন্নতায়, কক্ষের কেন্দ্রে পরিমাণটি তার গতির দিক দিয়ে, তার ঘরের প্রান্তে বেগ দ্বারা বহন করে চলেছে। অন্য কথায়, আপনি যদি ধারণা করেন যে আপনি অনুমোদিত পরিমাণ এবং আপনি কোষ কেন্দ্রে বসে আছেন, আপনাকে ঘরের প্রান্তে বেগ দিয়ে আপনার সামনে কোষে নিয়ে যাওয়া হবে। উভয় ঘরের প্রান্তে ঘনত্ব এবং বেগের পণ্য হিসাবে সেল প্রান্তে প্রবাহগুলি মূল্যায়ন সঠিক নয় এবং এটি অনুমোদিত পরিমাণ সংরক্ষণ করে না।
আগত এবং বহির্গামী ফ্লাক্সগুলি হিসাবে মূল্যায়ন করা হয়:
Φi+1/2=ui+1/2+|ui+1/2|2ρi+ui+1/2−|ui+1/2|2ρi+1
Φi−1/2=ui−1/2+|ui−1/2|2ρi−1+ui−1/2−|ui−1/2|2ρi
The above treatment of flux differencing ensures upstream-definiteness. In other words, it adjusts the differencing direction according to the sign of velocity.
The Courant-Friedrichs-Lewy (CFL) stability criterion, when doing time differencing with simple first order, forward Euler differencing is given as:
μ=uΔtΔx≤1
Note that in 2 dimensions, the CFL stability criterion is more strict:
μ=cΔtΔx≤12–√
where c is velocity magnitude, u2+v2−−−−−−√.
Some things to consider. This scheme may or may not be appropriate for your application depending on what kind of process you are simulating. This scheme is highly diffusive, and is appropriate for very smooth flows without sharp gradients. It is also more diffusive for shorter time steps. In the 1-D case, you will obtain an almost exact solution if the gradients are very small, and if μ=1. In the 2-D case, this is not possible, and diffusion is anisotropic.
If your physical system considers shock waves or high gradients of other sort, you should look into upstream differencing of higher order (e.g. 3rd or 5th order). Also, it may be worthwhile looking into the Flux Corrected Transport family of schemes (Zalesak, 1979, JCP); anti-diffusion correction for the above scheme by Smolarkiewicz (1984, JCP); MPDATA family of schemes by Smolarkiewicz (1998, JCP).
For time differencing, 1st order forward Euler differencing may be satisfactory for your needs. Otherwise, look into higher-order methods such as Runge-Kutta (iterative), or Adams-Bashforth and Adams-Moulton (multi-level).
It would be worthwhile looking into some CFD graduate-level textbook for a summary of above mentioned schemes and many more.