Dimulai dengan persamaan advection dalam bentuk konservasi.
di mana adalah kecepatan yang bergantung pada ruang, dan adalah konsentrasi spesies yang dilestarikan.u
Mendiskritkan fluks (di mana fluks , didefinisikan pada tepi sel di antara titik-titik mesh) memberi, u t = 1
Menggunakan urutan pertama melawan arah angin kami memperkirakan fluks sebagai,
Jika konstan maka ini akan mengurangi ke skema arah angin yang sudah dikenal yaitu, .u t = a
Pertanyaan saya adalah, bagaimana kita dapat memperlakukan koefisien tidak konstan dari persamaan advection? Kecepatan didefinisikan di pusat sel, jadi pendekatan sederhana adalah sebagai berikut,
Ini adalah pendekatan yang saya sukai karena sangat mudah diterapkan.
Namun, kita juga bisa menggunakan (saya menebak) skema rata-rata untuk menentukan kecepatan di tepi sel,
Dalam buku LeVeque dia berkata,
Sejauh ini kita telah mengasumsikan bahwa kecepatan variabel ditentukan oleh nilai konstan dalam sel grid j-th. Dalam beberapa kasus, lebih alami untuk mengasumsikan bahwa kecepatan ditentukan pada setiap antarmuka sel.a j a j - 1
Tapi dia tidak benar-benar menguraikan terlalu banyak setelah itu. Apa itu pendekatan umum?
Saya memecahkan masalah konservasi (saya menggunakan persamaan advection sebagai persamaan kontinuitas) jadi saya ingin memastikan bahwa setelah menerapkan diskritisasi, properti konservasi akan dilestarikan. Saya ingin menghindari kejutan tersembunyi mengenai koefisien variabel ini! Adakah yang punya komentar dan panduan umum?
Pembaruan Ada dua jawaban yang sangat bagus di bawah ini dan saya hanya bisa memilih satu :(
sumber
Yang lain mengatakan itu semua, tapi saya hanya ingin menambahkan poin yang sederhana, namun terkadang halus. Diskritisasi melawan angin Anda tetap konservatif selama Anda menggunakan interpolasi konsisten dari pada batas sel.a(x)
Yang saya maksud dengan konsisten adalah bahwa satu - satunya syarat yang perlu dipenuhi oleh interpolasi adalah
Dengan kata lain, selama metode interpolasi Anda terus menerus melintasi batas sel, diskritisasi Anda dijamin akan tetap konservatif.
Ini mungkin tidak tampak seperti masalah besar di sini dalam 1D (dan seharusnya tidak) tetapi dapat menyebabkan masalah pada antarmuka kasar-halus pada kisi-kisi AMR multi-level.
sumber
Anda dapat menggunakan semua jenis interpolasi untuk menentukan , dan metode ini akan tetap konservatif.a(xj−12)
Untuk melihat mengapa demikian, pertimbangkan bahwa definisi analitik konservatif adalah itu
di mana adalah domain masalah. Ini mengatakan bahwa perubahan dalam jumlah yang dilestarikan sama dengan fluks pada batas, dan segera mengikuti dari mengintegrasikan hukum konservasi.D
Jika diskritisasi kami berbentuk
di mana adalah titik kisi kami, , dan , , maka pernyataan konservasi setara diskrit adalahx1,…,xn D=[c,d] c=x12 d=xn+12
dan ini dapat dengan mudah diamati untuk dipegang dengan memperluas jumlah di sisi kiri. Perhatikan bahwa, dalam kasus Anda dengan upwinding, dan , meskipun harus dicatat bahwa skema ini hanya skema melawan angin yang tepat jika selalu positif.uj+1uj−12=uj−1 a(x)uuj+12=uj a(x)u
Untuk metode tingkat tinggi, asalkan halus, seseorang dapat dengan mudah memasukkan polinomial ke titik , dan mengevaluasi polinomial di .a ( x j - r ) , … , a ( x j + s ) a ( x j - 1a(x) a(xj−r),…,a(xj+s) a(xj−12)
sumber