Ketika beralih dari bentuk kuat PDE ke bentuk FEM, tampaknya orang harus selalu melakukan ini dengan terlebih dahulu menyatakan bentuk variasional. Untuk melakukan ini, Anda mengalikan formulir kuat dengan elemen di beberapa (Sobolev) ruang dan mengintegrasikan di wilayah Anda. Ini bisa saya terima. Yang tidak saya mengerti adalah mengapa seseorang juga harus menggunakan formula Green (satu atau beberapa kali).
Saya sebagian besar telah bekerja dengan persamaan Poisson, jadi jika kita mengambil itu (dengan kondisi batas Dirichlet homogen) sebagai contoh, yaitu
maka diklaim bahwa cara yang tepat untuk membentuk bentuk variasional adalah
Tapi apa yang menghentikan saya dari menggunakan ekspresi pada baris pertama, bukankah itu juga bentuk variasi yang dapat digunakan untuk mendapatkan formulir FEM? Bukankah itu sesuai dengan bilinear dan bentuk linear dan ? Apakah masalah di sini bahwa jika saya menggunakan fungsi basis linear (fungsi bentuk) maka saya akan berada dalam kesulitan karena matriks kekakuan saya akan menjadi matriks nol (tidak dapat dibalik)? Tetapi bagaimana jika saya menggunakan fungsi bentuk non-linear? Apakah saya masih harus menggunakan formula Green? Jika saya tidak harus: apakah itu disarankan? Jika tidak, apakah saya kemudian memiliki formulasi variasi-tetapi-tidak-lemah?l ( v ) = ( f , v )
Sekarang, katakanlah saya memiliki PDE dengan turunan orde tinggi, apakah itu berarti ada banyak kemungkinan bentuk variasi, tergantung pada bagaimana saya menggunakan rumus Green? Dan mereka semua mengarah pada perkiraan FEM (berbeda)?
sumber
Jawaban:
Jawaban singkat:
Tidak, Anda tidak perlu melakukan integrasi untuk FEM tertentu. Tetapi dalam kasus Anda, Anda harus melakukan itu.
Jawaban panjang:
Katakanlah adalah solusi elemen hingga. Jika Anda memilih sedikit polinomial linier sebagai basis Anda, maka mengambil di atasnya akan memberi Anda distribusi pesanan 1 (berpikir mengambil turunan pada fungsi langkah Heaviside), dan integrasi mengalikan dengan hanya akan masuk akal bila Anda menganggapnya sebagai pasangan dualitas daripada produk dalam . Anda tidak akan mendapatkan matriks nol, teorema representasi Riesz mengatakan bahwa ada elemen dalam dapat mengkarakterisasi pasangan dualitas dengan produk dalam di : Δ - Δ u h ∈ H - 1 v L 2 φ - Δ u h ∈ H 1 0 H 1 ⟨ - Δ u h , v ⟩ H - 1 , H 1 0 = ∫ ohm ∇ φ - Δ u h ⋅ ∇ v ⏟ produk dalam H 1 . u h Tuh Δ −Δuh∈H−1 v L2 φ−Δuh∈H10 H1
Integrasi oleh bagian dapat ditelusuri kembali ke teori Sobolev untuk elips pde menggunakan fungsi halus, di mana -spaces semua penutupan fungsi halus di bawah jenis norma integral. Lalu orang mengatakan apa keteraturan minimum di sini bahwa kita dapat melakukan produk dalam. Juga perlu diingat bahwa solusi lemah reguler dalam kondisi tertentu adalah solusi kuat (keteraturan elips). Tapi polinomial linear kontinu piecewise bukanlah , dari sudut pandang ini, tidak masuk akal untuk mengambil produk dalam menggunakan baik. W k , p H 1 H 2 H 2 Δ u hWk,p Wk,p H1 H2 H2 Δuh
Untuk FEM tertentu, Anda tidak perlu melakukan integrasi dengan bagian-bagian. Misalnya, elemen hingga kuadrat-terkecil. Tulis pde urutan kedua sebagai sistem urutan pertama: Maka Anda ingin meminimalkan fungsional kuadrat-terkecil: memiliki semangat yang sama dengan fungsional Ritz-Galerkin, formulasi elemen hingga dari meminimalkan fungsional di atas dalam suatu ruang elemen hingga tidak memerlukan integrasi oleh bagian-bagian.
sumber
Tidak ada yang menghentikan Anda melakukan hal itu secara teknis, tetapi ketika Anda mengintegrasikan bagian-bagiannya Anda mendapatkan lebih banyak fleksibilitas dengan ruang solusi karena mereka tidak perlu memiliki keteraturan (diperlukan untuk formulasi non IBP). Elemen-elemen linier yang Anda sarankan pada umumnya telah menegakkan kesinambungan antara elemen-elemen, dan karenanya tidak mungkin dalam . Formulasi IBP selanjutnya adalah simetris, yang memiliki beberapa kelebihannya sendiri.H 2H2 H2
sumber
Jawaban yang sangat baik sudah ada di halaman ini, tetapi masih ada (kecil) poin yang hilang.
OP bertanya:
Mengintegrasikan oleh bagian-bagian (dengan cara yang benar ) penting ketika Anda memiliki tipe kondisi batas Neumann. Sebenarnya oleh ibp Anda memperhitungkan Neumann bc dalam formulasi variasional Anda. Bentuk Neumann bc tergantung pada bagaimana Anda mengintegrasikan bagian, lih. jawaban ini pada integrasi oleh bagian-bagian dalam elastisitas linier. Jadi, bahkan untuk PDE eliptik orde dua, integrasi oleh bagian harus dilakukan dengan cara tertentu, untuk memulihkan formulasi variasi yang valid untuk Neumann atau kondisi batas campuran. (Dan ini tentu saja terlepas dari kenyataan bahwa Anda diskritkan oleh FEM).
Dalam fisika matematika, di mana Neumann bc memiliki makna yang terdefinisi dengan baik (fluks panas, stres ...), integrasi oleh bagian-bagian penting untuk menjaga interpretasi hasil yang benar. Bahkan untuk kondisi Dirichlet dan FEM yang homogen ini benar, karena jika kita menggunakan metode pengali Lagrange untuk memaksakan bc, pengali menjadi kuantitas fisik, seperti fluks atau gaya pekat.
sumber