Group Abstract Group Abstract

Message Boards Message Boards

0
|
4.2K Views
|
1 Reply
|
0 Total Likes
View groups...
Share
Share this post:

Divergence of numerical solution of discontinuous nonllinear ordinary ODE

POSTED BY: Zdenek Hurak

I think it helps to highlight the unitbox:

Show[StreamPlot[{v, 
   0.4` Sign[v] UnitBox[v, x] - 0.1` v - x}, {x, -0.6, 0.6}, {v, -0.6,
    0.6}, StreamStyle -> LightGray, 
  Prolog -> {LightBlue, Rectangle[{-0.5, -0.5}, {0.5, 0.5}]}], f3]

Blue Box

As you can see, the blue box is the source of the trouble. The curves are smooth until they hit the box. The part where it goes wrong is a very difficult part of the box. If you look at the stream, it is directing the curve back into the edge of the box.

enter image description here

It's getting pushed back into the discontinuity basically over and over. This part probably needs some very special treatment.

POSTED BY: Sean Clarke
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard