search

IANW-Projects / ConservationLaws

A set of tools for numerically solving different kinds of conservation or balance laws with OpenCL.
3 stars 5 forks source link

A Dipole as a new example #31

Closed simonius closed 5 years ago

simonius commented 5 years ago

I wrote a new example. A Dipole stays on the left side of the YZ plane and the magnetic field magnitude oscillates on the opposite border ob the Domain. The simulation is semi-stable as it still works as of t = 2 but crashes at t = 10. I would like to hear your comments about it. I run the example with N = 100 in every direction and X_MIN = 0.0 = Y_MIN = Z_MIN, X_MAX = 1.0 = Y_MAX = Z_MAX.

Kostaszki commented 5 years ago

I'll be away on a business trip until wednesday, so I am not sure when I'll be able look at your code. I'll try to squeeze in some time on sunday evening.

ranocha commented 5 years ago

I'm also on a business trip until easter monday (with few days back at home in between). I'll try to have a look at this asap.

ranocha commented 5 years ago

You should rebase on master (or merge) in order to resolve conflicts.

simonius commented 5 years ago

I used I_TI('cfl') = 0.3 in my tests and the simulation was at least stable until t = 1. I also tested I_TI('cfl') = 0.65 and the simulation crashed. I concluded that the reason for the crash on your computer should be the higher CFL number.

ranocha commented 5 years ago

It seems that you merged something from my fork that is totally unrelated to this PR. In particular, the rebase on master and the addition of examples/ideal_MHD_Dipole.m is not included. Could you please revert that merge (and merge the correct branch, if you like).