WolframAlpha.com
WolframCloud.com
All Sites & Public Resources...
Products & Services
Wolfram|One
Mathematica
Wolfram|Alpha Notebook Edition
Programming Lab
Finance Platform
SystemModeler
Wolfram Player
Wolfram Engine
WolframScript
Enterprise Private Cloud
Enterprise Mathematica
Wolfram|Alpha Appliance
Enterprise Solutions
Corporate Consulting
Technical Consulting
Wolfram|Alpha Business Solutions
Resource System
Data Repository
Neural Net Repository
Function Repository
Wolfram|Alpha
Wolfram|Alpha Pro
Problem Generator
API
Data Drop
Products for Education
Mobile Apps
Wolfram Player
Wolfram Cloud App
Wolfram|Alpha for Mobile
Wolfram|Alpha-Powered Apps
Services
Paid Project Support
Wolfram U
Summer Programs
All Products & Services »
Technologies
Wolfram Language
Revolutionary knowledge-based programming language.
Wolfram Cloud
Central infrastructure for Wolfram's cloud products & services.
Wolfram Science
Technology-enabling science of the computational universe.
Wolfram Notebooks
The preeminent environment for any technical workflows.
Wolfram Engine
Software engine implementing the Wolfram Language.
Wolfram Natural Language Understanding System
Knowledge-based broadly deployed natural language.
Wolfram Data Framework
Semantic framework for real-world data.
Wolfram Universal Deployment System
Instant deployment across cloud, desktop, mobile, and more.
Wolfram Knowledgebase
Curated computable knowledge powering Wolfram|Alpha.
All Technologies »
Solutions
Engineering, R&D
Aerospace & Defense
Chemical Engineering
Control Systems
Electrical Engineering
Image Processing
Industrial Engineering
Mechanical Engineering
Operations Research
More...
Finance, Statistics & Business Analysis
Actuarial Sciences
Bioinformatics
Data Science
Econometrics
Financial Risk Management
Statistics
More...
Education
All Solutions for Education
Trends
Machine Learning
Multiparadigm Data Science
Internet of Things
High-Performance Computing
Hackathons
Software & Web
Software Development
Authoring & Publishing
Interface Development
Web Development
Sciences
Astronomy
Biology
Chemistry
More...
All Solutions »
Learning & Support
Learning
Wolfram Language Documentation
Fast Introduction for Programmers
Wolfram U
Videos & Screencasts
Wolfram Language Introductory Book
Webinars & Training
Summer Programs
Books
Need Help?
Support FAQ
Wolfram Community
Contact Support
Premium Support
Premier Service
Technical Consulting
All Learning & Support »
Company
About
Company Background
Wolfram Blog
Events
Contact Us
Work with Us
Careers at Wolfram
Internships
Other Wolfram Language Jobs
Initiatives
Wolfram Foundation
MathWorld
Computer-Based Math
A New Kind of Science
Wolfram Technology for Hackathons
Student Ambassador Program
Wolfram for Startups
Demonstrations Project
Wolfram Innovator Awards
Wolfram + Raspberry Pi
Summer Programs
More...
All Company »
Search
Join
Sign In
Dashboard
Groups
People
Message Boards
Answer
(
Unmark
)
Mark as an Answer
GROUPS:
Staff Picks
Mathematics
Calculus
Equation Solving
Wolfram Language
Computer-Based Maths
Wolfram High School Summer Camp
5
Rauan Kaldybaev
[WSC20] Exploring staircase solutions to nonlinear wave equations
Rauan Kaldybaev
Posted
6 months ago
1438 Views
|
2 Replies
|
5 Total Likes
Follow this post
|
Introduction
F
o
r
s
o
m
e
f
[
u
]
,
n
o
n
l
i
n
e
a
r
w
a
v
e
e
q
u
a
t
i
o
n
s
o
f
t
h
e
f
o
r
m
2
∂
u
∂
2
t
=
2
∂
u
∂
2
x
+
f
[
u
]
h
a
v
e
s
o
l
u
t
i
o
n
s
t
h
a
t
l
o
o
k
l
i
k
e
s
t
a
i
r
c
a
s
e
s
,
e
.
g
.
f
e
a
t
u
r
e
s
e
v
e
r
a
l
d
i
s
t
i
n
c
t
l
e
v
e
l
s
s
e
p
a
r
a
t
e
d
b
y
s
t
e
e
p
i
n
c
l
i
n
e
s
.
T
h
e
f
o
r
m
a
t
i
o
n
o
f
t
h
e
s
e
"
s
t
a
i
r
c
a
s
e
s
"
d
e
p
e
n
d
s
o
n
i
n
i
t
i
a
l
c
o
n
d
i
t
i
o
n
s
a
n
d
t
h
e
f
o
r
m
o
f
f
[
u
]
.
T
h
e
b
o
r
d
e
r
b
e
t
w
e
e
n
t
w
o
l
e
v
e
l
s
p
r
o
p
a
g
a
t
e
s
w
i
t
h
t
h
e
s
p
e
e
d
a
p
p
r
o
x
i
m
a
t
e
l
y
e
q
u
a
l
t
o
1
,
a
n
d
g
e
t
s
m
o
r
e
a
n
d
m
o
r
e
s
h
a
r
p
w
i
t
h
t
i
m
e
.
T
h
e
r
e
g
i
o
n
o
f
t
h
e
s
t
a
i
r
c
a
s
e
s
o
l
u
t
i
o
n
t
h
a
t
i
s
i
n
s
i
d
e
t
h
e
b
o
r
d
e
r
s
t
h
e
d
e
v
i
a
t
i
o
n
o
f
u
f
r
o
m
t
h
e
n
e
a
r
e
s
t
e
q
u
i
l
i
b
r
i
u
m
i
s
s
m
a
l
l
,
s
o
t
h
e
r
e
t
h
e
e
q
u
a
t
i
o
n
c
a
n
b
e
a
p
p
r
o
x
i
m
a
t
e
d
a
s
l
i
n
e
a
r
.
General information about nonlinear wave equations
This chapter discusses some general properties of nonlinear wave equations.
Brief overview
Nonlinear wave equation are equations of the form
2
∂
u
∂
2
t
=
2
∂
u
∂
2
x
+
f
[
u
]
(
1
)
Where
f[u]
is a real-valued function and
x
goes from -
∞
to +
∞
The equation (1) is a nonlinear partial differential equation. It is described by the Lagrangian and conserves the energy function, which puts restrictions on its' solutions, and in particular means that
u
,
u
t
and
u
x
have to remain finite. For any given value of energy for a certain type of
f[u],
the value of
u
is bounded to lie in a certain region. The equation can exhibit stationary, oscillatory, unbounded, or singular behavior, depending on initial conditions and the form of
f[u]
. The equation allows waves to travel along the x-axis, and for the vast majority of
f[u]
the speed of the waves depends on the waves themselves.
Numerical solutions
Let's have a look at several numerical solutions to (1) to get the feel of it. I am going to use the NDSolve function and plot solutions using Plot3D.
Firstly, let us specify the initial conditions. The choice is arbitrary, however I'm going to start with a Gaussian. This function is defined for all x, symmetric with respect to
x = 0,
and decays quickly while being infinitely differentiable for all values of
x
.
i
n
i
t
C
o
n
d
s
=
u
[
0
,
x
]
E
x
p
[
-
x
^
2
]
,
D
e
r
i
v
a
t
i
v
e
[
1
,
0
]
[
u
]
[
0
,
x
]
0
;
I
n
[
]
:
=
G
e
n
e
r
a
l
l
y
,
u
[
t
,
x
]
d
e
n
o
t
e
s
t
h
e
v
a
l
u
e
o
f
u
a
t
t
i
m
e
t
a
n
d
p
o
s
i
t
i
o
n
x
,
a
n
d
D
e
r
i
v
a
t
i
v
e
[
n
,
m
]
[
u
]
[
t
,
x
]
i
s
t
h
e
s
a
m
e
a
s
2
∂
u
n
∂
t
m
∂
x
(
t
,
x
)
.
u
[
0
,
x
]
a
n
d
D
e
r
i
v
a
t
i
v
e
[
1
,
0
]
[
u
]
[
0
,
x
]
r
e
p
r
e
s
e
n
t
t
h
e
i
n
i
t
i
a
l
v
a
l
u
e
s
o
f
u
a
n
d
u
t
,
r
e
s
p
e
c
t
i
v
e
l
y
.
H
e
r
e
a
n
d
i
n
t
h
e
r
e
s
t
o
f
t
h
e
a
r
t
i
c
l
e
u
t
i
s
a
s
h
o
r
t
n
o
t
a
t
i
o
n
f
o
r
∂
u
∂
t
.
S
i
m
i
l
a
r
l
y
,
u
x
≡
∂
u
∂
t
,
u
t
t
≡
2
∂
u
∂
2
t
,
a
n
d
u
t
x
≡
2
∂
u
∂
t
∂
x
.
Below is the plot of initial values of u. As you can see, it decays very quickly, and even for x = 3 it’s already almost zero.
P
l
o
t
[
E
x
p
[
-
x
^
2
]
,
{
x
,
-
3
,
3
}
]
Now, let us define the boundaries. Ideally, (1) should be solved within an infinite interval of x. However, for practical reasons, let us instead solve the equation in the finite interval
x
∈
[ -R, R ]
and
t
∈
[ 0, T ],
where
R
and
T
some positive numbers. I will take
T
=
2
0
;
R
=
3
2
;
I
n
[
]
:
=
Finally, we have the boundary conditions. In the original problem, the domain of x is infinite, and we don't have to worry about those. In the numerical approximation, however, boundary conditions are compulsory. I will choose them to be
u
x
(
t
,
-
R
)
=
u
x
(
t
,
R
)
=
0
, for a number of reasons. Firstly, they are simple and familiar. They will not introduce any unexpected motion or disturbance. Secondly, they serve as a decent approximation to the missing semi-infinite string. Thirdly, as we will discover soon, these boundary conditions will prevent energy flux in or our of the studied interval.
b
o
u
n
d
a
r
y
C
o
n
d
s
=
D
e
r
i
v
a
t
i
v
e
[
0
,
1
]
[
u
]
[
t
,
-
R
]
0
,
D
e
r
i
v
a
t
i
v
e
[
0
,
1
]
[
u
]
[
t
,
R
]
0
;
I
n
[
]
:
=
Lastly, we have our equation. In general, it is written as
e
q
n
T
e
m
p
l
a
t
e
=
(
D
e
r
i
v
a
t
i
v
e
[
2
,
0
]
[
u
]
[
t
,
x
]
D
e
r
i
v
a
t
i
v
e
[
0
,
2
]
[
u
]
[
t
,
x
]
+
f
[
u
[
t
,
x
]
]
)
;
I
n
[
]
:
=
When we will need to study the case of a specific f[u], we can just use replacement like that:
e
q
n
T
e
m
p
l
a
t
e
/
.
f
S
i
n
Or like that:
e
q
n
T
e
m
p
l
a
t
e
/
.
f
F
u
n
c
t
i
o
n
[
u
,
u
-
u
^
3
]
The
Function[u, func]
thing means that we are defining a function called
func
with the argument
u
. Note:
Function[u, u - u^3]
gives the same result as
Function[x, x - x^3]
.
This is how we use NDSolve to obtain a numerical solution for (1):
s
o
l
n
1
=
N
D
S
o
l
v
e
[
J
o
i
n
[
i
n
i
t
C
o
n
d
s
,
b
o
u
n
d
a
r
y
C
o
n
d
s
,
{
e
q
n
T
e
m
p
l
a
t
e
/
.
f
F
u
n
c
t
i
o
n
[
u
,
-
u
]
}
]
,
u
,
{
t
,
0
,
T
}
,
{
x
,
-
R
,
R
}
]
(
*
t
h
e
c
u
r
l
y
b
r
a
c
k
e
t
s
a
r
o
u
n
d
f
[
u
[
t
,
x
.
.
.
a
r
e
h
e
r
e
b
e
c
a
u
s
e
o
f
t
h
e
r
u
l
e
s
o
f
W
o
l
f
r
a
m
l
a
n
g
u
a
g
e
;
t
h
e
y
d
o
n
'
t
c
h
a
n
g
e
t
h
e
m
e
a
n
i
n
g
*
)
u
I
n
t
e
r
p
o
l
a
t
i
n
g
F
u
n
c
t
i
o
n
D
o
m
a
i
n
:
{
{
0
.
,
2
0
.
}
,
{
-
3
2
.
,
3
2
.
}
}
O
u
t
p
u
t
:
s
c
a
l
a
r
O
u
t
[
]
=
Now, we can plot our solution, just like that:
P
l
o
t
3
D
[
u
[
t
,
x
]
/
.
s
o
l
n
1
,
{
t
,
0
,
T
}
,
{
x
,
-
R
,
R
}
,
P
l
o
t
R
a
n
g
e
A
l
l
,
P
l
o
t
P
o
i
n
t
s
5
0
]
For future use, I've defined a function that automatically uses NDSolve to give numerical solutions, removing the need to specify a ton of arguments every time we want to solve our beloved equation.
s
o
l
v
e
T
h
e
P
D
E
[
r
e
p
l
a
c
e
m
e
n
t
R
u
l
e
F
o
r
F
_
R
u
l
e
]
:
=
N
D
S
o
l
v
e
[
J
o
i
n
[
i
n
i
t
C
o
n
d
s
,
b
o
u
n
d
a
r
y
C
o
n
d
s
,
{
e
q
n
T
e
m
p
l
a
t
e
/
.
r
e
p
l
a
c
e
m
e
n
t
R
u
l
e
F
o
r
F
}
]
,
u
,
{
t
,
0
,
T
}
,
{
x
,
-
R
,
R
}
]
;
I
n
[
]
:
=
Here is how it works:
s
o
l
v
e
T
h
e
P
D
E
[
f
C
o
s
]
Also, here is the plotter function:
O
p
t
i
o
n
s
[
m
y
P
l
o
t
]
=
{
R
a
s
t
e
r
i
z
e
T
r
u
e
}
;
m
y
P
l
o
t
[
s
o
l
u
t
i
o
n
T
o
P
l
o
t
_
,
O
p
t
i
o
n
s
P
a
t
t
e
r
n
[
]
]
:
=
M
o
d
u
l
e
[
{
p
l
t
}
,
p
l
t
=
P
l
o
t
3
D
[
u
[
t
,
x
]
/
.
s
o
l
u
t
i
o
n
T
o
P
l
o
t
,
{
t
,
0
,
T
}
,
{
x
,
-
R
,
R
}
,
P
l
o
t
R
a
n
g
e
A
l
l
,
P
l
o
t
P
o
i
n
t
s
7
5
]
;
I
f
[
O
p
t
i
o
n
V
a
l
u
e
[
"
R
a
s
t
e
r
i
z
e
"
]
,
R
a
s
t
e
r
i
z
e
[
p
l
t
]
,
p
l
t
]
]
;
I
n
[
]
:
=
This "Rasterize" option controls whether the plot will be converted to an image or not. Converting dynamic plots to static images allows to reduce the load on the computer, however sometimes we might want to be able to spin and play around with the original plot.
Lastly, in the rest of this article I'll be working mainly not with
f[u]
but with
V[u]
, as it gives more insight into the problem. NDSolve, however, only accepts actual differential equations and thus needs
f[u]
. Thus, I've created a function to calculate
f[u]
from
V[u]
.
V
t
o
F
[
V
r
u
l
e
_
]
:
=
W
i
t
h
[
{
e
x
p
r
=
-
D
[
(
V
/
.
V
r
u
l
e
)
[
u
]
,
u
]
}
,
T
o
E
x
p
r
e
s
s
i
o
n
[
"
F
u
n
c
t
i
o
n
[
u
,
"
<
>
T
o
S
t
r
i
n
g
[
I
n
p
u
t
F
o
r
m
[
e
x
p
r
]
]
<
>
"
]
"
]
]
It looks complicated, but basically what it does is take the derivative and multiplies the result by -1. Here is an example.
V
t
o
F
[
V
F
u
n
c
t
i
o
n
[
u
,
u
^
3
]
]
The derivative of
3
u
with respect to
u
is
3
2
u
, and
f = -V' =
-
3
2
u
.
Lagrangian
As Stephen Wolfram has shown, equation (1) can be described by the Lagrangian with the density
ℒ
=
2
u
t
-
2
u
x
2
-
V
[
u
]
(
2
)
where
V
[
u
]
=
-
u
∫
0
f
[
u
'
]
d
u
'
(
3
)
Symmetry, Noether’s theorem, and conservation laws
Time translation symmetry
It can be seen that (2) is translationally symmetric with respect to t and x. The invariant corresponding to time symmetry, energy, is, as Stephen Wolfram has shown,
E
=
∞
∫
-
∞
ρ
d
x
,
d
E
d
t
=
0
,
ρ
=
2
u
t
+
2
u
x
2
+
V
[
u
]
Space translation symmetry
However, because (2) is also translationally symmetric with respect to
x
, there is one more energy-like invariant. It can be found using the same formula for Hamiltonian used to determine the energy function, except that instead of
t
there should be
x
.
G
=
∞
∫
0
g
d
t
,
d
G
d
x
=
-
u
t
u
x
∞
0
,
g
=
-
2
u
t
-
2
u
x
2
+
V
[
u
]
The physical meaning of
G
is not entirely clear to me. In addition, since it is not constant throughout the x-axis, it's hard to find a practical application for it.
Hyperbolic rotation symmetry
W
e
c
a
n
w
r
i
t
e
(
1
)
i
n
t
h
e
f
o
r
m
w
h
e
n
a
l
l
o
p
e
r
a
t
o
r
s
a
r
e
t
a
k
e
n
t
o
t
h
e
l
e
f
t
h
a
n
d
s
i
d
e
.
I
n
t
h
i
s
c
a
s
e
,
i
t
b
e
c
o
m
e
s
2
∂
t
-
2
∂
x
)
u
=
f
[
u
]
.
F
r
o
m
h
e
r
e
,
i
t
i
s
e
a
s
y
t
o
s
e
e
t
h
a
t
t
h
e
e
q
u
a
t
i
o
n
i
s
s
y
m
m
e
t
r
i
c
r
e
l
a
t
i
v
e
t
o
h
y
p
e
r
b
o
l
i
c
r
o
t
a
t
i
o
n
i
n
t
h
e
t
-
x
p
l
a
n
e
.
I
f
w
e
m
a
k
e
t
h
e
c
h
a
n
g
e
o
f
v
a
r
i
a
b
l
e
s
t
'
=
t
c
o
s
h
[
ε
]
+
x
s
i
n
h
[
ε
]
x
'
=
x
c
o
s
h
[
ε
]
+
t
s
i
n
h
[
ε
]
It can be shown that the equation (1), as well as
, are written in exactly the same way!
Energy flow
Let us derive an expression for the flow of energy, that I will denote
J
. In order to do that, let us use the law of conservation of energy.
∂
ρ
∂
t
+
∇
·
J
=
0
(
4
)
Since the problem is one-dimensional, del operator is identical to
∂
/
∂
x. Then, using the definition of
ρ
, we obtain
∂
∂
t
2
u
t
+
2
u
x
2
+
V
[
u
]
+
∂
J
∂
x
=
0
Which can be simplified to
u
t
u
t
t
+
u
x
u
x
t
+
V
'
[
u
]
u
t
=
-
∂
J
∂
x
Using the equations (1) and (3), we arrive at the following.
∂
J
∂
x
=
-
∂
∂
x
(
u
t
u
x
)
J
=
C
-
u
t
u
x
Where
C
is an arbitrary constant. From physical reasoning, the flow of energy must be zero if
u
and all its’ derivatives are zero. Then, the expression for the flow of energy is the following.
J
=
-
u
t
u
x
(
5
)
Stability of space-independent form of (1)
In this section, I'm studying the space-independent form of (1) to better understand the original equation. I've identified four types of behavior, one stable and two unstable.
Understanding the stability of (1) is an important stage in understanding the equation as a whole. But first, let us study the stability of space-independent form of (1), when
u[t, x +
Δ
x] = u[t, x]
for any
Δ
x
. As we will later see, the general form of (1) exhibits similar types of behavior to its' space-independent form.
I
f
u
[
t
,
x
+
Δ
x
]
=
u
[
t
,
x
]
,
t
h
e
n
,
u
s
i
n
g
t
h
e
d
e
f
i
n
i
t
i
o
n
o
f
d
e
r
i
v
a
t
i
v
e
u
x
=
l
i
m
h
0
u
[
t
,
x
+
h
]
-
u
[
t
,
h
]
h
i
t
'
s
e
a
s
y
t
o
s
h
o
w
t
h
a
t
u
x
=
u
x
x
=
0
f
o
r
a
l
l
t
a
n
d
x
.
A
l
s
o
s
i
n
c
e
u
n
o
l
o
n
g
e
r
d
e
p
e
n
d
s
o
n
x
i
t
c
a
n
b
e
w
r
i
t
t
e
n
a
s
u
=
u
[
t
]
.
S
o
,
t
h
e
e
q
u
a
t
i
o
n
(
1
)
t
h
e
n
c
a
n
b
e
w
r
i
t
t
e
n
a
s
2
d
u
2
d
t
=
f
[
u
]
(
6
)
This is the equation of motion of a unit mass experiencing a force
f
. It can exhibit four types of behavior: stationary, oscillatory, unbounded, and singular. Let us look at each of those in more detail.
Stationary
This type of behavior is probably the most boring of all. If at some time
t
0
d
u
d
t
[
t
0
]
=
0
f
[
u
[
t
0
]
]
=
0
2
d
u
2
d
t
[
t
0
]
=
0
then the value of
u
will always be equal to
u
[
t
0
]
.
The plot of such a solution looks just like a straight line.
P
l
o
t
[
1
,
{
t
,
0
,
1
}
]
O
u
t
[
]
=
Oscillatory
The equation (6) can be re-written as
d
1
2
2
d
u
d
t
=
-
d
V
Integrating this, we obtain
2
d
u
d
t
=
2
(
C
-
V
[
u
]
)
Where C is an arbitrary constant.
From this equation, it can be seen that the mass can only be moving if
V[u] < C
. If the mass is surrounded from both sides by potential barriers higher than
C
, it will never overcome any of them and instead will be oscillating somewhere in-between.
M
o
r
e
f
o
r
m
a
l
l
y
,
i
f
a
t
s
o
m
e
p
o
i
n
t
i
n
t
i
m
e
t
h
e
r
e
e
x
i
s
t
s
u
c
h
u
m
a
n
d
u
p
t
h
a
t
u
m
<
u
[
t
0
]
,
u
p
>
u
[
t
0
]
V
[
u
m
]
>
V
[
u
[
t
0
]
]
+
2
u
t
[
t
0
]
2
,
V
[
u
p
]
>
V
[
u
[
t
0
]
]
+
2
u
t
[
t
0
]
2
then
u
will forever remain bounded within that interval.
u
m
<
u
[
t
]
<
u
p
The plot of
u[t]
then looks something like this.
o
s
c
i
l
l
a
t
o
r
y
S
o
l
n
T
o
6
=
N
D
S
o
l
v
e
[
{
u
'
'
[
t
]
-
S
i
n
[
u
[
t
]
-
1
]
,
u
[
0
]
0
,
u
'
[
0
]
1
.
7
}
,
u
,
{
t
,
0
,
1
6
}
]
;
I
n
[
]
:
=
P
l
o
t
[
u
[
t
]
/
.
o
s
c
i
l
l
a
t
o
r
y
S
o
l
n
T
o
6
,
{
t
,
0
,
1
6
}
]
O
u
t
[
]
=
B
e
c
a
u
s
e
u
a
n
d
2
d
u
d
t
a
r
e
c
o
n
n
e
c
t
e
d
b
y
t
h
e
l
a
w
o
f
c
o
n
s
e
r
v
a
t
i
o
n
o
f
e
n
e
r
g
y
,
t
h
e
f
u
n
c
t
i
o
n
u
[
t
]
m
u
s
t
b
e
p
e
r
i
o
d
i
c
.
Unbounded
I
f
e
i
t
h
e
r
o
f
u
m
o
r
u
p
m
e
n
t
i
o
n
e
d
e
a
r
l
i
e
r
d
o
n
o
t
e
x
i
s
t
,
t
h
e
m
a
s
s
c
a
n
g
o
o
f
f
t
o
e
i
t
h
e
r
n
e
g
a
t
i
v
e
i
n
f
i
n
i
t
y
o
r
i
n
f
i
n
i
t
y
,
r
e
s
p
e
c
t
i
v
e
l
y
.
I
f
u
m
a
n
d
u
p
a
r
e
b
o
t
h
u
n
d
e
f
i
n
e
d
,
t
h
e
m
a
s
s
c
a
n
e
s
c
a
p
e
i
n
a
n
y
d
i
r
e
c
t
i
o
n
.
I
f
t
h
e
m
a
s
s
g
o
e
s
o
f
f
t
o
i
n
f
i
n
i
t
y
a
n
d
a
≤
l
i
m
u
∞
V
[
u
]
≤
b
w
h
e
r
e
a
a
n
d
b
a
r
e
s
o
m
e
n
u
m
b
e
r
s
,
t
h
e
n
d
u
d
t
w
i
l
l
b
e
b
o
u
n
d
e
d
.
I
n
o
t
h
e
r
w
o
r
d
s
,
t
h
e
v
a
l
u
e
o
f
u
i
n
c
r
e
a
s
e
s
i
n
d
e
f
i
n
i
t
e
l
y
b
u
t
d
o
e
s
s
o
a
t
a
r
e
s
t
r
i
c
t
e
d
r
a
t
e
.
T
h
e
p
l
o
t
o
f
a
n
u
n
b
o
u
n
d
e
d
s
o
l
u
t
i
o
n
l
o
o
k
s
s
o
m
e
w
h
a
t
l
i
k
e
t
h
i
s
.
u
n
b
o
u
n
d
e
d
S
o
l
n
T
o
6
=
N
D
S
o
l
v
e
[
{
u
'
'
[
t
]
-
S
i
n
[
u
[
t
]
-
1
]
,
u
[
0
]
0
,
u
'
[
0
]
1
.
7
8
}
,
u
,
{
t
,
0
,
1
6
}
]
;
I
n
[
]
:
=
P
l
o
t
[
u
[
t
]
/
.
u
n
b
o
u
n
d
e
d
S
o
l
n
T
o
6
,
{
t
,
0
,
1
6
}
]
O
u
t
[
]
=
Singular
If the point is not bounded by potential barriers and goes off, let’s say, to infinity, and
l
i
m
u
∞
V
[
u
]
=
-
∞
then the point will be moving infinitely fast. An example of a singular solution is plotted below.
s
i
n
g
u
l
a
r
S
o
l
n
T
o
6
=
N
D
S
o
l
v
e
[
{
u
'
'
[
t
]
u
[
t
]
,
u
[
0
]
0
,
u
'
[
0
]
1
.
7
8
}
,
u
,
{
t
,
0
,
1
6
}
]
;
I
n
[
]
:
=
P
l
o
t
[
u
[
t
]
/
.
s
i
n
g
u
l
a
r
S
o
l
n
T
o
6
,
{
t
,
0
,
1
6
}
]
O
u
t
[
]
=
Actually, there are two types of a singular solution: those that actually have a singularity and those that just rise quickly. To see why is that, let us use the following equation derived earlier.
2
d
u
d
t
=
2
(
C
-
V
[
u
]
)
I
f
u
g
o
e
s
o
f
f
t
o
i
n
f
i
n
i
t
y
,
t
h
e
v
a
l
u
e
o
f
d
u
d
t
w
i
l
l
a
l
w
a
y
s
r
e
m
a
i
n
p
o
s
i
t
i
v
e
.
T
h
e
n
,
d
u
d
t
=
√
(
2
(
C
-
V
[
u
]
)
)
From here,
d
t
d
u
=
1
√
(
2
(
C
-
V
[
u
]
)
)
And
t
[
u
]
=
1
√
2
u
∫
u
[
0
]
d
u
'
√
(
C
-
V
[
u
'
]
)
Let us find the time it takes for the solution to go to infinity. If the solution actually has a singularity,
u
will reach infinity in finite time. Else, the solution is simply a quickly increasing function.
t
∞
=
1
√
2
∞
∫
u
[
0
]
d
u
'
√
(
C
-
V
[
u
'
]
)
A
s
w
e
c
a
n
s
e
e
,
a
s
o
l
u
t
i
o
n
h
a
s
a
s
i
n
g
u
l
a
r
i
t
y
i
f
t
h
e
i
n
t
e
g
r
a
l
∞
∫
u
[
0
]
d
u
'
√
(
C
-
V
[
u
'
]
)
c
o
n
v
e
r
g
e
s
.
However, for our purposes quickly rising solutions and singular solutions are not too different from one another, so they both fell under one classification slot.
Stability of general form of (1)
When the derivative of u with respect to
x
is not zero, the equation (1) can no longer be reduced to (6). However, the stability classification system developed for (6) is still applicable. As earlier, there are four types of behavior, listed below.
Stationary
Here, everything is clear without words.
R
a
s
t
e
r
i
z
e
[
P
l
o
t
3
D
[
1
,
{
t
,
0
,
1
}
,
{
x
,
-
1
,
1
}
]
]
O
u
t
[
]
=
Oscillatory
In the case of (1), the criterion for oscillatory behavior is not as strict as that for (6). However, if the plot of
V[u]
looks like a pit, the behavior of (1) is likely to be oscillatory.
V
1
=
F
u
n
c
t
i
o
n
[
u
,
u
^
2
]
;
I
n
[
]
:
=
P
l
o
t
[
V
1
[
u
]
,
{
u
,
-
1
,
1
}
]
O
u
t
[
]
=
In the case of V[u] =
2
u
, for example, the solution looks like this.
m
y
P
l
o
t
[
s
o
l
v
e
T
h
e
P
D
E
[
f
V
t
o
F
[
V
V
1
]
]
]
O
u
t
[
]
=
Unbounded
Just like (6), the equation (1) can allow u to escape to infinity with a finite speed.
V
2
=
F
u
n
c
t
i
o
n
[
u
,
-
u
^
2
/
(
1
+
u
^
2
)
]
;
I
n
[
]
:
=
P
l
o
t
[
V
2
[
u
]
,
{
u
,
-
4
,
4
}
]
O
u
t
[
]
=
If the plot of V[u] looks somewhat like this, then the solution increases with a finite speed.
m
y
P
l
o
t
[
s
o
l
v
e
T
h
e
P
D
E
[
f
V
t
o
F
[
V
F
u
n
c
t
i
o
n
[
u
,
-
u
^
2
/
(
1
+
u
^
2
)
]
]
]
]
O
u
t
[
]
=
You can see that the solution growth with a steady rate bounded from above by the law of conservation of energy.
Singular
Some functions V[u] just yield singularities. If the value of u manages to escape to infinity and the limit of V[u] is negative infinity, u starts increasing infinitely fast.
V
3
=
F
u
n
c
t
i
o
n
[
u
,
-
u
^
2
]
;
I
n
[
]
:
=
P
l
o
t
[
V
3
[
u
]
,
{
u
,
-
3
,
3
}
]
O
u
t
[
]
=
s
i
n
g
u
l
a
r
S
o
l
n
=
s
o
l
v
e
T
h
e
P
D
E
[
f
V
t
o
F
[
V
V
3
]
]
;
I
n
[
]
:
=
m
y
P
l
o
t
[
s
i
n
g
u
l
a
r
S
o
l
n
]
O
u
t
[
]
=
An important thing about singular behavior is that if the u-function develops a singularity in one particular point, u also explodes in the neighboring points. To see why that's the case, let us plot
J =
-
u
t
u
x
at time 20.
J
f
o
r
S
i
n
g
u
l
a
r
S
o
l
n
[
x
_
]
:
=
-
F
i
r
s
t
[
D
e
r
i
v
a
t
i
v
e
[
1
,
0
]
[
u
]
[
2
0
,
x
]
*
D
e
r
i
v
a
t
i
v
e
[
0
,
1
]
[
u
]
[
2
0
,
x
]
/
.
s
i
n
g
u
l
a
r
S
o
l
n
]
;
I
n
[
]
:
=
P
l
o
t
[
J
f
o
r
S
i
n
g
u
l
a
r
S
o
l
n
[
x
]
,
{
x
,
-
R
,
R
}
,
P
l
o
t
L
e
g
e
n
d
s
{
"
J
"
}
,
P
l
o
t
R
a
n
g
e
A
l
l
,
P
l
o
t
S
t
y
l
e
O
r
a
n
g
e
]
J
O
u
t
[
]
=
As you can see, the flux of energy is positive to the right from the center and negative to the left. What this means is that to the right of the center the energy is flowing to the right, and to the left of the center it's flowing to the left, which basically implies that energy is being drained from the center.
V
e
c
t
o
r
P
l
o
t
[
{
J
f
o
r
S
i
n
g
u
l
a
r
S
o
l
n
[
x
]
,
0
}
,
{
x
,
-
2
0
,
2
0
}
,
{
y
,
-
1
,
1
}
,
V
e
c
t
o
r
P
o
i
n
t
s
T
a
b
l
e
[
{
x
,
0
}
,
{
x
,
-
7
,
7
,
1
}
]
,
A
s
p
e
c
t
R
a
t
i
o
1
/
6
,
I
m
a
g
e
S
i
z
e
L
a
r
g
e
]
O
u
t
[
]
=
This figure shows the direction of flux of energy around the center.
Upper and lower bounds for
u
The law of conservation of energy determines the maximal and minimal possible value
u
can have for a given energy
E
. In this subsection, I will prove the latter statement and find the bounds for
u
. For the sake of shortness, I will only find the upper bound and not the lower bound, as the process of their derivation is practically identical.
Why are there any bounds for
u
?
Due to the law of conservation of energy, the value of
u
is restricted to some finite interval. To see why that is, let us examine the expression for the energy of the system.
E
=
∞
∫
-
∞
2
u
t
+
2
u
x
2
+
V
[
u
]
d
x
O
n
e
o
b
v
i
o
u
s
w
a
y
t
o
m
a
k
e
u
l
a
r
g
e
i
s
t
o
m
a
k
e
i
t
l
a
r
g
e
a
t
o
n
e
p
o
i
n
t
w
h
i
l
e
m
a
k
i
n
g
i
t
z
e
r
o
e
v
e
r
y
w
h
e
r
e
e
l
s
e
.
H
o
w
e
v
e
r
,
a
s
i
t
'
s
e
a
s
y
t
o
s
h
o
w
,
t
h
a
t
w
o
u
l
d
m
a
k
e
t
h
e
i
n
t
e
g
r
a
l
∞
∫
-
∞
2
u
x
d
x
i
n
f
i
n
i
t
e
.
S
p
r
e
a
d
i
n
g
u
o
v
e
r
a
l
a
r
g
e
i
n
t
e
r
v
a
l
w
o
u
l
d
n
'
t
b
e
a
p
a
r
t
i
c
u
l
a
r
l
y
g
o
o
d
s
o
l
u
t
i
o
n
t
o
o
,
a
s
t
h
a
t
w
o
u
l
d
m
a
k
e
∞
∫
-
∞
V
[
u
]
d
x
l
a
r
g
e
.
T
h
e
n
,
w
e
c
o
n
c
l
u
d
e
t
h
a
t
t
h
e
o
p
t
i
m
a
l
a
r
r
a
n
g
e
m
e
n
t
l
i
e
s
s
o
m
e
w
h
e
r
e
i
n
-
b
e
t
w
e
e
n
t
h
o
s
e
t
w
o
.
T
h
e
f
i
g
u
r
e
b
e
l
o
w
s
h
o
w
s
t
h
e
q
u
a
l
i
t
a
t
i
v
e
p
l
o
t
s
o
f
t
h
e
t
h
r
e
e
s
o
l
u
t
i
o
n
s
:
t
o
o
s
h
a
r
p
,
t
o
o
s
m
o
o
t
h
,
a
n
d
o
p
t
i
m
a
l
.
P
l
o
t
[
{
I
f
[
A
b
s
[
x
]
<
0
.
1
,
0
.
5
,
0
]
,
0
.
7
5
E
x
p
[
-
x
^
2
/
8
]
,
E
x
p
[
-
A
b
s
[
x
]
]
}
,
{
x
,
-
3
,
3
}
,
P
l
o
t
S
t
y
l
e
T
h
i
c
k
,
A
x
e
s
O
r
i
g
i
n
{
-
2
,
0
}
,
P
l
o
t
L
e
g
e
n
d
s
{
"
T
o
o
s
h
a
r
p
"
,
"
T
o
o
w
i
d
e
"
,
"
O
p
t
i
m
a
l
"
}
,
E
x
c
l
u
s
i
o
n
s
N
o
n
e
]
T
o
o
s
h
a
r
p
T
o
o
w
i
d
e
O
p
t
i
m
a
l
O
u
t
[
]
=
Finding the bounds
Now, the question is how to find
an
optimal function
u[t, x]
. The word
an
here is important, as due to translation symmetry we can get as many optimal functions
u[t, x]
from one of them by shifting it either by either the
t
or
x
axis. To avoid ambiguity, let us restrict ourselves to such
u[t, x]
that have their global maximum at
t = 0, x = 0
.
Energy is the limiting factor in the process of increasing
u
. Thus, the optimal solution will be such that uses the minimal amount of energy.
F
r
o
m
t
h
e
e
x
p
r
e
s
s
i
o
n
f
o
r
e
n
e
r
g
y
,
i
t
i
s
o
b
v
i
o
u
s
t
h
a
t
a
t
t
=
0
,
u
t
m
u
s
t
b
e
z
e
r
o
e
v
e
r
y
w
h
e
r
e
:
s
i
n
c
e
u
t
o
n
l
y
e
n
t
e
r
s
t
h
e
i
n
t
e
g
r
a
l
u
n
d
e
r
t
h
e
s
q
u
a
r
e
s
i
g
n
,
m
a
k
i
n
g
u
t
a
n
y
t
h
i
n
g
o
t
h
e
r
t
h
a
n
z
e
r
o
w
o
u
l
d
i
n
c
r
e
a
s
e
t
h
e
e
n
e
r
g
y
w
h
i
l
e
n
o
t
i
n
c
r
e
a
s
i
n
g
t
h
e
m
a
x
i
m
a
l
v
a
l
u
e
o
f
u
.
L
a
t
e
r
i
n
t
h
i
s
s
e
c
t
i
o
n
w
e
w
i
l
l
p
r
a
c
t
i
c
a
l
l
y
f
o
r
g
e
t
a
b
o
u
t
t
h
e
t
i
m
e
d
e
p
e
n
d
e
n
c
y
o
f
u
a
n
d
,
f
o
r
c
o
n
v
e
n
i
e
n
c
e
,
c
a
l
l
u
[
x
]
=
u
[
0
,
x
]
.
W
i
t
h
u
a
n
d
u
x
t
h
i
n
g
s
a
r
e
m
o
r
e
c
o
m
p
l
i
c
a
t
e
d
:
w
e
n
e
e
d
t
o
f
i
n
d
a
b
a
l
a
n
c
e
b
e
t
w
e
e
n
m
a
k
i
n
g
t
h
e
s
o
l
u
t
i
o
n
t
o
o
s
h
a
r
p
a
n
d
g
i
v
i
n
g
t
o
o
m
u
c
h
e
n
e
r
g
y
t
o
t
h
e
V
[
u
]
t
e
r
m
.
F
o
r
t
h
a
t
,
w
e
c
a
n
u
s
e
t
h
e
E
u
l
e
r
-
L
a
g
r
a
n
g
e
e
q
u
a
t
i
o
n
.
S
i
n
c
e
E
=
∞
∫
-
∞
ρ
d
x
,
t
h
e
e
q
u
a
t
i
o
n
w
i
l
l
b
e
w
r
i
t
t
e
n
i
n
t
e
r
m
s
o
f
ρ
.
d
d
x
∂
ρ
∂
u
x
=
∂
ρ
∂
u
S
i
n
c
e
ρ
=
2
u
t
+
2
u
x
2
+
V
[
u
]
,
∂
ρ
∂
u
x
=
u
x
a
n
d
∂
ρ
∂
u
=
d
V
d
u
=
-
f
[
u
]
.
T
h
u
s
,
d
d
x
[
u
x
]
=
-
f
[
u
]
So, finally, we get an equation describing the optimal solution.
2
d
u
2
d
x
=
-
f
[
u
]
I
t
c
a
n
b
e
s
i
m
p
l
i
f
i
e
d
i
f
w
e
i
n
t
r
o
d
u
c
e
a
n
a
u
x
i
l
i
a
r
y
v
a
r
i
a
b
l
e
y
≡
d
u
d
x
.
T
H
e
l
a
t
t
e
r
c
a
n
t
h
e
n
b
e
r
e
-
w
r
i
t
t
e
n
a
s
f
o
l
l
o
w
s
.
y
d
y
d
u
=
-
f
[
u
]
Or
y
d
y
=
-
f
[
u
]
d
u
B
u
t
,
s
i
n
c
e
y
d
y
=
d
2
y
2
a
n
d
-
f
[
u
]
d
u
=
d
V
,
d
2
y
2
=
d
V
2
y
=
2
(
V
+
C
)
Where
C
is the constant of integration. Substituting the definition of
y
, we obtain
2
d
u
d
x
=
2
(
V
[
u
]
+
C
)
W
e
c
a
n
u
n
d
e
r
s
t
a
n
d
t
h
e
m
e
a
n
i
n
g
o
f
C
b
y
s
e
t
t
i
n
g
u
t
o
0
.
S
i
n
c
e
V
[
0
]
=
-
0
∫
0
f
[
u
]
d
u
=
0
,
C
=
1
2
2
d
u
d
x
u
=
0
.
H
e
r
e
2
d
u
d
x
u
=
0
m
e
a
n
s
t
h
e
s
q
u
a
r
e
o
f
d
u
d
x
a
t
t
h
e
p
o
i
n
t
w
h
e
n
u
=
0
.
S
o
,
C
b
a
s
i
c
a
l
l
y
m
e
a
s
u
r
e
s
h
o
w
s
t
e
e
p
t
h
e
p
l
o
t
o
f
u
i
s
w
h
e
n
i
t
c
r
o
s
s
e
s
t
h
e
u
=
0
s
u
r
f
a
c
e
.
B
y
t
h
e
i
n
i
t
i
a
l
c
o
n
d
i
t
i
o
n
s
,
t
h
e
v
a
l
u
e
o
f
u
a
n
d
u
x
a
t
i
n
f
i
n
i
t
y
e
q
u
a
l
s
z
e
r
o
.
O
n
e
w
a
y
t
o
s
a
t
i
s
f
y
t
h
e
s
e
c
o
n
d
i
t
i
o
n
s
i
s
t
o
s
e
t
C
t
o
z
e
r
o
.
T
h
i
s
,
h
o
w
e
v
e
r
,
i
s
n
o
t
a
l
w
a
y
s
v
i
a
b
l
e
,
a
s
t
h
e
s
h
a
p
e
o
f
t
h
e
s
o
l
u
t
i
o
n
w
o
u
l
d
t
h
e
n
b
e
g
o
v
e
r
n
e
d
b
y
t
h
e
e
q
u
a
t
i
o
n
2
d
u
d
x
=
2
V
[
u
]
.
T
h
i
s
e
q
u
a
t
i
o
n
o
n
l
y
a
l
l
o
w
s
t
h
e
v
a
l
u
e
o
f
u
t
o
c
h
a
n
g
e
w
i
t
h
x
i
f
V
[
u
]
>
0
.
W
h
a
t
t
h
i
s
m
e
a
n
s
i
s
t
h
a
t
t
h
e
v
a
l
u
e
o
f
u
c
a
n
n
o
t
s
u
r
p
a
s
s
p
o
i
n
t
s
w
h
e
r
e
V
[
u
]
i
s
z
e
r
o
.
B
u
t
s
i
n
c
e
w
e
a
l
s
o
w
a
n
t
t
o
c
o
m
p
u
t
e
a
n
o
p
t
i
m
a
l
s
o
l
u
t
i
o
n
f
o
r
a
r
b
i
t
r
a
r
y
p
e
a
k
v
a
l
u
e
o
f
u
,
w
e
n
e
e
d
t
o
f
i
n
d
a
d
i
f
f
e
r
e
n
t
a
p
p
r
o
a
c
h
.
One possible solution to the problem is to make
C
nonzero and forcefully equate
u
to zero at some two points and further away from the origin. This would satisfy the initial conditions and allow us to reach any values of
u
. Below you can see an illustration explaining how this would work. The solution is smooth everywhere except the points
x = 0
,
x =
π
/ 2
and
x = -
π
/ 2
.
P
l
o
t
[
U
n
i
t
S
t
e
p
[
P
i
/
2
-
A
b
s
[
x
]
]
(
P
i
/
8
+
C
o
s
[
x
]
-
A
b
s
[
x
]
/
4
)
,
{
x
,
-
3
,
3
}
,
P
l
o
t
S
t
y
l
e
T
h
i
c
k
]
O
u
t
[
]
=
In this case, the value of
C
should be chosen so that
u
would reach a certain value.
Anyway, the shape of this ideal solution is described by the equation
d
u
d
x
=
±
√
(
2
(
V
[
u
]
+
C
)
)
The sign of
u
x
depends on which side you are at: on the left it's positive and on the right it's negative. It is clear that the solution will have a mirror symmetry around the origin, so let us focus on the right half of the solution where
u
x
is negative.
d
u
d
x
=
-
√
(
2
(
V
[
u
]
+
C
)
)
S
i
n
c
e
d
u
d
x
=
1
d
x
d
u
,
w
e
c
a
n
f
i
n
d
x
a
s
a
f
u
n
c
t
i
o
n
o
f
u
b
y
f
l
i
p
p
i
n
g
t
h
e
l
a
t
t
e
r
e
x
p
r
e
s
s
i
o
n
a
n
d
i
n
t
e
g
r
a
t
i
n
g
.
x
[
u
]
=
u
∫
u
m
a
x
d
u
'
-
√
(
2
(
V
[
u
'
]
+
C
)
)
=
1
√
2
u
m
a
x
∫
u
d
u
'
√
(
V
[
u
'
]
+
C
)
After that, we can get
u[x]
as an inverse function of the latter integral. After that, energy as a function of maximal
u
can be found by integrating the energy density. Finally,
u
m
a
x
can be found as an inverse function of
E
[
u
m
a
x
]
.
Speed of signal propagation
For the vast majority of f[u], the speed at which waves travel depends on the shape of the waves themselves. In this section the latter is demonstrated on the example of linear
f[u]
. To prove the point, I perform mathematical manipulations and create some numerical solutions.
For the vast majority of f[u], the speed at which waves travel depends on the shape of the waves themselves. To see this, we can study the case of linear
f[u]
, when
f
[
u
]
=
-
c
u
And (1) takes the form
2
∂
u
∂
2
t
2
∂
u
∂
2
x
-
c
u
Since this is a linear equation, its' solution will be a sum of exponential functions of the form
u[t, x] = C Exp[i(
ω
t - k x)]
, where
i
is the imaginary unit. Solving the characteristic equation, we arrive at the following nonlinear relation between
ω
and k.
2
ω
=
2
k
+
c
A
s
i
t
'
s
k
n
o
w
n
,
n
o
n
l
i
n
e
a
r
d
i
s
p
e
r
s
i
o
n
r
e
l
a
t
i
o
n
s
c
a
u
s
e
d
i
f
f
e
r
e
n
t
t
y
p
e
s
o
f
w
a
v
e
s
t
o
t
r
a
v
e
l
w
i
t
h
d
i
f
f
e
r
e
n
t
s
p
e
e
d
s
.
B
e
l
o
w
y
o
u
c
a
n
s
e
e
t
w
o
p
l
o
t
s
,
b
o
t
h
s
h
o
w
i
n
g
n
u
m
e
r
i
c
a
l
s
o
l
u
t
i
o
n
s
t
o
t
h
e
s
a
m
e
P
D
E
2
∂
u
∂
2
t
2
∂
u
∂
2
x
-
2
4
u
o
n
t
h
e
s
a
m
e
i
n
t
e
r
v
a
l
o
f
v
a
l
u
e
s
o
f
t
a
n
d
x
.
T
h
e
h
o
r
i
z
o
n
t
a
l
a
x
i
s
i
s
x
(
f
r
o
m
l
e
f
t
t
o
r
i
g
h
t
)
,
a
n
d
t
h
e
v
e
r
t
i
c
a
l
a
x
i
s
i
s
t
(
f
r
o
m
b
o
t
t
o
m
t
o
t
o
p
)
.
T
h
e
o
n
l
y
d
i
f
f
e
r
e
n
c
e
b
e
t
w
e
e
n
t
h
e
t
w
o
p
l
o
t
s
a
r
e
t
h
e
i
n
i
t
i
a
l