Penggunaan optimal Strang splitting (untuk persamaan difusi reaksi)

Saya membuat pengamatan aneh sambil menghitung solusi untuk persamaan difusi reaksi 1D sederhana:

$\frac{\partial}{\partial t}a=\frac{\partial^2}{\partial x^2}a-ab$

$\frac{\partial}{\partial t}b=-ab$

$\frac{\partial}{\partial t}c = a$

Nilai awal adalah konstan ( ), dan aku hanya tertarik terpisahkan atas dari ke ( ). Tujuan dari dan persamaan hanya untuk mengevaluasi integral ini. $b$ $b(0,x)=b_0$ $a$ $0$ $1$ $\int_0^1a(t,x)dt$ $c$ $\frac{\partial}{\partial t}c = a$

Saya menggunakan skema pemisahan Strang untuk penggandaan antara difusi dan reaksi (reaksi setengah langkah, kemudian difusi langkah penuh, dan kemudian lagi setengah langkah reaksi), skema Crank Nicholson untuk difusi, dan solusi analitis untuk reaksi ( termasuk persamaan ). $\frac{\partial}{\partial t}c = a$

Karena satu langkah dari solusi analitis lebih dari faktor 3 lebih lambat dari satu langkah skema Crank Nicholson, saya mencoba membuat lebih dari satu langkah Crank Nicholson untuk setiap langkah reaksi. Saya berharap untuk bertahan dengan langkah yang lebih sedikit dari skema pemisahan Strang, sehingga saya akan lebih cepat secara keseluruhan.

Namun, efek sebaliknya dapat diamati, yaitu bahwa lebih banyak langkah untuk skema pemisahan Strang diperlukan jika lebih dari satu langkah Crank Nicholson digunakan. (Saya hanya peduli dengan akurasi integral atas , yang tampaknya konvergen lebih cepat dari itu sendiri.) Setelah bertanya-tanya untuk beberapa waktu, saya melihat bahwa efek yang sama juga terjadi untuk , dan saya bahkan mengerti mengapa untuk kasus ini. Intinya adalah bahwa jika saya membuat tepat satu langkah Crank Nicholson, maka skema keseluruhan berubah menjadi aturan trapesium (jika ). $a$ $a$ $b(t,x)=b_0=0$ $b(t)=0$

Jadi jika saya akan memperlakukan sebagai bagian dari langkah difusi, meningkatkan jumlah langkah Crank Nicholson (mungkin) tidak akan menyebabkan berkurangnya akurasi keseluruhan (seperti yang diamati). Tapi itu tampaknya mengalahkan tujuan menggunakan solusi analitis untuk bagian reaksi (non-linear dan berpotensi sangat kaku) dari sistem. $\frac{\partial}{\partial t}c = a$

Jadi inilah pertanyaan saya: Apakah ada cara yang lebih baik untuk memperlakukan dalam konteks pemisahan Strang, daripada memperlakukannya sebagai bagian dari langkah reaksi, atau untuk memperlakukan itu sebagai bagian dari langkah difusi. Saya ingin menghindari "dipaksa" untuk menggunakan tepat satu langkah Crank Nicholson untuk difusi. (Misalnya dalam 3D, saya lebih suka menyelesaikan difusi secara analitik dengan FFT daripada menggunakan Crank Nicholson. Tentu saja saya juga bisa menggabungkan FFT dengan Crank Nicholson, jadi ini bukan masalah besar.) $\frac{\partial}{\partial t}c = a$

parabolic-pde operator-splitting

— Thomas Klimpel
sumber

Di people.maths.ox.ac.uk/dellar/OperatorLB.html , efek yang serupa tampaknya dijelaskan. Kesimpulannya adalah sangat penting untuk menggunakan Crank Nicholson alih-alih solusi yang tepat. Jadi mungkin jawaban untuk pertanyaan saya adalah tidak.

— Thomas Klimpel

Tampaknya ada yang salah dengan persamaan Anda. tidak muncul dalam dua yang pertama membuat kopling satu arah dan artinya Anda dapat menghitung pada sebagai langkah pasca-pemrosesan.

c

$c$

c

$c$

t

$t$

— Bill Barth

@ BillBarth Saya mengubah pertanyaan untuk memperjelas peran . Jadi hanyalah alat untuk menghitung . Tolong beri tahu saya jika Anda memiliki saran cara menghitung integral ini lebih akurat (daripada apa yang saya dapatkan dari kombinasi pemisahan Strang dan Crank Nicholson yang dijelaskan di atas), berpotensi menggunakan langkah pasca-pemrosesan.

c

$c$

c

$c$

\int_{0}^{1} a (t, x) d t

$\int_0^1a(t,x)dt$

— Thomas Klimpel

Ini sekarang sudah lama berlalu, tetapi apakah Anda menyadari bahwa sistem persamaan ini dapat ditulis sebagai PDE parabola dalam dengan istilah reaksi eksponensial? Saya kira saya ingin tahu apakah Anda benar-benar ingin menyelesaikan sistem 3-variabel ini daripada yang disederhanakan.

c

$c$

— Bill Barth

@BillBarth Saya akan tertarik untuk mempelajari bagaimana sistem ini dapat ditulis sebagai PDE parabola dengan istilah reaksi eksponensial. Kecepatan solusi model ini adalah faktor pembatas selama kalibrasi model (yang bisa memakan waktu beberapa jam), meskipun demikian akurasi yang digunakan sehubungan dengan integrasi waktu cukup jauh dari konvergensi penuh.

