Bagaimana cara menggabungkan syarat batas dengan metode Galerkin?

21

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:

L[u(x,y)]=0on(x,y)Ω,S[u]=0on(x,y)Ω

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:LΩR2ΩSu(x,y)

u(x,y)i=1Naigi(x,y)

di mana adalah sekumpulan fungsi yang akan kita gunakan untuk memperkirakan . Mengganti dalam BVP:giu

iaiL[gi(x,y)]=R(a1,...,aN,x,y)

Karena perkiraan kami tidak tepat, sisa tidak tepat nol. Dalam metode Galerkin-Ritz-Raleigh kami meminimalkan sehubungan dengan set fungsi perkiraan dengan memerlukan . KarenanyaRRR,gi=0

R,gi=j=1NajL[gj],gi=0

Oleh karena itu, untuk menemukan koefisien , kita harus menyelesaikan persamaan matriks:ai

(L[g1],g1L[gN],g1L[g1],gNL[gN],gN)(a1aN)=0

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.S[u]

becko
sumber
1
becko, selamat datang di 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.
Geoff Oxberry
Namun, umumnya dianggap perilaku kasar untuk lintas pos. Jika Anda melihat daftar situs beta di Area 51, banyak dari mereka masih dalam versi beta publik setelah satu tahun. Kami masih akan berada di sana untuk sementara waktu (setidaknya pada skala waktu yang lebih lama daripada yang dibutuhkan sebagian besar pertanyaan di situs ini untuk dijawab). Selain itu, kecuali jika mathpengguna yang menjawab pertanyaan Anda juga adalah scicomppengguna, mereka tidak akan mendapatkan kredit atau atribusi yang sesuai untuk jawaban mereka di sini scicompjika Anda menyalin-tempel dari math, dan sebaliknya.
Geoff Oxberry
1
Apakah Anda yakin bahwa juga operator urutan kedua? Secara umum, itu bukan masalah yang diajukan dengan baik. Misalnya, jika , maka Anda meminta di yang memiliki solusi sangat unik (mis. Solusi PDE pada domain yang lebih besar dengan kondisi batas apa pun juga merupakan solusi). Biasanya kami meminta untuk menjadi operator urutan pertama (mungkin nonlinier). S = L L u = 0 ˉ Ω SSS=LLu=0Ω¯S
Jed Brown
2
Bahkan jika , Anda masih mencari solusi yang tidak unik. Pertimbangkan apakah adalah operator laplace dan adalah operator linear urutan kedua lainnya. Maka setiap sedemikian rupa sehingga untuk beberapa vektor konstan dapat ditambahkan ke solusi untuk membuat solusi lain. L S u u = k kSLLSuu=kk
Dan
1
@ GeoffOxberry Demi ketenangan pikiran Anda, pertanyaan duplikat di mathtelah dihapus. Jelas Anda benar tentang menyimpan pertanyaan di sini. Saya telah menerima tanggapan yang sangat membantu.
becko

Jawaban:

15

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+uNu0 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+τbpBbpττ

  • 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.

Allan P. Engsig-Karup
sumber
1
Saya kira adalah λ , bukan? τλ
shuhalo
Ya. dikoreksi.
Allan P. Engsig-Karup
1
Haruskah ia membaca "Metode Galerkin mengharuskan Anda memilih satu set fungsi dasar yang memenuhi BC masalah"?
knl
@knl: Saya pikir juga begitu, sedangkan kalimat lainnya tidak masuk akal. Saya akan melakukan edit.
davidhigh
7

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 .AbAbb

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.pu¯u¯

Bernardo MR
sumber
6

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)=0uf(x,y)δΩ

  • Anda dapat memperhalus elemen Anda sehingga sepenuhnya berada di batas berbagai elemenδΩ

sebenarnya sangat sederhana. Persamaan Anda:

perlu diganti dengan (

(L[g1],g1L[gN],g1L[g1],gNL[gN],gN)(a1aN)=0
(L[g1],g1L[gN],g1L[g1],gNL[gN],gN)(a1aN)=b

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 ΣbuδΩL[gj],gigjgiasebagai persamaan matriks, dan nilai-nilai dari elemenbharus pop keluar yang tepat sebagai produk batinLoperasi secara interior Anda dengan unsur-unsur Anda dasar batas.

R,gi=j=1NajL[gj],gi=0
bL
Dan
sumber
b
a2u/xy+b2u/x2=0)a,b
becko
LS
1
Saya mengedit pertanyaan sedikit untuk membuatnya lebih jelas, saya pikir. Saya tidak ingin memposting masalah yang sebenarnya saya coba selesaikan karena saya ingin menjaga pertanyaannya sedapat mungkin. Saya pikir saya bisa memahami metode itu dengan lebih baik.
becko
@becko: Kami mungkin ingin memindahkan ini untuk mengobrol , karena sudah agak lama.
Dan
2

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 .:

Lu=f
B(x)v(x)g(x)
u(x)v(x)+B(x)g(x)f(x)LB(x)
Lv=g
v(x)

B(x)

Berikutnya adalah penjelasan saya sendiri:

  • xu(x,y)|x=x0=1.

    B(x)

    xu(x,y)|x=x0=0.
  • u(x,y) = ijaijϕi(x)φj(y)
    xu(x,y) = ijaijϕi(x)φj(y)
    x=x0
  • ϕi(x)ϕi(x)|x=x0=0i

  • =1ϕi(x)|x=x0=1

    xu(x,y)|x=x0 = ijaijφj(y)
    1yaij

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.

Davidhigh
sumber