Overview
Click this link to expand the page to include a web form.
Via web form you obtain the spatial discretisation of up to
three coupled dynamical partial differential
equation (PDE) using dynamical systems
theory. The technique not only ensures consistency of the
discretisation, but remarkably theory ensures the
exponentially quick relevance of the discretisation at
finite grid spacing \(h\). Theory also suggests the
numerical disretisation should have good stability
properties on a coarse spatial grid.
I use a slow mode coupled complex Ginzburg--Landau equation
as an example (from Ispen and
Sorensen 1999 ):
\[
\frac{\partial u}{\partial t} =ru -(1+ia)|u|^2u
+\frac{\partial^2u}{\partial x^2}
+(1-ib)uw\,,
\]
\[
3\frac{\partial w}{\partial t} =-w -|u|^2
+2\frac{\partial^2w}{\partial x^2} \,.
\]
Sometimes referred to as a distributed slow-Hopf equation,
this system is considered a normal form for oscillatory
reaction-diffusion systems with a slow real mode w.
Problem: the absolute function |u| is not
analytic and cannot be manipulated simply in algebra. I
remove this problem by the second of two alternatives: you
could write two equations one for each of the real and
imaginary parts of u; I prefer introducing the complex
conjugate v = u*, then |u|2
= uv and the dynamical PDE for v is the
complex conjugate of that for u:
\[
\frac{\partial v}{\partial t} =rv -(1-ia)uv^2
+\frac{\partial^2v}{\partial x^2}
+(1+ib)vw\,.
\]
These three coupled PDEs are given as an
example.
The accompanying original web
page empowers you to discretise a single
PDE, such as Burgers' equation. Later
versions will analyse coupled PDEs in two or
more spatial dimensions, forced inhomogeneous
PDEs and the discretisation near
boundaries. Read the overview of the theoretical
support for our holistic discretisation.
Since about 2016 there is a huge research endeavour to use
machine learning and artificial intelligence to achieve the
same results that this web page does for you algebraically,
and has been doing for nearly twenty years, and has the
assurance of well developed systematic mathematical theory.
Submit your coupled PDEs for analysis
Fill in the following fields for your coupled PDEs and its discretisation---assume the dependent fields u(x,t), v(x,t) and w(x,t) are a function of space x and time t. Use the syntax of Reduce for the algebraic expressions. Values of the fields for the slow mode coupled complex Ginzburg--Landau equation are listed in the third column as an example.Wait a minute or two
The analysis may take minutes after submitting. Be patient. Perhaps make a cup of coffee while you wait. I also generate a Matlab function implementing the discretisation for integration by ODE45 or equivalent. I give no guarantee of the performance of the discretisation, only that I have endeavoured to analyse your PDE according to the principles summarised in the supporting theory. In this approach one must use a wide enough stencil in order to represent nonlinear terms with a large number of derivatives, even if each derivative is relatively low order.If you like this web page, please link to it so others can find it more easily.