Pembentukan Pola

Tujuan Pembelajaran

Setelah menyelesaikan bab ini, Anda akan mampu melakukan hal-hal berikut.

Menjelaskan mengapa himpunan parameter yang berbeda dalam satu model generik dapat menghasilkan jenis pola yang berbeda.
Menganalisis kestabilan linear solusi homogen pada sistem pembentuk pola.
Memprediksi panjang gelombang pola yang muncul di atas ambang.
Mengenali bahwa suku nonlinear yang berbeda menghasilkan jenis pola yang berbeda.
Menyusun kembali model sederhana untuk pembentukan pola vegetasi.

Ilustrasi empat panel pola garis di alam: riak pasir, rusuk kaktus, pita pada ikan, dan deret awan. Deskripsi panjang tersedia. — Gambar 10.1. Pola garis di alam: riak pasir, rusuk kaktus saguaro, pita warna pada tubuh ikan, dan struktur gulungan pada awan. *Deskripsi Gambar*

Pola (lihat Gambar 10.1), seperti garis dan bintik pada kulit hewan atau riak pasir di pantai, sangat umum dijumpai di alam dan dalam eksperimen laboratorium yang dikendalikan dengan cermat. Pola biasanya muncul dalam sistem yang didorong jauh dari kesetimbangan oleh gaya luar atau sumber energi. Ketika gaya penggerak lebih besar daripada yang diperlukan untuk mengimbangi gaya inersia, sistem menanggapinya dengan menata ulang dirinya menjadi struktur periodik yang disebut pola. Ada banyak jenis pola, seperti gulungan, persegi, atau heksagon, dan pola tersebut dapat bersifat stasioner atau dinamis. Kajian tentang sistem pembentuk pola telah menjadi bidang penelitian aktif dalam fisika, kimia, optika nonlinear, dan matematika terapan selama lebih dari setengah abadCatatan: Lihat, misalnya, buku karya P. Ball yang berjudul The Self-made Tapestry: Pattern formation in nature (Oxford, New York, 1999).. Selanjutnya, kita akan membahas secara singkat pola Turing, mengkaji sebuah model kanonik pembentuk pola yang dikenal sebagai persamaan Swift-Hohenberg kompleks, dan membahas model reaksi-difusi untuk pola vegetasi.

Pola Turing

Pada tahun 1952, A.M. TuringCatatan: A.M. Turing, The chemical basis of morphogenesis, Phil. Trans. R. Soc. London B 237, 37-72 (1952). mengajukan gagasan bahwa pola diferensiasi, seperti pola yang menentukan daerah tempat hidra menumbuhkan tentakelnya, sebenarnya merupakan pola kimia, khususnya pola yang berupa struktur periodik stasionerCatatan: Turing juga menyebutkan kemungkinan adanya pola gelombang, tetapi pola tersebut bukan fokus utama makalahnya.. Turing menjelaskan bagaimana pola semacam itu dapat muncul akibat ketakstabilan solusi homogen pada suatu sistem persamaan reaksi-difusi terkopel. Ketakstabilan ini disebut ketakstabilan yang memecahkan simetri karena struktur periodik yang tumbuh akibat ketakstabilan tersebut memecahkan invariansi translasi solusi homogen awal. Baru pada dasawarsa 1990-an pola Turing diamati dalam eksperimen kimiaCatatan: V. Castets, E. Dulos, J. Boissonade, dan P. De Kepper, Experimental evidence of a sustained standing Turing-type nonequilibrium chemical pattern, Phys. Rev. Lett. 64, 2953-2956 (1990).Catatan: Q. Ouyang dan H.L. Swinney, Transition from a uniform state to hexagonal and striped Turing patterns, Nature 352, 610-612 (1991)., berupa heksagon dan garis, tetapi juga sebagai struktur tak beraturan yang berkaitan dengan keadaan turbulensi kimiaCatatan: Q. Ouyang dan H.L. Swinney, Transition to chemical turbulence, Chaos 1, 411-420 (1991)..

Persamaan Swift-Hohenberg kompleks

Simulasi numerik pola dua dimensi. Deskripsi panjang tersedia. — Gambar 10.2. Simulasi numerik persamaan Swift-Hohenberg kompleks dengan syarat batas periodik yang memperlihatkan pola garis dan heksagon. *Deskripsi Gambar*

Kita akan mengkaji pembentukan pola melalui satu model yang disebut persamaan Swift-Hohenberg kompleks. Bergantung pada apakah koefisien model ini real atau kompleks, pola stasioner atau pola yang merambat akan teramati. Gambar 10.2 memperlihatkan pola garis dan heksagon yang dihasilkan model ini. Teori umum pembentukan polaCatatan: Untuk tinjauan, lihat misalnya M.C. Cross dan P.C. Hohenberg, Pattern formation outside of equilibrium, Rev. Mod. Phys. 65, 851-1112 (1993), serta A. C. Newell, T. Passot, dan J. Lega, Order parameter equations for patterns, Ann. Rev. Fluid Mech. 25, 399-453 (1993). menjelaskan mengapa model yang berbeda dapat menghasilkan pola serupa dan, sebagai akibatnya, mengapa sistem kimia, fisika, atau biologi yang berbeda dapat memperlihatkan pola yang tampak serupa. Teori ini juga membenarkan penggunaan model generik pembentuk pola untuk memahami mekanisme yang terlibat dalam pembentukan pola, seperti yang kita lakukan pada bagian berikutnya.

