Saya telah membaca beberapa sumber di web tentang metode Galerkin untuk menyelesaikan PDE, tapi saya tidak jelas tentang sesuatu. Berikut ini adalah catatan saya sendiri tentang apa yang saya pahami.
Pertimbangkan masalah nilai batas (BVP) berikut:
di mana adalah operator diferensiasi linier ke-2, adalah domain dari BVP, adalah batas domain, dan adalah operator diferensial linier pertama. Expess sebagai aproximation dari form:
di mana adalah sekumpulan fungsi yang akan kita gunakan untuk memperkirakan . Mengganti dalam BVP:
Karena perkiraan kami tidak tepat, sisa tidak tepat nol. Dalam metode Galerkin-Ritz-Raleigh kami meminimalkan sehubungan dengan set fungsi perkiraan dengan memerlukan . Karenanya
Oleh karena itu, untuk menemukan koefisien , kita harus menyelesaikan persamaan matriks:
Pertanyaan saya adalah: Bagaimana cara saya memasukkan syarat batas ke dalam ini?
EDIT: Awalnya pertanyaan mengatakan bahwa adalah operator diferensial linear urutan 2. Saya mengubahnya menjadi operator diferensial linear urutan pertama.
scicomp
! Kebijakan kami pada lintas-posting mengikuti kebijakan dari situs-situs Stack Exchange lainnya . Diperbolehkan melakukan posting silang jika Anda menyesuaikan pertanyaan yang sama (kurang lebih) untuk audiens yang berbeda. Anda diperbolehkan mengajukan pertanyaan untuk dimigrasikan ke situs lain setelah beberapa waktu, jika Anda merasa bahwa pertanyaan Anda tidak dijawab dengan memuaskan (atau bahkan tidak sama sekali) di situs tempat awalnya diposting.math
pengguna yang menjawab pertanyaan Anda juga adalahscicomp
pengguna, mereka tidak akan mendapatkan kredit atau atribusi yang sesuai untuk jawaban mereka di siniscicomp
jika Anda menyalin-tempel darimath
, dan sebaliknya.math
telah dihapus. Jelas Anda benar tentang menyimpan pertanyaan di sini. Saya telah menerima tanggapan yang sangat membantu.Jawaban:
Jawaban cepat dan umum tanpa abstraksi matematis. Ada beberapa opsi untuk memaksakan kondisi batas, misalnya
Sebenarnya metode Galerkin mengharuskan Anda memilih satu set fungsi basis yang memenuhi BC masalah (misalnya melalui basis rekombinasi dan / atau pemisahan pendekatan denganu 0uh=u0+uN u0 bertanggung jawab untuk solusi inhomogenous dan jumlah parsial yang bergantung pada fungsi dasar yang memenuhi kondisi homogen)uN
Metode hukuman / Lagrange mengalikan di mana orang pada dasarnya menambahkan istilah hukuman yang memasukkan kondisi batas, misalnya A + mana B adalah matriks yang bertanggung jawab untuk kondisi batas diskrit dan b p bertanggung jawab untuk istilah tidak homogen. Dalam batas τ → ∞ ketentuannya diberlakukan dengan kuat dan sebaliknya diberlakukan secara lemah. Pilihan τ mempengaruhi pengkondisian sistem.τ∗B=b+τ∗bp B bp τ→∞ τ
Metode Tau di mana sejumlah persamaan dipertukarkan (modifikasi baris dalam sistem Galerkin) dengan versi terpisah dari kondisi batas yang kemudian ditegakkan secara eksplisit. Catatan: satu opsi juga untuk membuat sistem overdetermined dengan kondisi batas tambahan.
Sebelum diskritisasi (Metode Ritz) menulis ulang formulasi Galerkin melalui teorema divergensi Gauss untuk mengubah integral volume ke integral batas dan kemudian menggabungkan (tepat atau kurang-lebih) kondisi batas langsung dalam formulasi sebelum diskritisasi.
Akhirnya, dengan mengeksploitasi koneksi antara ekspansi nodal / modal juga dimungkinkan untuk mendapatkan metode Galerkin nodal di mana solusi untuk sistem adalah koefisien dari basis Lagrange daripada basis modal.
sumber
Satu kemungkinan adalah untuk merakit sistem matriks dan vektor sisi kanan b , dengan derajat kebebasan yang ditentukan sebagai tidak diketahui, seperti tingkat kebebasan lainnya. Kemudian, A dan b dimodifikasi dengan memusatkan baris dan kolom yang terkait dengan dof yang ditentukan, dan menempatkan satu ke entri diagonal yang sesuai, dan memodifikasi vektor rh dengan tepat b .A b A b b
Ketika Anda nol baris, masukkan satu ke dalam diagonal dan ubah rhs sehingga Anda menegakkan nilai yang ditentukan, sistem tidak lagi simetris. Itu sebabnya Anda nol kolom dan memodifikasi vektor rhs untuk memperhitungkan nilai yang ditentukan.b
Kemungkinan lain adalah menambahkan angka sangat besar (biasanya 1e10) ke diagonal dof yang ditentukan dan kemudian mengatur entri rhs ke p * ˉ u , di mana ˉ u adalah nilai yang ditentukan dari dof itu.p u¯ u¯
sumber
Masalah umum berurusan dengan kondisi batas dengan metode elemen hingga bisa sangat rumit. Tapi jika:
sedemikian rupa sehingga satu-satunya pengenaan S ( u ) = 0 dibuat pada bentuk u adalah sama dengan beberapa f ( x , y ) pada δ Ω .S(u) S(u)=0 u f(x,y) δΩ
Anda dapat memperhalus elemen Anda sehingga sepenuhnya berada di batas berbagai elemenδΩ
sebenarnya sangat sederhana. Persamaan Anda:
perlu diganti dengan ( ⟨
di mana vektor sisi kanan mewakili kondisi batas.b
Untuk menentukan , mengatur unsur-unsur dasar Anda yang menentukan nilai u di delta ohm apa pun nilai-nilai mereka harus memenuhi kondisi batas. Dalam ⟨ L [ g j ] , g i ⟩ , Anda harus mengecualikan mereka dari g j tetapi bukan g i (unsur-unsur yang bersesuaian bahwa untuk fungsi-fungsi ini telah ditentukan, sehingga mereka tidak harus dimasukkan dalam matriks persamaan). Kemudian, mengatur ⟨ R , g i ⟩ = N Σb u δΩ ⟨L[gj],gi⟩ gj gi a sebagai persamaan matriks, dan nilai-nilai dari elemenbharus pop keluar yang tepat sebagai produk batinLoperasi secara interior Anda dengan unsur-unsur Anda dasar batas.
sumber
Berikut adalah metode yang dikenal sebagai basis rekombinasi , yang belum disebutkan dalam utas ini. Saya mengutip dari buku JP Boyd, "Chebyshev and Fourier Spectral Methods", 2nd Ed., Bab 6.5 .:
Berikutnya adalah penjelasan saya sendiri:
Yang menyenangkan tentang keseluruhan pendekatan ini adalah bahwa ia bekerja pada tingkat yang relatif abstrak. Bahan-bahan yang diperlukan hanya linieritas dari operator BC dan ansatz dalam hal fungsi basis produk. Dengan demikian, ini juga berlaku untuk metode perkiraan.
sumber