Apa gagasan umum metode Nitsche dalam analisis numerik?

17

Saya tahu bahwa metode Nitsche adalah metode yang sangat menarik karena memungkinkan untuk memperhitungkan kondisi batas tipe Dirichlet atau kontak dengan kondisi batas gesekan dengan cara yang lemah tanpa menggunakan pengganda Lagrange. Dan keuntungannya, yaitu mengubah kondisi batas Dirichlet menjadi istilah lemah seperti halnya kondisi batas Neumann, dibayar oleh fakta bahwa implementasi tergantung pada model.

Namun, sepertinya terlalu umum untukku. Bisakah Anda memberi saya ide yang lebih spesifik tentang metode ini? Contoh sederhana akan dihargai.

Anh-Thi DINH
sumber
Saya kira saya tidak mengerti pertanyaan Anda. Anda mengidentifikasi dengan benar mengapa metode ini ditemukan (untuk menangani kondisi Dirichlet dalam bentuk yang lemah). Apa yang Anda maksud dengan "Namun, sepertinya terlalu umum untuk saya. Bisakah Anda memberi saya gagasan yang lebih spesifik tentang metode ini? Contoh sederhana itu mahal."?
Wolfgang Bangerth
@ WolfgangBangerth: Saya perlu contoh (sederhana) untuk ide ini. Ini sangat abstrak bagi saya.
Anh-Thi DINH
2
@Liver: Saya berasumsi maksud Anda "mahal" seperti dalam "sayang", "berharga", yaitu, "dihargai"? Saya telah mengambil kebebasan untuk mengubah kata; jika Anda tidak setuju, silakan memutar kembali hasil edit.
Christian Clason

Jawaban:

24

Metode Nitsche terkait dengan metode Galerkin diskontinyu (memang, seperti yang ditunjukkan Wolfgang, itu adalah pendahulu dari metode ini), dan dapat diturunkan dengan cara yang serupa. Mari kita pertimbangkan masalah paling sederhana, persamaan Poisson: Kami sekarang mencari formulasi variasional itu

(1){-Δkamu=fdi Ω,kamu=gdi Ω.
  1. puas dengan solusi (lemah) (yaitu, konsisten),kamuH1(Ω)
  2. simetris dalam kamu dan v ,
  3. mengakui solusi unik (yang berarti bahwa bentuk bilinear adalah koersif).

Kita mulai seperti biasa dengan mengambil bentuk kuat dari persamaan diferensial, mengalikannya dengan fungsi tes dan mengintegrasikan dengan bagian-bagian. Dimulai dengan sisi kanan, kita memperoleh di mana dalam persamaan terakhir kami telah menambahkan nol produktif pada batas. Mengatur ulang istilah untuk memisahkan bentuk linear dan bilinear sekarang memberikan persamaan variasional untuk bentuk bilinear simetris yang dipenuhi untuk solusi dari .( f , v ) = ( - Δ u , v )vH1(Ω) 0=u-guH1(Ω)(1)

(f,v)=(-Δkamu,v)=(kamu,v)-Ωνkamuvds=(kamu,v)-Ωνkamuvds-Ω(kamu-g)νvds
0=kamu-gkamuH1(Ω)(1)

Namun, bentuk bilinear bukan paksaan, karena Anda tidak dapat mengikatnya dari bawah untuk oleh (karena kami tidak memiliki syarat batas untuk arbitrary , kita tidak dapat menggunakan ketidaksetaraan Poincaré seperti biasa - ini berarti kita dapat membuat bagian dari norma menjadi besar secara sewenang-wenang tanpa mengubah bentuk bilinear). Jadi kita perlu menambahkan istilah lain (simetris) yang hilang untuk solusi yang benar: untuk beberapa cukup besar. Ini mengarah pada formulasi lemah (simetris, konsisten, koersif): Cari sedemikian rupa sehingga c v 2 H 1 v H 1 ( Ω ) L 2 η Ω ( u - g ) vkamu=vcvH12vH1(Ω)L.2η > 0 u H 1 ( Ω ) ( u , v ) - Ω ∂ ∂ ν u vηΩ(kamu-g)vdsη>0kamuH1(Ω)

(u,v)ΩνuvdsΩuνvds+ηΩuvds=Ωgνvds+ηΩgvds+Ωfvdxfor all vH1(Ω).

Mengambil alih-alih pendekatan terpisah menghasilkan perkiraan Galerkin biasa. Perhatikan bahwa karena tidak sesuai karena kondisi batas (kami mencari solusi diskrit dalam ruang yang lebih besar dari yang kami cari solusi kontinu), seseorang tidak dapat menyimpulkan posisi yang baik dari masalah diskrit dari masalah terus menerus. Nitsche sekarang menunjukkan bahwa jika dipilih sebagai untuk cukup besar, masalah diskrit sebenarnya stabil (sehubungan dengan norma yang bergantung pada mesh yang cocok).kamu,vH1(Ω)kamuh,vhVhH1(Ω)ηch-1c>0

(Ini bukan derivasi asli Nitsche, yang mendahului metode Galerkin terputus-putus dan dimulai dari masalah minimisasi yang setara. Bahkan, makalah aslinya tidak menyebutkan bentuk bilinear yang sesuai sama sekali, tetapi Anda dapat menemukannya di, misalnya, Freund dan Stenberg, Pada kondisi batas yang diberlakukan secara lemah untuk masalah tingkat kedua , Prosiding Elemen Hingga Kesembilan Int. Elemen Hingga di Cairan, Venesia 1995. M. Morandi Cecchi dkk., Eds. Hlm. 327-336 .)

Christian Clason
sumber
1
Kalimat pertama Anda tidak salah, tetapi secara historis tidak akurat: Gagasan Nitsche didahulukan dan mengilhami pengembangan metode Galerkin yang terputus-putus. Yang mengatakan, ini tidak mengambil dari jawaban yang sangat baik.
Wolfgang Bangerth
@WolfgangBangerth Anda tentu saja benar; tidak ada kausalitas yang tersirat, hanya korelasi. Tetapi penting untuk memberikan atribusi yang tepat, terutama kepada orang-orang yang jika tidak mendapat giliran kerja pendek. Saya akan mengedit untuk memperjelasnya.
Christian Clason
Pertanyaan: 1. Bisakah Anda menguraikan lebih lanjut tentang masalah koersivitas sebelum menambahkan syarat batas tambahan? 2. Apa yang dimaksud dengan "ketidaksesuaian" di sini? 3. Saya pikir saya membaca bahwa stabilitas adalah hasil otomatis dari koersivitas bentuk bilinear ..? Meskipun penjelasan ini cukup bagus (satu-satunya penjelasan yang saya dapat temukan sebenarnya), adakah yang bisa menautkan ke penjelasan keseluruhan lain dari metode (dan / atau turunannya) hanya untuk perbandingan? Bahkan jika saya dapat menemukan kertas aslinya, tidak yakin itu akan banyak membantu. Makalah Freund dan Stenberg hanya memberikan sinopsis singkat dan spesifik pasangan
Malam
Ketidaksesuaian: ruang solusi diskrit bukan subruang dari ruang solusi kontinu - karena kondisi batas Dirichlet ditegakkan hanya dalam arti yang lemah. Berikut adalah tautan yang berpotensi bermanfaat . VhHg1(Ω)
GoHokies
@Malam Saya sudah mengedit jawaban untuk menjawab poin Anda (kecuali yang ada di paragraf kedua, jelas).
Christian Clason