Perhatikan persamaan diferensial parsial berikut

$\displaystyle \begin{align*}\frac{\partial \psi}{\partial t} = & (\mu + i \nu) \psi - (\alpha + i \beta) \left(\Omega + \frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2}\right)^2 \psi \\ & + i \eta \left(\frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2}\right) \psi - (\gamma + i \delta) |\psi|^2 \psi - \zeta |\psi|^2, \qquad (10.1)\end{align*}$

dengan semua parameter bernilai real, $\alpha > 0$ , $\gamma > 0$ , dan $\psi(x,y,t)$ secara a priori bernilai kompleks. Persamaan ini merupakan suatu versi persamaan Swift-HohenbergCatatan: J. Swift dan P.C. Hohenberg, Hydrodynamic fluctuations at the convective instability, Phys. Rev. A 15, 319-328 (1977). dengan koefisien kompleksCatatan: J. Lega, J.V. Moloney, dan A.C. Newell, Swift-Hohenberg equation for lasers, Phys. Rev. Lett. 73, 2978-2981 (1994). dan suku kuadratik. Mudah dilihat bahwa $\psi=0$ merupakan suatu solusi, dan wajar untuk menanyakan dalam kondisi apa solusi itu juga stabil. Kita melanjutkan dengan cara yang sama seperti pada sistem persamaan diferensial biasa, yaitu melinearkan Persamaan (10.1) di sekitar solusi $\psi=0$ . Persamaan hasil linearisasinya adalah

$\displaystyle \begin{align*} \frac{\partial \psi}{\partial t} = & (\mu + i \nu) \psi - (\alpha + i \beta) \left(\Omega + \frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2}\right)^2 \psi \\ & + i \eta \left(\frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2}\right) \psi. \qquad (10.2)\end{align*}$

Pada umumnya, kita tertarik pada sistem yang membentang dalam ruang dan cukup besar untuk memungkinkan pola berkembang (makna "cukup besar" akan diperjelas sebentar lagi). Karena itu, kita hanya meninjau ketakstabilan yang berkembang di bagian dalam sistem, dan menganggap bahwa $\psi$ dapat diuraikan menjadi mode-mode Fourier. Mode Fourier $u_{\vec k} \equiv a_{\vec k}\, \exp(i \vec k \cdot \vec r)$ dengan vektor gelombang $\vec k = (k_x, k_y)$ merupakan fungsi eigen operator linear, dengan nilai eigen $\lambda_{\vec k}$ yang diberikan oleh

$\lambda_{\vec k} = \mu + i \nu - (\alpha + i \beta) (\Omega - k_x^2 - k_y^2)^2 - i \eta (k_x^2 + k_y^2).$

Dengan demikian, laju pertumbuhan mode ini adalah

$\sigma_k = \operatorname{Re}(\lambda_{\vec k}) = \mu - \alpha (\Omega - k^2)^2,$

dengan $k = \sqrt{k_x^2 + k_y^2} = ||\vec k||$ . Jadi, setiap mode Fourier tumbuh (atau meluruh) dengan laju yang hanya bergantung pada magnitudo $k$ dari vektor gelombang $\vec k$ yang bersangkutan, bukan pada arahnya. Jika $\mu < 0$ , semua mode Fourier meluruh dan solusi $\psi = 0$ stabil secara linear. Ketika $\mu$ meningkat, beberapa mode Fourier akan menjadi tak stabil.

Grafik laju pertumbuhan. Tersedia deskripsi panjang. — Gambar 10.3. Sketsa grafik laju pertumbuhan σ_k sebagai fungsi k ≥ 0, untuk 0 < μ < α Ω². (a) Ω < 0; (b) Ω > 0. *Deskripsi Gambar*

Jika $\Omega < 0$ , grafik $\sigma_k$ sebagai fungsi $k \ge 0$ bersifat monoton (lihat Gambar 10.3.a), dan mode Fourier pertama yang menjadi tak stabil adalah mode $k = 0$ . Ketakstabilan ini terjadi pada $\mu = \alpha\, \Omega^2$ . Sebaliknya, jika $\Omega > 0$ , maka $\sigma_k$ mencapai maksimum pada $k = \sqrt \Omega \equiv k_c$ (lihat Gambar 10.3.b), dan mode-mode dengan $k = k_c$ menjadi tak stabil ketika $\mu$ meningkat melewati nol. Dalam kasus ini, kita memperkirakan terbentuknya suatu pola di atas ambang dengan panjang karakteristik $l_c = 2 \pi / k_c$ . Karena itu, ukuran sistem tempat pola tersebut hendak diamati harus jauh lebih besar daripada $l_c$ .