— Thomas Klimpel

Saya akan menulis ini sebagai jawaban walaupun tidak langsung menjawab pertanyaan.

Memasukkan persamaan kedua dan persamaan ketiga ke yang pertama, dan memasukkan yang ketiga ke yang kedua, bersama-sama memberikan: Menyusun ulang keduanya memberi: Sekarang, kita dapat mengintegrasikan keduanya sekaligus dalam , pergi ke persamaan pertama:

\begin{aligned} \frac{\partial^{2} c}{\partial t^{2}} & = \frac{\partial^{2}}{\partial x^{2}} \frac{\partial c}{\partial t} + \frac{\partial b}{\partial t} \\ \frac{\partial b}{\partial t} & = - (\frac{\partial c}{\partial t}) b \end{aligned}

$\begin{align} \frac{\partial^2 c}{\partial t^2} &= \frac{\partial^2}{\partial x^2}\frac{\partial c}{\partial t} + \frac{\partial b}{\partial t} \\ \frac{\partial b}{\partial t} &= - \left(\frac{\partial c}{\partial t}\right) b \end{align}$

\begin{aligned} \frac{\partial}{\partial t} (\frac{\partial c}{\partial t} - \frac{\partial^{2} c}{\partial x^{2}} - b) & = 0 \\ \frac{1}{b} (\frac{\partial b}{\partial t}) & = - \frac{\partial c}{\partial t} \end{aligned}

$\begin{align} \frac{\partial}{\partial t}\left(\frac{\partial c}{\partial t} - \frac{\partial^2 c}{\partial x^2} - b \right) &= 0 \\ \frac{1}{b}\left(\frac{\partial b}{\partial t}\right) &= -\frac{\partial c}{\partial t} \end{align}$

t

$t$

\frac{\partial c}{\partial t} - \frac{\partial^{2} c}{\partial x^{2}} = b + A (x)

$\frac{\partial c}{\partial t} - \frac{\partial^2 c}{\partial x^2} = b + A(x)$ Dengan menggunakan persamaan ketiga, kita dapat mengekspresikan "konstan" integrasi sebagai . Persamaan kedua sedikit lebih rumit. Menulis ulang sedikit, kita memiliki: Ini mengarah ke solusi Atau eksponensial memberi: Dan, akhirnya, menghubungkan ini ke PDE untuk memberi

A (x) = a_{0} - \frac{\partial^{2} c_{0}}{\partial x^{2}} - b_{0}

$A(x)=a_0-\frac{\partial^2 c_0}{\partial x^2}-b_0$

\int_{0}^{t} \frac{1}{b (x, t^{'})} (\frac{\partial b (x, t^{'})}{\partial t^{'}}) d t^{'} = - \int_{0}^{t} \frac{\partial c (x, t^{'})}{\partial t^{'}} d t^{'}

$\int_0^t \frac{1}{b(x,t')}\left(\frac{\partial b(x,t')}{\partial t'}\right) dt' = -\int_0^t \frac{\partial c(x,t')}{\partial t'} dt'$

\ln b (x, t) - \ln b_{0} (x) = - c (x) + c_{0} (x)

$\ln b(x,t) - \ln b_0(x) = - c(x) + c_0(x)$

\ln \frac{b}{b_{0}} = - c + c_{0}

$\ln \frac{b}{b_0} = -c + c_0$

b = b_{0} e^{c_{0} - c}

$b = b_0 e^{c_0-c}$

c

$c$

\frac{\partial c}{\partial t} - \frac{\partial^{2} c}{\partial x^{2}} = b_{0} e^{c_{0} - c} + A (x)

$\frac{\partial c}{\partial t} - \frac{\partial^2 c}{\partial x^2} = b_0 e^{c_0-c}+A(x)$

Mengganti dengan , atau menggunakan kondisi awal , persamaan ini menyederhanakan menjadi Sekarang, Anda harus dapat menemukan literatur yang cukup tentang cara terbaik untuk menyelesaikan persamaan ini. Saya tidak tahu apakah Crank-Nicholson adalah pilihan yang baik untuk istilah eksponensial, tetapi tampaknya masuk akal. Beberapa perawatan mungkin harus diambil untuk menjamin bahwa mana-mana atau solusi dapat meledak dengan cepat. $c$ $c-c_0$ $c_0(x)=0$

\frac{\partial c}{\partial t} = \frac{\partial^{2} c}{\partial x^{2}} + a_{0} - (1 - e^{- c}) b_{0}

$\frac{\partial c}{\partial t} = \frac{\partial^2 c}{\partial x^2} + a_0 - (1-e^{-c}) b_0$

c > c_{0}

$c > c_0$

Saya hanya berjalan melalui derivasi ini dua kali, jadi mungkin ada satu atau dua kesalahan di dalamnya, tetapi itu terasa tepat bagi saya. Jika mana-mana, maka ini jelas merupakan solusi yang tepat, dan memiliki aroma masuk akal sebaliknya. $b_0 = 0$

Memecahkan ini sampai dan mengevaluasi akan memberi Anda jawaban yang Anda cari. $t=1$ $c(x,1)$

— Bill Barth
sumber

c

$c$

a

$a$

Tidak masalah. Saya mengomel setelah pertukaran komentar awal kami. Saya harap ini membantu.

— Bill Barth