Tekanan sebagai Pengganda Lagrange

12

Dalam persamaan Navier-Stokes yang tidak dapat dimampatkan, istilah tekanan sering disebut sebagai pengganda Lagrange yang menegakkan kondisi tidak tertahankan.

ρ(ut+(u)u)=p+μΔu+fu=0

Dalam hal apa ini benar? Apakah ada rumusan persamaan Navier-Stokes yang tidak dapat dikompres sebagai masalah optimasi yang tunduk pada kendala ketidakterkompresan? Jika demikian, apakah ada analog numerik di mana persamaan aliran fluida yang tidak dapat dipecahkan diselesaikan dalam kerangka kerja optimasi?

Ben
sumber

Jawaban:

18

Ini paling mudah dilihat dengan mempertimbangkan persamaan Stokes stasioner yang setara dengan masalah Jika Anda menuliskan Lagrangian dan kemudian kondisi optimalitas dari masalah optimasi ini, Anda akan menemukan bahwa memang tekanannya adalah pengganda Lagrange.

μΔu+p=fu=0
minuμ2u2(f,u)so thatu=0.

Kesetaraan di antara masalah-masalah ini tidak dieksploitasi dalam skema numerik apa pun (yang saya ketahui) tetapi ini merupakan alat penting dalam analisis karena ini menunjukkan bahwa persamaan Stokes pada dasarnya adalah persamaan Poisson pada subruang linear. Hal yang sama berlaku untuk persamaan Stokes yang tergantung waktu (yang sesuai dengan persamaan panas pada subruang) dan dapat diperluas ke persamaan Navier-Stokes.

Wolfgang Bangerth
sumber
Terima kasih atas jawaban yang bagus. Apakah Anda tahu jika formulasi ini dapat diperluas ke kasus tergantung waktu?
Ben
1
Ya, seperti yang saya katakan itu mengarah pada persamaan panas pada subruang fungsi divergensi bebas.
Wolfgang Bangerth
1
Maaf, saya seharusnya lebih jelas. Apakah ada cara untuk menyusun kembali persamaan Stokes yang bergantung waktu (atau Navier-Stokes) sebagai masalah optimisasi, mungkin dari fungsional yang terintegrasi dari waktu ke waktu?
Ben
1
Bukan sebagai masalah optimisasi - solusi persamaan panas tidak meminimalkan apa pun (meskipun itu adalah titik stasioner dari fungsi Lagrangian). Tetapi Anda dapat merumuskan persamaan Stokes sebagai berikut: Temukan sehingga untuk semua tunduk pada batasan yang . Perhatikan bahwa saya telah memilih ruang uji yang lebih kecil daripada ruang uji coba sehingga sisi kiri dan kanan persamaan variasional tidak akan sama. Perbedaannya adalah tekanan. uHdiv(ut,φ)+(u,φ)=(f,φ)φ{vHdiv:v=0}u=0
Wolfgang Bangerth