Di atas ambang, suku-suku nonlinearlah yang menentukan pola mana yang terpilih. Biasanya, jika terdapat suku kuadratik (yaitu jika $\zeta \ne 0$ dalam Persamaan (10.1)), maka pola heksagonal diamati di dekat ambang. Sebaliknya, jika suku kubik dominan, gulungan (atau garis-garis) lebih banyak muncul. Jika bagian imajiner koefisien-koefisien dalam (10.1) dibuat nol, pola yang berkembang di atas ambang bersifat stasioner. Jika tidak, pola tersebut bergantung pada waktu. Aspek-aspek ini dan aspek lain dari dinamika Persamaan (10.1), termasuk ketakstabilan sekunder dan ketakteraturan ruang-waktu, dapat dieksplorasi dengan antarmuka GUI MATLAB bernama Patterns. Antarmuka ini memungkinkan pengguna memilih parameter dan memulai simulasi dari kondisi awal acak yang kecil. Cuplikan berkode warna dari bagian real solusi $\psi$ ditampilkan pada waktu-waktu berturut-turut selama simulasi berlangsung. Setelah simulasi selesai, parameter baru dapat dimasukkan dan simulasi baru kemudian dapat dimulai kembali dari solusi terakhir.

Teori pembentukan polaCatatan: Lihat dua artikel tinjauan yang disebutkan di atas. menyediakan cara untuk mendeskripsikan dinamika nonlinear suatu sistem pembentuk pola di dekat dan di atas ambang. Namun, pembahasan lengkap mengenai topik ini berada di luar cakupan catatan ini. Pada bagian berikutnya, kita akan menyinggung secara singkat model reaksi-difusi untuk mendeskripsikan pola vegetasi di wilayah semiarid.

Pola vegetasi

Di kawasan semiarid, tumbuhan (semak, rumput, dan sebagainya) cenderung tersusun menjadi pita-pita yang sejajar dengan garis kontur ketinggian di medan berbukit, dan menjadi petak-petak di tanah datar. Model reaksi-difusi telah diajukan untuk menjelaskan pengamatan iniCatatan: C.A. Klausmeier, Regular and irregular patterns in semiarid vegetation, Science 284, 1826-1828 (1999).Catatan: R. Lefever, O. Lejeune dan P. Couteron, Generic modelling of vegetation patterns. A case study of Tiger bush in sub-saharian Sahel, dalam Mathematical Models for Biological Pattern Formation, disunting oleh P.K. Maini dan H.G. Othmer, hlm. 83-112, Springer, New York, 2001.Catatan: J. von Hardenberg, E. Meron, M. Shachak dan Y. Zarmi, Diversity of vegetation patterns and desertification, Phys. Rev. Lett. 87, 198101 (2001).. Dalam salah satu artikel tersebut, C.A. Klausmeier memperkenalkan sebuah model sederhana dengan dua variabel terikat, yaitu biomassa tumbuhan $n$ dan jumlah air $w$ . Dalam bentuk tak berdimensi, model tersebut ditulis sebagai

$\begin{align*} \displaystyle \frac{\partial w}{\partial t} & = a - w - w\,n^2 + v \frac{\partial w}{\partial x}, \\ \displaystyle \frac{\partial n}{\partial t} & = w\,n^2 - m\, n + \left(\frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2}\right)n,\end{align*}$

di mana $a$ menyatakan masukan air melalui hujan, $v$ adalah laju pengangkutan air menuruni lereng (ke arah $-x$ ), $w\, n^2$ menyatakan pertambahan biomassa akibat konsumsi air, dan $m$ adalah laju kematian tumbuhan. Selain itu, air menguap dengan laju $w$ , dan difusi $n$ merepresentasikan penyebaran petak-petak vegetasi.

Setelah menentukan titik-titik kesetimbangan (solusi stasioner dan seragam) model ini, penulis menelaah kestabilan linearnya. Gambar 2 dalam artikel Klausmeier tahun 1999 menampilkan himpunan aras panjang gelombang pola yang diperkirakan tumbuh di atas ambang pada bidang parameter $(a,m)$ , untuk suatu nilai tetap $v$ . Lebih tepatnya, bidang $(a,m)$ dapat dibagi menjadi tiga daerah: satu daerah tempat titik tetap yang bersesuaian dengan ketiadaan vegetasi bersifat stabil, satu daerah tempat titik tetap yang bersesuaian dengan vegetasi homogen bersifat stabil, dan daerah di antaranya tempat analisis kestabilan linear memprediksi keberadaan pola vegetasi. Simulasi numerik model dalam rezim terakhir ini mengonfirmasi adanya pita-pita vegetasi yang merambat di medan berbukit dan petak-petak vegetasi di tanah datar. Laju $v$ aliran air menuruni lereng menjadi ukuran kemiringan lereng.

Ringkasan

Berbagai jenis pola dapat diperoleh dari model generik yang sama dengan memilih kumpulan parameter yang berbeda. Kita telah melihat cara menggunakan analisis kestabilan linear untuk memprediksi apakah suatu pola dapat berkembang dari solusi homogen ketika sebuah parameter kendali ( $\mu$ dalam pembahasan persamaan Swift-Hohenberg) diubah dan, jika dapat, cara menentukan skala panjang khas struktur yang akan muncul. Walaupun analisis tersebut bersifat linear, hasilnya tetap memberi petunjuk tentang jenis pola yang dapat diharapkan. Metodologi yang sama, yaitu menentukan solusi homogen dan mempelajari kestabilan linearnya, dapat diterapkan pada setiap sistem pembentuk pola, sebagaimana diilustrasikan secara singkat pada model Klausmeier.

Deskripsi Gambar

Gambar 10.1: Ilustrasi vektor independen dengan empat panel yang menunjukkan struktur berkala serupa pada objek yang berbeda. Panel kiri atas menampilkan lima riak pasir bergelombang; panel kanan atas menampilkan rusuk vertikal pada kaktus saguaro; panel kiri bawah menampilkan pita-pita horizontal pada tubuh ikan; dan panel kanan bawah menampilkan tiga deret awan putih dengan latar langit biru. Ilustrasi ini tidak memakai piksel dari kolase foto sumber. [Kembali ke Gambar 10.1]

Gambar 10.2: Simulasi numerik (dengan syarat batas periodik) persamaan Swift-Hohenberg kompleks yang memperlihatkan pola garis dua dimensi (panel kiri) dan pola heksagonal (panel kanan). Orientasi garis-garis bervariasi di seluruh gambar sebelah kiri. Cacat (ketakteraturan pada pola) tampak pada kedua gambar. Cacat pada pola garis mencakup sebuah disklinasi (yang tampak seperti lengkungan), batas-batas butir pendek (garis tempat pola-pola garis dengan orientasi berbeda bertemu), dan dislokasi (titik tempat satu garis bercabang menjadi dua). Cacat pada pola heksagonal berbentuk pasangan penta-hepta (bercak-bercak bertetangga yang memiliki 5 dan 7 tetangga). [Kembali ke Gambar 10.2]

Gambar 10.3: Dua grafik, (a) dan (b), yang menampilkan laju pertumbuhan $y = \sigma_k$ terhadap bilangan gelombang $k$ . Grafik (a) memperlihatkan kurva yang menurun secara ketat dari suatu nilai negatif, yang menunjukkan sistem stabil tempat semua gangguan meluruh. Grafik (b) memperlihatkan kurva yang naik melewati sumbu $x$ dan mencapai laju pertumbuhan maksimum $\mu$ pada $k = \sqrt{\Omega}$ , yang menunjukkan bahwa frekuensi spasial tertentu akan tumbuh dan membentuk pola. [Kembali ke Gambar 10.3]

Bahan Renungan

Soal 1

Perhatikan Persamaan (10.2) dengan $\mu = 0$ . Carilah solusi persamaan ini dalam bentuk

$\psi(x,y,t) = \exp[i(\omega t + q x + p y)], \qquad \omega,\, q,\,p\, \in \mathbb{R}.$

Cobalah menyatakan syarat keberadaan serta $\omega$ , $q$ , dan $p$ dalam parameter-parameter persamaan tersebut. Apakah solusinya tunggal?

Soal 2

Carilah solusi dalam bentuk

$\psi(x,y,t) = \exp[i(\omega t + q x + p y)], \qquad \omega,\, q,\,p\, \in \mathbb{R}.$

untuk Persamaan (10.2) dengan $\mu \ne 0$ . Cobalah menyatakan syarat keberadaan serta $\omega$ , $q$ , dan $p$ dalam parameter-parameter persamaan tersebut. Apakah solusinya tunggal?

Soal 3

Carilah solusi dalam bentuk

$\psi(x,y,t) = \exp[\lambda t + i(\omega t + q x + p y)], \qquad \lambda,\,\omega,\, q,\,p\, \in \mathbb{R}$

untuk Persamaan (10.2) dengan $\mu \ne 0$ . Cobalah menyatakan $\lambda$ , $\omega$ , $q$ , dan $p$ dalam parameter-parameter persamaan tersebut. Apakah solusinya tunggal?

Soal 4

Bacalah artikel C.A. Klausmeier tahun 1999. Jelaskan cara memperoleh bentuk tak berdimensi (2) dari sistem persamaan (1) dalam artikel tersebut.

Soal 5

Bacalah artikel C.A. Klausmeier tahun 1999. Tentukan dimensi parameter $A$ , $L$ , $R$ , $V$ , $J$ , $M$ , dan $D$ , serta variabel $W$ dan $N$ dalam sistem persamaan (1) pada artikel tersebut.

Soal 6

Bacalah artikel C.A. Klausmeier tahun 1999.

Jelaskan cara menemukan titik-titik tetap sistem persamaan (2) dalam artikel tersebut.
Tentukan matriks Jacobi sistem (2).

Petunjuk, pemeriksaan, dan pembahasan

Bagian tambahan ini ditulis untuk edisi Bahasa Indonesia. Bukalah seperlunya setelah Anda berusaha menyelesaikan soal secara mandiri.

Dukungan untuk Soal 1

Petunjuk

Tetapkan $Q=q^2+p^2$ . Operator Laplace mengalikan gelombang bidang dengan $-Q$ , sedangkan turunan waktunya adalah $i\omega\psi$ . Samakan bagian real dan imajiner dari faktor pengalinya; modulus tetap menuntut bagian real sama dengan nol.

Periksa jawaban

Pemeriksaan akhir: Dengan $Q=q^2+p^2$ , substitusi memberi $i\omega=\mu+i\nu-(\alpha+i\beta)(\Omega-Q)^2-i\eta Q$ . Untuk $\mu=0$ dan $\alpha>0$ , bagian real lenyap tepat ketika $Q=\Omega$ . Karena itu harus berlaku $\Omega\ge0$ , dan frekuensinya adalah $\omega=\nu-\eta\Omega$ . Jika $\Omega>0$ , pasangan $(q,p)$ dapat berupa titik mana pun pada lingkaran $q^2+p^2=\Omega$ , sehingga solusi tidak unik terhadap arah. Jika $\Omega=0$ , hanya $q=p=0$ . Jika $\Omega<0$ , tidak ada pasangan real $q,p$ dan tidak ada mode bermodulus satu dalam bentuk yang diminta.
Cek cepat: Periksa bahwa $q^2+p^2$ tidak pernah negatif. Pada $Q=\Omega$ , seluruh suku yang mengandung $\alpha$ dan $\beta$ harus lenyap, sehingga hasil akhir tidak bergantung pada keduanya.

Pembahasan atau rubrik

Langkah penyelesaian

Penjelasan
Definisikan kuadrat magnitudo vektor gelombang sebagai satu besaran tak negatif.
Rumus
$Q=q^2+p^2\ge0$
Penjelasan
Gunakan aksi operator diferensial pada gelombang bidang.
Rumus
$\partial_t\psi=i\omega\psi,\qquad(\partial_x^2+\partial_y^2)\psi=-Q\psi$
Penjelasan
Substitusikan ke persamaan linear dan bagi dengan $\psi\ne0$ .
Rumus
$i\omega=\mu+i\nu-(\alpha+i\beta)(\Omega-Q)^2-i\eta Q$
Penjelasan
Untuk $\mu=0$ , bagian real hanya dapat nol jika kuadratnya nol.
Rumus
$0=-\alpha(\Omega-Q)^2\quad\Longrightarrow\quad Q=\Omega$
Penjelasan
Ambil bagian imajiner setelah syarat netral dipenuhi.
Rumus
$\omega=\nu-\beta(\Omega-Q)^2-\eta Q=\nu-\eta\Omega$

Simpulan

Mode dengan $q,p\in\mathbb R$ ada hanya untuk $\Omega\ge0$ ; lingkaran arah membuatnya tidak unik ketika $\Omega>0$ , sedangkan $\Omega=0$ hanya memberi mode spasial seragam.

Unduh notebook Python terbuka untuk Soal 1

Dukungan untuk Soal 2

Petunjuk

Gunakan kembali $Q=q^2+p^2\ge0$ . Karena ansatz tidak memiliki faktor pertumbuhan $e^{\lambda t}$ , bagian real nilai eigen harus nol. Selesaikan kuadrat $(\Omega-Q)^2=\mu/\alpha$ , lalu saring setiap cabang dengan syarat $Q\ge0$ .

Periksa jawaban

Pemeriksaan akhir: Mode bermodulus satu menuntut $0=\mu-\alpha(\Omega-Q)^2$ . Karena $\alpha>0$ dan $\mu\ne0$ , solusi dengan $q,p\in\mathbb R$ hanya mungkin untuk $\mu>0$ . Tuliskan $s=\sqrt{\mu/\alpha}$ ; cabang radialnya $Q=\Omega\pm s$ , dengan hanya $Q\ge0$ yang diperbolehkan. Jika $\Omega<-s$ , tidak ada solusi. Jika $\Omega=-s$ , satu-satunya cabang ialah $Q=0$ , sehingga $(q,p)=(0,0)$ dan solusinya unik. Untuk setiap $Q>0$ , pasangan $(q,p)$ mengisi lingkaran $q^2+p^2=Q$ , sehingga solusi tidak unik; hal ini mencakup $\Omega=s$ , yang juga memiliki cabang $Q=0$ selain lingkaran positif. Pada setiap cabang sah, $\omega=\nu-\beta\mu/\alpha-\eta Q$ . Untuk $\mu<0$ , tidak ada gelombang bidang bermodulus satu dengan $q,p\in\mathbb R$ dalam bentuk ini.
Cek cepat: Masukkan setiap kandidat $Q$ ke $(\Omega-Q)^2$ : hasilnya harus $\mu/\alpha$ . Buang cabang dengan $Q<0$ , dan jangan menyatakan adanya solusi untuk $\mu<0$ .

Pembahasan atau rubrik

Langkah penyelesaian

Penjelasan
Pisahkan faktor pengali menjadi bagian real dan imajiner.
Rumus
$i\omega=\underbrace{\mu-\alpha(\Omega-Q)^2}_{\text{real}}+i\underbrace{[\nu-\beta(\Omega-Q)^2-\eta Q]}_{\text{imajiner}}$
Penjelasan
Modulus satu mengharuskan laju pertumbuhan real nol.
Rumus
$(\Omega-Q)^2=\frac{\mu}{\alpha}$
Penjelasan
Akar real memerlukan $\mu>0$ , lalu menghasilkan dua cabang aljabar yang harus diuji terhadap $Q\ge0$ .
Rumus
$Q_{\pm}=\Omega\pm\sqrt{\frac{\mu}{\alpha}},\qquad Q_{\pm}\ge0$
Penjelasan
Gunakan syarat real tadi untuk menyederhanakan bagian imajiner.
Rumus
$\omega_{\pm}=\nu-\frac{\beta\mu}{\alpha}-\eta Q_{\pm}$
Penjelasan
Untuk radius positif, parameter sudut memperlihatkan keluarga arah yang kontinu.
Rumus
$q=\sqrt Q\cos\theta,\qquad p=\sqrt Q\sin\theta,\qquad0\le\theta<2\pi$

Simpulan

Keberadaan ditentukan oleh tanda $\mu$ dan ketaknegatifan setiap cabang $Q$ . Solusi unik hanya terjadi pada batas $\mu>0,\ \Omega=-\sqrt{\mu/\alpha}$ , ketika $Q=0$ ; setiap cabang positif menghasilkan ketidakunikan arah, dan dua cabang sah dapat menambah ketidakunikan radial.

Unduh notebook Python terbuka untuk Soal 2

Dukungan untuk Soal 3

Petunjuk

Sekarang turunan waktu menghasilkan $(\lambda+i\omega)\psi$ , sehingga bagian real tidak perlu dibuat nol. Tetapkan $Q=q^2+p^2$ dan cocokkan langsung bagian real serta imajiner dari nilai eigen Fourier.

Periksa jawaban

Pemeriksaan akhir: Untuk $Q=q^2+p^2\ge0$ , substitusi memberi $\lambda+i\omega=\mu+i\nu-(\alpha+i\beta)(\Omega-Q)^2-i\eta Q$ . Jadi $\lambda=\mu-\alpha(\Omega-Q)^2$ dan $\omega=\nu-\beta(\Omega-Q)^2-\eta Q$ . Setiap $q,p\in\mathbb R$ mendefinisikan satu mode; untuk $Q>0$ , semua arah dengan radius $\sqrt Q$ memiliki $\lambda$ dan $\omega$ yang sama. Mode tumbuh jika $\lambda>0$ , netral jika $\lambda=0$ , dan meluruh jika $\lambda<0$ .
Cek cepat: Bagian real harus persis sama dengan kurva pertumbuhan bab, $\sigma_Q=\mu-\alpha(\Omega-Q)^2$ . Menetapkan $\lambda=0$ harus mengembalikan syarat Soal 1 dan 2.

Pembahasan atau rubrik

Langkah penyelesaian

Penjelasan
Definisikan kuadrat magnitudo vektor gelombang.
Rumus
$Q=q^2+p^2$
Penjelasan
Hitung turunan waktu dan aksi Laplasian.
Rumus
$\partial_t\psi=(\lambda+i\omega)\psi,\qquad\Delta\psi=-Q\psi$
Penjelasan
Operator Swift-Hohenberg kuadrat menjadi pengali skalar.
Rumus
$(\Omega+\Delta)^2\psi=(\Omega-Q)^2\psi$
Penjelasan
Samakan faktor kompleks pada kedua ruas.
Rumus
$\lambda+i\omega=\mu+i\nu-(\alpha+i\beta)(\Omega-Q)^2-i\eta Q$
Penjelasan
Pisahkan bagian real dan imajiner.
Rumus
$\lambda=\mu-\alpha(\Omega-Q)^2,\qquad\omega=\nu-\beta(\Omega-Q)^2-\eta Q$

Simpulan

Relasi dispersi bergantung hanya pada besar vektor gelombang, bukan arahnya; $\lambda$ menentukan pertumbuhan amplitudo dan $\omega$ menentukan perubahan fase.

Unduh notebook Python terbuka untuk Soal 3

Dukungan untuk Soal 4

Petunjuk

Kartu persamaan dan penskalaan—dukungan baru, dengan fakta model dirujuk ke Klausmeier, DOI 10.1126/science.284.5421.1826: $W_T=A-LW-RWN^2+VW_X$ , $N_T=JRWN^2-MN+D\Delta N$ . Gunakan $t=LT$ , $(x,y)=\sqrt{L/D}(X,Y)$ , $n=N\sqrt{R/L}$ , $w=JW\sqrt{R/L}$ , serta $a=AJ\sqrt R/L^{3/2}$ , $m=M/L$ , dan $v=V/\sqrt{LD}$ . Tuliskan aturan rantai sebelum menyederhanakan setiap suku.

Periksa jawaban

