Mengikuti dari pertanyaan saya sebelumnya, saya mencoba menerapkan kondisi batas pada jaring volume terbatas yang tidak seragam ini,
Saya ingin menerapkan kondisi batas tipe Robin ke lhs domain ( , sehingga,
di mana adalah nilai batas; adalah koefisien yang didefinisikan pada batas, adveksi dan difusi masing-masing; , adalah turunan dari dievaluasi pada batas dan adalah variabel yang sedang kita pecahkan. a , d u x = ∂ u uu
Kemungkinan pendekatan
Saya dapat memikirkan dua cara untuk menerapkan kondisi batas ini pada jala volume terbatas di atas:
Pendekatan sel hantu.
Tulis sebagai perbedaan hingga termasuk sel hantu. σ L = d u 1 - u 0
A. Kemudian gunakan interpolasi linier dengan poin dan untuk menemukan nilai perantara, .x 1 u ( x L )
B. Atau temukan dengan rata-rata di atas sel,u ( x L ) = 1
Dalam kedua kasus, ketergantungan pada sel hantu dapat dihilangkan dengan cara biasa (melalui substitusi ke dalam persamaan volume hingga).
Pendekatan ekstrapolasi.
Cocokkan fungsi linier (atau kuadratik) dengan dengan menggunakan nilai pada titik ( ). Ini akan memberikan nilai pada u ( x L ) . Fungsi linear (atau kuadrat) kemudian dapat dibedakan untuk menemukan ekspresi untuk nilai turunannya, u x ( x L ) , pada batas. Pendekatan ini tidak menggunakan sel hantu.x 1 , x 2 x 3
Pertanyaan
- Pendekatan mana dari ketiganya, (1A, 1B atau 2) yang "standar" atau Anda sarankan?
- Pendekatan mana yang memperkenalkan kesalahan terkecil atau paling stabil?
- Saya pikir saya dapat menerapkan pendekatan sel hantu sendiri, namun, bagaimana pendekatan ekstrapolasi dapat diterapkan, apakah pendekatan ini memiliki nama?
- Apakah ada perbedaan stabilitas antara pemasangan fungsi linier atau persamaan kuadrat?
Persamaan khusus
Saya ingin menerapkan batas ini pada persamaan advection-difusion (dalam bentuk konservasi) dengan istilah sumber non-linear,
Mendiskritkan persamaan ini pada mesh di atas menggunakan metode memberikan,
Namun untuk titik batas ( ) saya lebih suka menggunakan skema sepenuhnya implisit ( θ = 1 ) untuk mengurangi kompleksitas,
Perhatikan titik hantu , ini akan dihapus dengan menerapkan syarat batas.
Koefisien memiliki definisi,
Semua variabel " " didefinisikan seperti pada diagram di atas. Akhirnya, Δ t yang merupakan langkah waktu ( NB ini adalah disederhanakan kasus dengan konstan suatu dan d koefisien, dalam praktek " r " koefisien sedikit lebih rumit untuk alasan ini).
sumber
Jawaban:
Ini bukan komentar umum tentang FVM daripada jawaban atas pertanyaan konkret. Dan pesannya adalah bahwa seharusnya tidak ada kebutuhan untuk diskritisasi adhoc seperti kondisi batas.
Tidak seperti dalam metode FE- atau FD, di mana titik awalnya adalah solusi diskrit untuk solusi, pendekatan FVM membuat solusi tidak tersentuh (pada awalnya) tetapi rata-rata pada segmentasi domain. Diskritisasi solusi hanya berlaku ketika sistem persamaan keseimbangan yang diperoleh diubah menjadi sistem persamaan aljabar dengan cara memperkirakan fluks di seluruh antarmuka.
Dalam hal ini, mengingat kondisi batas, saya menyarankan untuk tetap menggunakan bentuk solusi yang berkelanjutan selama mungkin dan untuk memperkenalkan pendekatan diskrit hanya di bagian paling akhir.
But now, to turn this into an algebraic equation, one typically assumes that on cellCi the function u is constant in space, i.e. u(t,x)|Ci=ui(t) . Thus, having associated u(xi)≈ui , one can express ux|hi at the cell boarders via the difference quotient in ui and ui+1 . To express u at the cell boarders one can use interpolation (i.e. central differences or upwind schemes).
What to do at the boundary? In the example, it is all about approximating(−au+dux)|x=0 , no matter what has been done to u so far.
Givenu|x=0=gD one can introduce a ghost cell and the condition that an interpolant between u0 and u1 is equal to gD at the boarder.
Givenux|x=0=gN one can introduce a ghost cell and the condition that an approximation to the derivative between u0 and u1 matches gN at the boarder
If the flux itself is prescribed:(−au+dux)|x=0=gR , there is no need for a discretization.
However, I am not sure, what to do in the case that there are Robin type bc's that do not match the flux directly. This, will need some regularization because of the discontinuity of the advection and diffusion parameters.
===> Some personal thoughts on FVM <===
sumber