Pemeriksaan akhir: Kartu mandiri yang dirujuk ke Klausmeier (DOI 10.1126/science.284.5421.1826) menggunakan $W_T=A-LW-RWN^2+VW_X$ dan $N_T=JRWN^2-MN+D\Delta_{X,Y}N$ . Tetapkan $s=\sqrt{R/L}$ , lalu $t=LT$ , $x=\sqrt{L/D}X$ , $y=\sqrt{L/D}Y$ , $n=sN$ , dan $w=JsW$ . Dengan $a=AJs/L=AJ\sqrt R/L^{3/2}$ , $m=M/L$ , dan $v=V/\sqrt{LD}$ , aturan rantai menghasilkan tepat $w_t=a-w-wn^2+vw_x$ dan $n_t=wn^2-mn+\Delta_{x,y}n$ . Variabel $t,x,y,n,w$ dan parameter $a,m,v$ yang baru semuanya tak berdimensi.
Cek cepat: Koefisien $w$ , $wn^2$ , dan $\Delta n$ harus menjadi satu. Koefisien adveksi harus menyederhana menjadi $V/\sqrt{LD}$ , bukan $V\sqrt{L/D}$ .

Pembahasan atau rubrik

Langkah penyelesaian

Penjelasan
Kartu sumber-fakta untuk perhitungan ini adalah sistem dimensional Klausmeier; uraian dan penurunan berikut merupakan dukungan baru.
Rumus
$W_T=A-LW-RWN^2+VW_X,\qquad N_T=JRWN^2-MN+D\Delta_{X,Y}N$
Penjelasan
Perkenalkan singkatan skala amplitudo dan balikkan definisi variabel terikat.
Rumus
$s=\sqrt{R/L},\qquad W=\frac{w}{Js},\qquad N=\frac n s$
Penjelasan
Transformasi koordinat memberi faktor turunan yang diperlukan.
Rumus
$\partial_T=L\partial_t,\qquad\partial_X=\sqrt{L/D}\partial_x,\qquad\Delta_{X,Y}=\frac LD\Delta_{x,y}$
Penjelasan
Substitusi ke persamaan air dan perkalian dengan $Js/L$ menormalkan peluruhan serta konsumsi.
Rumus
$w_t=\frac{AJs}{L}-w-wn^2+\frac{V}{\sqrt{LD}}w_x$
Penjelasan
Substitusi ke persamaan biomassa dan perkalian dengan $s/L$ menormalkan pertumbuhan serta difusi.
Rumus
$n_t=wn^2-\frac ML n+\Delta_{x,y}n$
Penjelasan
Namai tiga kombinasi parameter yang tersisa.
Rumus
$a=\frac{AJ\sqrt R}{L^{3/2}},\qquad m=\frac ML,\qquad v=\frac V{\sqrt{LD}}$

Simpulan

Penskalaan mengubah sistem dimensional menjadi $w_t=a-w-wn^2+vw_x$ , $n_t=wn^2-mn+\Delta n$ , dengan hanya tiga parameter tak berdimensi $a,m,v$ .

Unduh notebook Python terbuka untuk Soal 4

Dukungan untuk Soal 5

Petunjuk

Kartu persamaan—dukungan baru, dengan fakta model dirujuk ke Klausmeier, DOI 10.1126/science.284.5421.1826: $W_T=A-LW-RWN^2+VW_X$ dan $N_T=JRWN^2-MN+D\Delta N$ . Secara fisik, $W$ adalah massa air per satuan luas dan $N$ adalah massa kering biomassa tumbuhan per satuan luas. Misalkan $[X]=\mathsf L$ dan $[T]=\mathsf T$ , lalu samakan dimensi setiap suku dengan dimensi ruas kirinya.

Periksa jawaban

Pemeriksaan akhir: Dengan $\mathsf L$ sebagai dimensi panjang dan $\mathsf T$ sebagai dimensi waktu, $W$ berdimensi massa air per satuan luas (misalnya kg H₂O m⁻²) dan $N$ berdimensi massa kering biomassa tumbuhan per satuan luas (misalnya kg massa kering m⁻²). Hasil parameter ialah $[A]=[W]\mathsf T^{-1}$ , $[L]=[M]=\mathsf T^{-1}$ , $[R]=[N]^{-2}\mathsf T^{-1}$ , $[V]=\mathsf L\mathsf T^{-1}$ , $[J]=[N][W]^{-1}$ , dan $[D]=\mathsf L^2\mathsf T^{-1}$ .
Cek cepat: Setiap suku pada persamaan $W$ harus berdimensi $[W]/\mathsf T$ , dan setiap suku pada persamaan $N$ harus berdimensi $[N]/\mathsf T$ . Kombinasi $a,m,v,n,w$ pada penskalaan Soal 4 semuanya harus tak berdimensi.

Pembahasan atau rubrik

Langkah penyelesaian

Penjelasan
Suku sumber pada persamaan air harus mempunyai dimensi yang sama dengan $W_T$ .
Rumus
$[A]=[W]\mathsf T^{-1}$
Penjelasan
Peluruhan linear pada kedua persamaan menetapkan dimensi $L$ dan $M$ .
Rumus
$[LW]=[W]\mathsf T^{-1},\ [MN]=[N]\mathsf T^{-1}\quad\Longrightarrow\quad[L]=[M]=\mathsf T^{-1}$
Penjelasan
Samakan suku konsumsi air dengan $W_T$ .
Rumus
$[R][W][N]^2=[W]\mathsf T^{-1}\quad\Longrightarrow\quad[R]=[N]^{-2}\mathsf T^{-1}$
Penjelasan
Turunan ruang pertama dan kedua menetapkan dimensi adveksi dan difusi.
Rumus
$[V]\frac{[W]}{\mathsf L}=\frac{[W]}{\mathsf T},\qquad[D]\frac{[N]}{\mathsf L^2}=\frac{[N]}{\mathsf T}$
Penjelasan
Gunakan suku produksi biomassa untuk menentukan faktor hasil $J$ .
Rumus
$[J][R][W][N]^2=[N]\mathsf T^{-1}\quad\Longrightarrow\quad[J]=[N][W]^{-1}$
Penjelasan
Kumpulkan hasil adveksi dan difusi.
Rumus
$[V]=\mathsf L\mathsf T^{-1},\qquad[D]=\mathsf L^2\mathsf T^{-1}$

Simpulan

Dengan $W$ sebagai massa air per luas dan $N$ sebagai massa kering biomassa per luas, dimensi tersebut membuat kedua persamaan homogen dan membuktikan bahwa semua kelompok pada penskalaan Soal 4 benar-benar tanpa satuan.

Unduh notebook Python terbuka untuk Soal 5

Dukungan untuk Soal 6

Petunjuk

Untuk keadaan spasial homogen, turunan ruang lenyap. Tetapkan $f(w,n)=a-w-wn^2$ dan $g(w,n)=wn^2-mn=n(wn-m)$ . Pisahkan cabang $n=0$ dari cabang $wn=m$ , kemudian hitung Jacobi dalam urutan variabel $(w,n)$ .

Periksa jawaban

Pemeriksaan akhir: Titik tetap tanpa vegetasi adalah $(w,n)=(a,0)$ . Untuk vegetasi positif, $wn=m$ dan $mn^2-an+m=0$ , sehingga—untuk $a,m>0$ —titik vegetasi ada tepat ketika $a\ge2m$ : $n_{\pm}=[a\pm\sqrt{a^2-4m^2}]/(2m)$ dan $w_{\pm}=m/n_{\pm}$ . Matriks Jacobi bagian reaksi dalam urutan $(w,n)$ adalah $J(w,n)=\begin{pmatrix}-(1+n^2)&-2wn\\n^2&2wn-m\end{pmatrix}$ . Pada titik vegetasi, $wn=m$ , jadi $J_{\pm}=\begin{pmatrix}-(1+n_{\pm}^2)&-2m\\n_{\pm}^2&m\end{pmatrix}$ .
Cek cepat: Substitusi $(a,0)$ harus menolkan kedua reaksi. Untuk setiap akar vegetasi, hasil kali $n_+n_-=1$ dan $w_{\pm}n_{\pm}=m$ ; diskriminannya harus $a^2-4m^2$ .

Pembahasan atau rubrik

Langkah penyelesaian

Penjelasan
Hilangkan turunan spasial dan definisikan medan reaksi.
Rumus
$f(w,n)=a-w-wn^2,\qquad g(w,n)=wn^2-mn=n(wn-m)$
Penjelasan
Cabang tanpa vegetasi berasal dari $n=0$ .
Rumus
$n=0,\quad f=0\Longrightarrow w=a$
Penjelasan
Pada cabang vegetasi, gunakan $wn=m$ dalam persamaan air.
Rumus
$w=\frac mn,\qquad0=a-\frac mn-mn\quad\Longrightarrow\quad mn^2-an+m=0$
Penjelasan
Selesaikan persamaan kuadrat dan nyatakan syarat keberadaan real positif.
Rumus
$n_{\pm}=\frac{a\pm\sqrt{a^2-4m^2}}{2m},\qquad w_{\pm}=\frac m{n_{\pm}},\qquad a\ge2m$
Penjelasan
Diferensiasikan kedua fungsi reaksi dalam urutan $(w,n)$ .
Rumus
$J(w,n)=\begin{pmatrix}f_w&f_n\\g_w&g_n\end{pmatrix}=\begin{pmatrix}-(1+n^2)&-2wn\\n^2&2wn-m\end{pmatrix}$
Penjelasan
Pada titik vegetasi, sederhanakan dengan hubungan $wn=m$ .
Rumus
$J_{\pm}=\begin{pmatrix}-(1+n_{\pm}^2)&-2m\\n_{\pm}^2&m\end{pmatrix}$

Simpulan

Sistem mempunyai satu keadaan tanpa vegetasi untuk semua $a>0$ . Cabang vegetasi ada ketika $a\ge2m$ : terdapat dua keadaan berbeda untuk $a>2m$ , sedangkan pada $a=2m$ kedua akar vegetasi berimpit.

Unduh notebook Python terbuka untuk Soal 6